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ABSTRACT 


This document describes several navigation studies for low- 
altitude Earth satellites. The use of Global Positioning 
System Navigation Package data for Landsat-5 orbit determi- 
nation is evaluated. In addition, a navigation analysis for 
the proposed Tracking and Data Acquisition System is pre- 
sented. This analysis, based on simulations employing one- 
way Doppler data, is used to determine the agreement between 
the Research and Development Goddard Trajectory Determi- 
nation System and the Sequential Error Analysis Program re- 
sults. Properties of several geopotential error models are 
studied as they apply to orbit navigation error anlaysis, 
and an exploratory study of orbit smoother process noise is 
presented . 
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SECTION 1 


INTRODUCTION 


This document describes several studies performed for the 
National Aeronautics and Space Administration/Goddard Space 
Flight Center (NASA/GSFC) relating to low-altitude Earth 
satellite navigation analysis. The first two studies ana- 
lyzed the accuracy of satellite orbit determination using 
advanced navigation systems: the Global Positioning System 

(GPS) , currently under development, and the Tracking and 
Data Acquisition System (TDAS) , under study as a successor 
to the Tracking and Data Relay Satellite System (TDRSS ) . 

The latter two studies analyzed lumped geopotential modeling 
and orbit smoother process noise, two major components 
associated with orbit determination error analysis. 

The GPS navigation study analyzed the orbit determination 
capabilities of the Landsat-5 (705-kilometer (km) altitude, 
98-degree (deg) inclination) satellite. Landsat-5 uses 
navigation messages from GPS satellites that are extracted 
by an onboard Global Positioning System Navigation Package 
(GPSPAC) . This study extended the effort of a previous 
analysis^ - that considered Landsat-4 GPSPAC data. The re- 
sults of the Landsat-5 study and its comparison with the 
Landsat-4 study are presented in Section 2. 

The TDAS navigation study, presented in Section 3, analyzed 
the performance of a 600-km, 28-deg inclination satellite 
using simulated TDAS one-way Doppler data. The study used a 


1&. T. Fang and E. Seifert, "An Evaluation of Global Posi- 
tioning System Data for Landsat-4 Orbit Determination," 
paper delivered at AIAA 23rd Aerospace Science Meeting, 
January 1985; also Computer Sciences Corporation, 
CSC/TM-84/6077 , Tracking and Data Acquisition System/ 
Global Positioning System (TDAS/GPS) Navigation Analysis , 
September 1984. 
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number of TDAS enhancements that have been made to the Re- 
search and Development Goddard Trajectory Determination 
System (R&D GTDS) . It served as a verification of a previous 
study that considered the same satellite scenario but used 
the TDAS enhanced Sequential Error Analysis (SEA) program. 

The analysis of lumped geopotential modeling, presented in 
Section 4, studied several error models. All models use 
spherical-harmonic representations but with different sets 
of coefficients representing the errors. Global gravity 
error maps and orbit propagation errors corresponding to 
different error models are presented. 

The orbit ^smoother process noise study investigated the new 
smoother capacity to the SEA program. Limited investiga- 
tions using this capacity have been performed to further 
explore orbit smoothers and also to validate the fading mem- 
ory process noise option of the SEA smoother capability. 

The results of this study are presented in Section 5. 

Appendix A describes the validation of the SEA program 
smoother/TDAS capability used in the orbit smoother process 
noise study. Appendix B describes utilities used m con- 
junction with the lumped geopotential error model study. 
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SECTION 2 


LANDSAT-5 ORBIT DETERMINATION USING GPSPAC DATA 


This section discusses the accuracies of GPSPAC data for 
computing Landsat-5 orbits. The GPSPAC data were extracted 
from telemetry and processed by a large, sophisticated, 
batch orbit determination program (R&D GTDS) to produce 
16-hour arc orbit solutions. The accuracies of the GPSPAC 
data and the orbit solutions were inferred from the follow- 
ing: 

• Observation residuals 

• Overlap ephemeris comparisons 

• Comparisons with independent solutions derived from 
ground tracking data 

As with the Landsat-4 GPSPAC data, it was found that pseudo- 
range data appear inferior to delta pseudorange (Doppler) 
data. Efforts were made to determine whether pseudorange 
data are corrupted by several possible preprocessing sources 
of errors or whether data from individual GPS satellites may 
be defective. The latter effort involved computing 
Landsat-5 orbits from data from a single GPS satellite and 
was also a subject of interest by itself. 

Section 2.1 presents an overview of the data studied and the 
orbit determination scenarios used. Section 2.2 describes 
Landsat-5 orbit solutions computed from GPSPAC data; Sec- 
tion 2.3, solutions based on data originating from single 
GPS satellites; and Section 2.4, solutions obtained from 
ground tracking data and inquiries into their accuracy 
degradation in comparison with Landsat-4 solutions. Sec- 
tion 2.5 investigates several possible sources of preproc- 
essing errors and Section 2.6 summarizes the conclusions of 
the study. 
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2.1 OVERVIEW OF STUDY DATA 


The Landsat-4 spacecraft (705-km altitude, 98-deg inclina- 
tion orbit), launched in July 1982, carried an experimental 
NAVSTAR (Reference 2-1) GPSPAC and was the first satellite 
to use the GPS. The GPSPAC has a receiver/processor as- 
sembly that receives and decodes GPS tracking measurements 
and estimates the Landsat-4 position and velocity by using 
an onboard Kalman filter. The tracking measurements pro- 
viaea by GPSPAC are called pseudoranges and delta-pseudo 
ranges. The former are decoded from pseudorandom noise 
ranging code and represent measured signal transit times 
from the GPS satellites to Landsat-4. The latter are com- 
puted from Doppler shifts of carrier signals, integrated 
over a nominal interval of 0.6 second. The term "pseudo" is 
used because these measurements were derived from the Landsat 
clock, which was not synchronized with the GPS clocks. The 
quoted precisions of these data are 1.5 meters (m) for the 
pseudorange and 2 centimeters (cm) for the delta pseudorange. 

Landsat-5, launched in July 1984, is a replacement for the 
disabled. Landsat-4 spacecraft. It has similar mission ob- 
jectives and orbit characteristics. It also carries a 
GPSPAC with minor software modifications to avoid some known 
anomalies (Reference 2-2) that occurred in the Landsat-4 
GPSPAC. 

During the previous task assignment (42100) , an investiga- 
tion of the Landsat-4 GPSPAC data quality (Reference 2-3) 
was undertaken. The present study extended that effort to 
Landsat-5 GPSPAC data. As in the Landsat-4 study, the R&D 
GTDS batch orbit determination differential correction (DC) 
program was used to process both the GPSPAC telemetry to the 
ground and Ground Spaceflight Tracking and Data Network 
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(GSTDN) data. R&D GTDS allows extensive force modeling that 
is unavailable in the onboard GPSPAC Kalman filter. 

The data in this study were obtained over four separate time 
periods: m April and August of 1984, and in May and July of 
1985. The baseline force model and the orbit determination 
scenario used for the GSTDN and GPSPAC data solutions are as 
follows: 


GPSPAC data solutions 


21-by-21 Goddard Earth Model 
(GEM) -9 geopotential model 

Jacchia-Roberts atmospheric 
density model 

Cowell integrator 

16-hour data arc length 

Solve for orbital elements, 
drag parameter, clock bias, 
clock drift 

Select every 10th pair of 
observations 

Use observation standard de- 
viations = 1000 m (pseudo- 
range) and 0.8 cm (delta 
pseudorange) 


GSTDN data solutions 

21-by-21 GEM-9 geopotential 
model 

Jacchia-Roberts atmospheric 
density model 

Cowell integrator 

30-hour data arc length 

Solve for orbital elements, 
drag parameter 

Select all ab&ervations 

Use range-rate observations 
only for orbit solution 


2.2 ANALYSIS OF LANDSAT- 5 GPSPAC DATA ORBITS 

Eight 16-hour data arcs using the NAVSTAR system's GPS data, 
as received by Landsat-5 GPSPAC, were processed with the R&D 
GTDS DC program and the same baseline force model employed 
in analyzing the Landsat-4 GPSPAC orbits. (These orbits 
will henceforth be referred to as the baseline set of GPSPAC 
data orbits.) As will be discussed in Section 2.4, the 
Landsat-5 GSTDN solutions displayed ephemeris overlap dif- 
ferences consistently larger than those observed for 
Landsat-4. In comparison, the Landsat-4 GSTDN solutions 
displayed enough accuracy to be used as definitive solu- 
tions. Because of this degradation, the Landsat-5 GSTDN 
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orbits could not be used as definitive solutions for com- 
parison with corresponding GPSPAC solutions. Therefore, the 
present Landsat-5 analysis used the comparison among over- 
lapping GPSPAC data orbit ephemerides and observation re- 
siduals to provide an idea of the GPSPAC accuracies. 

The maximum ephemeris position differences for the four 
Landsat-5 GPSPAC data overlap periods vary from 64 to 19 m, 
compared to maxima with an average if 21 m, and a standard 
deviation of 14 m for six Landsat-4 GPSPAC overlap periods. 
As shown in Table 2-1, the 1984 GPSPAC overlap comparisons 
of 61 and 64 m are larger than the 1985 GPSPAC overlap com- 
parisons of 19 and 32 m. (It is interesting to note that 
studied GPSPAC data solutions in 1985 involve six GPS satel- 
lites whereas those in 1984 involve five; nonetheless, ref- 
erence to GPSPAC data solutions by the year of observation 
is for convenience only.) The improvement in the overlap 
comparisons for the 1985 GPSPAC data solutions, but no 
improvement in the corresponding 1985 GSTDN data (which is 
summarized in Section 2.4) , appears to rule out the possi- 
bility that some unexplained dynamic errors common to GSTDN 
and GPSPAC data solutions are responsible for poor 1984 
GPSPAC comparisons. 

An examination of data residuals of the Landsat-5 GPSPAC 
orbit solutions shows that large pseudorange measurement 
residuals exist and vary substantially from GPS satellite to 
GPS satellite. These characteristics are in agreement with 
those observed in the Landsat-4 GPSPAC data solutions. Fig- 
ure 2-1, which is a partial listing of a typical Landsat-5 
GPSPAC observation residual report illustrating the large 
pseudorange residuals, displays no correlation between 
pseudorange residuals and delta-pseudorange measurements. 
This rules out systematic time tag errors as possible 
sources for the systematic GPSPAC range errors observed. 
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Table 2-1. Landsat-5 GPSPAC Data Orbit Solution Characteristics and Comparisons 
With GSTDN Solutions 


RUN 

ID 

ARC START 
TIME 

ARC 

LENGTH 

NO OF 
OBSERVA 

SOLVE FOR PARAMETERS 

RESIDUAL 

STATISTICS 

STANDARD 

DEVIATIONS 

MAXIMUM POSITION DIFFERENCES 
OVER 4 HOUR OVERLAP OF ADJACENT 
ARCS (meters! 

MAXIMUM POSITION DIFFERENCES 
FROM GSTDN ORBITS (meters) 

REMARKS 

(date hour) 

(hours) 

TIONS 

Cl 

B 

(10 - 8 sec) 

• 

B 

(10" 8 sec/ 
sec) 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

RADIAL 

C 

ALONG 

TRACK 

RSS 

RADIAL 

C 

ALONG 

TRACK 

RSS 

G4A 

840413 

16 

1297 

1 592 

-3 127 

-3 262 

100 1 

2 39 





13 4 

306 

45 1 

543 



2 h 








569 

14 4 

61 0 

61 3 






G4B 

840413 
14 h 

16 

1192 

2 900 

-4 772 

-3 257 

71 50 

2 27 





10 6 

35 2 

70 3 

77 1 

ONLY EVERY 10TH 
PAIR OF OBSERVATION 
DATA PROCESSED 

G8A 

840807 

16 

1103 

2 125 

1 036 

-2 283 

84 25 

2 25 





14 0 

16 9 

49 1 

504 



3 h 








16 6 

76 

638 

639 






G8B 

840807 

16 

1110 

3 011 

-0 197 

-2 278 

83 02 

2 38 





24 3 

789 

680 

683 



15 h 

















G5A 

850524 

16 

1190 

4 913 

-2 228 

-0 6494 

76 44 

2 21 





13 5 

11 9 

52 2 

533 



19 h 








494 

7 67 

17 9 

19 4 






G5B 

850525 

16 

1146 

-6 119 

-2 508 

-0 6467 

66 33 

2 19 





18 9 

3 12 

77 2 

773 



7 h 

















G7A 

850726 

16 

1142 

2 493 

- 0 6720 

-0 3611 

84 39 

2 18 






_ 

_ 

_ 



0 h 








508 

19 4 

295 

32 1 






G7B 

850726 

16 

1000 

-2 259 

-0 8280 

-0 3571 

53 93 

2 18 





_ 

_ 

— 




12 h 


















NOTES SATELLITE AND TIME PERIOD LANDSAT 5 APRIL 1984 - JULY 1985 

DATA GPS RANGE AND DELTA RANGE WITH RANGE WEIGHTED LIGHTLY 
FORCE MODEL BASELINE 

SOLVE FOR PARAMETERS DRAG SCALE FACTOR C) = ACp/Cp, CLOCK BIAS B ANO DRIFT 8 AT EPOCH 
C L AND RSS ARE THE ACROSS TRACK ALONG TRACK AND ROOT SUM SQUARED COMPARE DIFFERENCES RESPECTIVELY 


0112 176a* )/85 




GTDS DC 


PROGRAM 


GTDS OBSERVATION RESIDUAL REPORT 


END OF ITERATION NUMBER 


M 

I 


TIME OF OBS 


GPS NO 

EDI T 

OBS 

O 


o-c 


RAT 10 

C 


ELV 

TO T 

A 


ob: 

Y YMMDD 

HHMM 

ss 

SSSS 



TYPE 





TO SIGMA 



GPS 



NO 

8404 14 

0335 

3 1 

329 1 

3 


RGPS 

20555 

9966 

57 

7255 

0 

0577 

20555 

9389 

26 

47 

0 

0 

1 099 

8404 1 4 

0335 

3 1 

329 1 

3 


DEL R 

0 

1936 

0 

47 29 

0 

59 1 1 

0 

1936 

26 

47 

0 

0 

1 100 

8404 1 4 

0336 

50 

6492 

3 


RGPS 

20608 

7939 

56 

0749 

0 

056 1 

20608 

7379 

25 

89 

0 

0 

110 1 

840414 

0336 

50 

6492 

3 


DELR 

0 

6073 

0 

5600 

0 

7000 

0 

6073 

25 

88 

0 

0 

1 102 

84 04 14 

0337 

56 

6295 

3 


RGPS 

20694 

0352 

56 

750 1 

0 

0568 

20693 

9784 

24 

99 

0 

0 

1 103 

840414 

0337 

56 

6295 

3 


DELR 

0 

9446 

-0 

8825 

1 

103 1 

0 

9447 

24 

98 

0 

0 

1 1 04 

8404 1 4 

0339 

1 1 

4 1 00 

3 


RGPS 

20834 

7875 

54 

9823 

0 

0550 

20834 

7325 

23 

54 

0 

0 

1 1 05 

8404 14 

0339 

1 1 

4 1 00 

3 


DELR 

1 

3 144 

2 

2863 

2 

8579 

I 

3 144 

23 

53 

0 

0 

1 1 06 

8404 1 4 

0340 

45 

4795 

2 


RGPS 

23799 

3095 

47 

80 1 9 

0 

0478 

23799 

26 1 7 

- 1 

3 1 

0 

0 

1 107 

8404 1 4 

0340 

45 

4795 

2 


DELR 

4 

3860 

- 1 

3533 

1 

69 16 

4 

3860 

- 1 

34 

0 

0 

1 108 

8404 1 4 

0342 

32 

4520 

3 


RGPS 

2 1428 

8755 

53 

50 1 7 

0 

0535 

2 1428 

8220 

1 7 

76 

0 

0 

1 109 

8404 1 4 

0342 

32 

4520 

3 


DELR 

2 

2053 

- 4 

56 15 

5 

70 18 

2 

2053 

1 7 

■>4 

0 

0 

1110 

8404 1 4 

0344 

02 

3532 

3 


RGPS 

2 1 784 

7896 

50 

3040 

0 

0503 

2 1 784 

7393 

1 4 

48 

0 

0 

till 

8404 14 

0344 

02 

3532 

3 


DELR 

2 

5404 

- 2 

4399 

3 

0499 

2 

5404 

1 4 

45 

0 

0 

1112 

8404 1 4 

0345 

24 

4544 

3 


RGPS 

22 150 

9298 

48 

98 1 6 

0 

0490 

22 150 

8808 

1 1 

2 1 

0 

0 

1113 

8404 1 4 

0345 

24 

4544 

3 


OELR 

2 

8066 

* 1 

5 1 24 

1 

0906 

2 

8066 

1 1 

18 

0 

0 

1114 

840414 

034 7 

27 

2582 

4 


RGPS 

23359 

5568 

83 

6347 

0 

0836 

23359 

47 3 1 

2 

07 

0 

0 

1115 

8404 1 4 

0347 

27 

2582 

4 


DELR 

-3 

0905 

3 

065 1 

3 

83 14 

-3 

0906 
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Figure 2-1. 


Partial Observation Residual Report for GPSPAC Data Orbit 
Solution (Run G4B) 
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Studies (discussed more fully in Section 2.5) that examine 
the effects of correcting assumed GPSPAC time tag errors 
give results that support this conclusion about no time tag 
errors. GPSPAC data observation residuals will be discussed 
further in Section 2.3. 

In keeping with the methods employed in the Landsat-4 GPSPAC 
analysis, the baseline GPSPAC data orbit determination 
scenario was repeated with two different a priori observa- 
tion weight specifications. When the weights were set to 
10 m for pseudorange observations and 1000 cm for delta- 
pseudorange observations (effectively producing a range-only 
solution) , the orbit comparison with the independent GSTDN 
data orbit produced larger differences than those seen in 
the baseline delta-range-only solutions in Table 2-1. When 
the weights were set to 1.5 m and 0.2 cm, respectively (in 
accordance with quoted data precision) , the orbit solutions 
resulted in larger overlap differences than the baseline 
delta-range-only solutions but were typically smaller than 
the corresponding range-only solution. These trends in the 
variously weighted Landsat-5 GPSPAC solutions are in agree- 
ment with the trends observed in the Landsat-4 data. 

Table 2-2 summarizes the Landsat-5 results. 

2.3 ORBIT SOLUTIONS USING DATA FROM SINGLE GPS SATELLITES 

As discussed in Section 2.2, the observed large pseudorange 
residuals present in the GPSPAC data solutions do not appear 
to be specific to any single GPS satellite. Nevertheless, 
there was interest in whether intercomparisons of GPSPAC 
solutions, each using data from only a single GPS NAVSTAR 
satellite, would divulge erratic data associated with a GPS 
or would provide insight into the cause of the observed large 
pseudorange residuals. It was uncertain whether single GPS 
relay solutions would even converge because of the lack of 
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Table 2-2. Landsat-5 GPSPAC Data Range-Rate and Evenly Weighted Orbit Solution 
Characteristics and Comparisons With GSTDN Solutions 


RUN 

ID 

ARC START 
TIME 

ARC 

LENGTH 

NO OF 
OBSERVA 

SOLVE FOR PARAMETERS 

RESIDUAL 

STATISTICS. 

STANDARD 

DEVIATIONS 

MAXIMUM POSITION DIFFERENCES 
OVER 4 HOUR OVERLAP OF ADJACENT 
ARCS (meters) 

MAXIMUM POSITION DIFFERENCES 
FROM GSTDN ORBITS (meters) 

REMARKS 

(date, hour) 

(hours) 

TIONS 

ei 

B 

<10— 3 sec) 

B 

(10 - ® sec/ 
sec) 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

RADIAL 

C 

ALONG 

TRACK 

RSS 

RADIAL 

C 

ALONG 

TRACK 

RSS 


R4A 

840413 

16 

1295 

-0 209 

-3 128 

-3 261 

41 85 

266 





25 

25 

70 

73 



2 h 








72 

18 

48 

51 





RANGE SOLUTION 

R4B 

840413 

16 

1191 

1 316 

-4 772 

-3 256 

50 35 

280 





30 

13 

112 

113 



14 h 

















R8A 

840807 

16 

1102 

19 757 

1 036 

-2283 

64 52 

3 23 





43 

15 

133 

134 


R8B 

3 h 

840807 

16 

1109 

1556 

-0 197 

-2 278 

68 52 

3 12 

45 

23 

171 

172 

31 

16 

99 

101 

ONLY EVERY 10TH 
PAIR OF OBSERVATION 
DATA PROCESSED 


15 h 

















P4A 

840413 

16 

1278 

-0 867 

-3 128 

-3 261 

36 79 

262 





24 

28 

64 

68 



2 h 








29 

14 

36 

39 






P4B 

840413 

16 

1183 

0 969 

-4 772 

-3 256 

48 01 

262 





26 

18 

101 

103 



14 h 
















EVENLY WEIGHTED 
SOLUTION 

P8A 

840807 

16 

665 

67 56 

1 035 

-2 283 

55 26 

308 





49 

33 

188 

191 


P8B 

3 h 

840807 
15 h 

16 

1103 

36 79 

-0 197 

-2 277 

63 93 

3 13 

48 

24 

177 

178 

32 

21 

122 

125 

ONLY EVERY 10TH 
PAIR OF OBSERVATION 
DATA PROCESSED 


NOTES SATELLITE AND TIME PERIOD LANDSAT 5 APRIL-AUGUST 1384 
FORCE MODEL BASELINE 

SOLVE FOR PARAMETERS DRAG SCALE FACTOR e, = ACp/Cp CLOCK BIAS B AND DRIFT B AT EPOCH 
C L AND RSS ARE THE ACROSS TRACK ALONG TRACK AND ROOT SUM SQUARED COMPARE DIFFERENCES RESPECTIVELY 
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sufficient geometry. First, all observations, except for 
those associated with the requested GPS satellite, were 
edited from the GPSPAC data. Next, the resulting single GPS 
satellite data were input into an unmodified R&D GTDS run to 
process GPSPAC data. These single relay orbit solutions-- 
except as otherwise specified — used the same dynamics, 
modeling, etc., employed in the corresponding composite re- 
lay solutions. 

Contrary to geometrical considerations, the single GPS 
satellite solutions displayed robustness in their ability to 
recover the clock parameters and to converge using very 
short observation arcs. The single GPS solutions can con- 
verge in as little as two Landsat orbital periods provided 
that the observation timespan's placement is chosen so as to 
minimize GPS-Landsat nonvisibility periods. Also, the 
single GPS satellite solutions, even those of short dura- 
tion, can recover the clock parameters (clock bias and 
drift) when no a priori information is provided. Table 2-3 
summarizes the characteristics of several single GPS satel- 
lite GPSPAC data solutions. 

Recovery of the clock parameters to a high degree of agree- 
ment, together with the single relay solution's ability to 
converge even with no a priori clock information, indicates 
that the clocks in the GPSPAC data are highly observable. 
This high clock observability provides the ability to re- 
solve small variations in the clock parameters. This clock 
variability is evident, in Table 2-3, between range and 
delta-range solutions over the same observation timespan. 

The clock parameters can be expressed in terms of their re- 
sulting corrections to pseudorange measurements to yield an 
effective range offset. Using this technique, the variabil- 
ity of the observed clock for the ensemble of single GPS 
satellite solutions over the same observation timespan is 
presented in Table 2-4. 
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Table 2-3. Characteristics of GPSPAC Data Solutions Using Observations From Single 
GPS Satellites 


RUN 3 

OBSERVATIONS 
FROM GPS 

ARC 

LENGTH 

(hours) 

NO OF 
OBSERVA- 
TIONS 

RESIDUAL 

STATISTICS, 

STANDARD 

DEVIATIONS 

SOLVE-FOR PARAMETERS 

COMMENTS ON RUN b 

ID 

SATELLITE 

NO 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

Cl 

B 

I10 -8 sec ) 

• 

B 

(10 -8 sec/ 
sec) 

G42A 

2 

3 

577 

2 35 

2 01 

-0 290 

r 3 1277 

-3 2631 

COMBINED SOLUTION 

G42B 

2 

3 

577 

2 35 

2 01 

-0 290 

-3 1277 

-3 2631 

COMBINED SOLUTION, ZERO A PRIORI 
CLOCK PARAMETERS 

G43A 

3 

2 5 

430 

41 86 

2 01 

0 001 

-3 1276 

-3 2640 

DELTA-RANGE SOLUTION 

G43B 

3 

2 5 

423 

1 53 

2 02 

-0 145 

-3 1280 

-3 2606 

COMBINED SOLUTION 

G43C 

3 

3 

592 

12 91 

2 80 

0 007 

-3 1269 

-3 2668 

RANGE SOLUTION 

G43D 

3 

3 

592 

24 34 

2 09 

-0 003 

-3 1279 

-3 2615 

DELTA-RANGE SOLUTION 

G43E 

3 

3 

585 

2 12 

2 11 

0 557 

-3 1276 

-3 2631 

COMBINED SOLUTION 

G43F 

3 

3 

585 

2 12 

2 11 

0 557 

-3 1276 

-3 2631 

COMBINED SOLUTION, ZERO A PRIORI 
CLOCK PARAMETERS 


NOTES SATELLITE AND TIME PERIOD LANDSAT-5, APRIL 1984 

DATA ALL AVAILABLE PAIRS OF OBSERVATIONS PROCESSED 
FORCE MODEL BASELINE 

SOLVE-FOR PARAMETERS DRAG SCALE FACTOR Cl = AC D /C D , CLOCK BIAS B, AND DRIFT B AT EPOCH 
a ALL SOLUTION ARC TIMES START ON 840413, 2 HOURS 

b ALL SOLUTIONS INVOLVE PAIRS OF GPS RANGE AND DELTA-RANGE OBSERVATIONS IN COMBINED SOLUTIONS, BOTH DATA TYPES HAVE COMPARABLE 
WEIGHTINGS IN DELTA-RANGE SOLUTIONS, RANGE IS WEIGHTED LIGHTLY, IN RANGE SOLUTIONS, DELTA-RANGE IS WEIGHTED LIGHTLY 
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Table 2-4. Intercomparisons of GPSPAC Data Solutions Using Single and Composite 
GPS Observation Sets 


RUN 3 

ID 

OBSERVATIONS 
FROM GPS 
SATELLITE 
NO 

NO OF 
OBSERVA- 
TIONS 13 

RESIDUAL 

STATISTICS, 

STANDARD 

DEVIATIONS 

SOLVE-FOR PARAMETERS 

RANGE EQUIVALENT 
CLOCK DIFFERENCE 
FROM COMPOSITE 0 
(meters) 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

Cl 

B 

(10 - 3 sec ) 

• 

B 

(10 - ® sec/ 
sec) 

G4A 

ALL GPSs— COMPOSITE 

1297 

100 1 

2 39 

1 591 

-3 12752 

-3 26234 

- 

G41 

1 

2650 

67 0 

2 20 • 

-3 280 

-3 12855 

-3 25963 

76 

G42 

2 

2446 

152 6 

2 20 

0 487 

-3 12708 

-3 26375 

-10 

G43 

3 

2581 

185 3 

2 27 

1 621 

-3 12681 

-3 26434 

-38 

G44 

4 

3000 

180 4 

2 49 

9 872 

-3 12924 

-3 25785 

130 

G45 

5 

2080 

55 4 

2 29 

-0 331 

-3 12828 

-3 26020 

45 


NOTES SATELLITE AND TIME PERIOD LANDSAT-5, APRIL 1984 


DATA GPS RANGE AND DELTA-RANGE, WITH RANGE WEIGHTED LIGHTLY 
FORCE MODEL BASELINE 

SOLVE-FOR PARAMETERS DRAG SCALE FACTOR ci = AC D /C D , CLOCK BIAS B, AND DRIFT B AT EPOCH 
a ALL SOLUTIONS HAVE DATA ARC LENGTHS OF 16 HOURS STARTING ON 840413, 2 HOURS 

b FOR COMPOSITE SOLUTION, ONLY EVERY 10TH PAIR OF OBSERVATION DATA IS PROCESSED, FOR SINGLE GPS SATELLITE SOLUTIONS, 
ALL PAIRS OF OBSERVATION DATA ARE PROCESSED 

C THIS VALUE IS THE DIFFERENCE BETWEEN THE RANGE EQUIVALENT CLOCK ERROR FOR RUN ID G4A AND THE RANGE EQUIVALENT 
CLOCK ERROR FOR THE GPS SATELLITE SOLUTION CORRESPONDING TO THE TABLE ENTRY THE RANGE EQUIVALENT OF CLOCK ERROR 
IS CALCULATED AS dR = (B + B At) C, WHERE At IS CHOSEN AS HALF THE SOLUTION ARC LENGTH, OR 8 HOURS, AND C IS THE SPEED 
OF LIGHT 
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Possibly related to the observed clock variations are char- 
acteristic gross structure signatures of the pseudorange re- 
sidual in the GPSPAC solution plots. The typical pseudorange 
residual plots can be fit fairly well with a linear, nonzero 
slope function. This is most pronounced in single GPS 
satellite solutions but is also present in composite GPS 
solutions. Representative pseudorange residual plots for 
both single GPS satellite and composite GPS solutions are 
shown in Figure 2-2. One possible interpretation of the 
nonzero slope structure in the pseudorange residuals is that 
an inconsistency occurs in the observed clock drift rate 
between the pseudorange and delta-pseudorange GPSPAC data 
measurements. 

In the search for information pertaining to the anomalous 
large pseudorange residuals, the inconsistency between 
pseudorange and delta-pseudorange measurements was redis- 
covered. This was previously discussed in connection with 
comparisons of GPSPAC data delta-range solutions to range 
solutions. In those instances, the inconsistency was most 
notable in the Landsat-4 analysis, when the delta-range 
solutions gave better comparison with GSTDN solutions than 
the corresponding range or evenly weighted solutions. In 
its present manifestation, the inconsistency can be inter- 
preted as a mismatch between the observed clock parameters 
for delta pseudorange and pseudorange data. Alternately, 
the inconsistencies could be attributed to uncorrected 
signal propagation errors or difficulties with electronics. 

Finally, intercomparisons between single GPS satellite solu- 
tions and between single satellite and composite solutions 
did not reveal any highly erratic single GPS satellite. All 
single satellite solutions agree within a maximum of ^ 80 m, 
and typically less than ^ 40 m. Vvhen observations corre- 
sponding to that GPS satellite for which single GPS satel- 
lite solution intercomparisons show the largest deviations 
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MINUTES AFTER MM 13 0 MOORS 


Figure 2-2. Representative Observation Residual Plots for GPSPAC Data Solutions 
Employing Composite and Single GPS Satellite Observation Sets 

(a) Range residuals for a solution that employs observations from all available GPS 
NAVSTAR satellites - a composite solution (Run G7B) . 

(b) Composite solution's corresponding delta-range residuals. Delta- range residual 
plots for single GPS satellite solutions display similar signature and residuals. 

(c) Range residuals for a GPSPAC data solution using observations from a single GPS 
satellite (Run G-43) . 





were edited from a composite solution, the large pseudorange 
residuals are still seen. 

2.4 ASSESSMENT OF LANDSAT- 5 ORBITS COMPUTED FROM GSTDN DATA 

2.4.1 COMPARISON OF LARGE OVERLAP DIFFERENCES TO LANDSAT-4 
GSTDN SOLUTIONS 

GSTDN data for this task were processed in 30-hour overlap- 
ping arcs, with overlap periods of 6 hours each. The dif- 
ference between two adjacent orbital solutions in the 
overlap period serves as a measure of the GSTDN data- 
determined orbital accuracy. In this study, Landsat-5 GSTDN 
data from two different calendar years were processed. 
Overlapping arcs, sets of two runs _each, were processed for 
April and August 1984 and for May 1985. The first two sets 
correspond to available GPSPAC data sets. The May 1985 
GSTDN data were chosen, based on optimal GSTDN and GPSPAC 
data availability, from 13 timespans in 1985. 

The maximum ephemeris position differences for the three 
overlapping periods were 106, 59, and 110 m. The root- 
mean-square (rms) ephemeris position differences for the two 
overlap periods were 64, 33, and 83 m. These differences 
are considerably greater than those in the Landsat-4 study. 
For Landsat-4, the maximum ephemeris position differences of 
seven overlap comparison periods averaged 29 m, with a 
standard deviation of 13 m. The rms ephemeris position dif- 
ferences for the Landsat-4 overlap periods averaged 18 m, 
with a standard deviation of 7 m. Tables 2-5 and 2-6 sum- 
marize the results of Landsat-5 GSTDN characteristics. 

Although this study processed only a limited number of GSTDN 
tracking intervals, independent results indicate that 1984 
and 1985 Landsat-5 solutions give consistently larger over- 
lap differences when compared to the previous Landsat-4 
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Table 2-5 


1984 GSTDN Data Orbit Solution Characteristics 







COMPUTED 
DRAG SCALE 
FACTOR 

RESIDUAL 

MAXIMUM POSITION DIFFERENCES 


RUN 

ID 

ARC START 
TIME 

ARC 

LENGTH 

NO OF 
TRACKING 
STATIONS 
INVOLVED 

NO OF 
OBSERVA- 

STATISTICS, 

RANGE-RATE 

STANDARD 

DEVIATIONS 

OVER 6-HOUR OVERLAP OF ADJACENT 
ARCS (meters) 

REMARKS 

(date, hour) 

(hours) 

TIONS 

Cl = ac d/Cd 

RADIAL 

CROSS- 

ALONG- 

RSS 








(cm/sec) 

TRACK 

TRACK 

TOTAL 


GDA 

840413 

30 

3 

507 

-1 472 

3 048 





BASELINE MODEL JACCHIA 


0 h 






10 6 

13 8 

105 7 

105 8 

ROBERTS ATMOSPHERIC 
DENSITY MODEL, A PRIORI 

GDB 

840414 

30 

5 

702 

-0 919 

4 114 





DRAG COEFFICIENT C D = 


0 h 



399 







2 0, A DRAG SCALE FAC- 
TOR SOLVED FOR, 21 x 21 
GEM 9 GEOPOTENTIAL 

GDM 

840807 

30 

3 

-8 833 

3 478 





MODEL 


0 h 






8 1 

87 

584 

59 2 


GDN 

840808 

30 

4 

444 

-7 493 

2 956 







0 h 











S4A 

840413 

30 

3 

507 

-1 965 

3 035 





HARRIS-PRIESTER F100 


0 h 






10 2 

13 8 

101 1 

101 3 

ATMOSPHERIC DENSITY 
MODEL, A PRIORI DRAG 

S4B 

840414 

30 

5 

702 

-0 674 

4 100 





COEFFICIENT C D = 2 0, A 


0 h 



399 

-3 858 

3 454 





DRAG SCALE FACTOR 
SOLVED FOR, 21 x 21 GEM 
9 GEOPOTENTIAL MODEL 

S8A 

840807 

30 

3 






0 h 






77 

8 1 

69 8 

70 3 


S8B 

840808 

30 

4 

444 

-4 348 

2 960 







0 h 











GMA 

840413 

30 

3 

507 

-1 549 

3 238 





BASELINE ATMOSPHERE 


0 h 






10 7 

13 6 

101 2 

101 7 

MODEL, WITH 36 x 36 GEM 
10B GEOPOTENTIAL MODEL 

GMB 

840414 

30 

5 

702 

-1 092 

3 860 







0 h 





* 






GMC 

840413 

30 

3 

507 

-1 136 

2 775 





BASELINE MODEL WITH 36 


0 h 






95 

98 

83 9 

842 

x 36 GEM 10B GEOPOTEN- 
TIAL MODEL, MARSH STA- 

GMD 

840414 

30 

5 

703 

-0 959 

2 822 





TION COORDINATES 


0 h 












NOTES SATELLITE AND TIME PERIOD LANDSAT-5, APRIL-AUGUST 1984 
DATA GSTDN RANGE-RATE DATA 
RSS, ROOT-SUM-SQUARED 
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Table 2-6 


1985 GSTDN Data Orbit Solution Characteristics 








RESIDUAL 

STATISTICS, 

STANDARD 

DEVIATIONS 

MAXIMUM POSITION DIFFERENCES 


RUN 

ID 

ARC START 
TIME 

ARC 

LENGTH 

NO OF 
TRACKING 
STATIONS 
INVOLVED 

NO OF 
OBSERVA- 

COMPUTED 
DRAG SCALE 
FACTOR 
el = a Cq/Cq 

OVER 6-HOUR OVERLAP OF 
ADJACENT ARCS (meters) 

REMARKS 

(date, hour) 

(hours) 

TIONS 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

RADIAL 

CROSS- 

TRACK 

ALONG- 

TRACK 

RSS 

TOTAL 


GS1 

850524 

30 

3 

239 

-11 07 

— 

1 319 





BASELINE MODEL 


6 h 







11 1 

26 4 

107 3 

1102 


GS2 

850525 

30 

5 

270 

440 

- 

1 396 





RANGE-RATE DATA ONLY 


6 h 












GSA 

850524 

30 

4 

400 

-8 97 

10 47 

5 361 





NO OBSERVATION PROCESS- 


6 h 







16 2 

19 7 

69 0 

69 1 

ING PERFORMED, BOTH 
RANGE AND RANGE-RATE 

GSB 

850525 

30 

6 

387 

5 73 

4 322 

5 775 





DATA INCLUDED 


6 h 













NOTES SATELLITE AND TIME PERIOD LANDSAT-5, MAY 1985 
DATA GSTDN RANGE-RATE DATA 
RSS IS ROOT-SUM-SQUARED 
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GSTDN solutions. Tracking statistics do indicate an in- 
creasing amount of TDRS-East tracking data available for 
Landsat-5 in 1985. Overlap differences of TDRS-only solu- 
tions, associated with the single-TDRS S-band certification 
results, report 82, 71, 35, and 74 m (Reference 2-4). In 
addition, private conversations with investigators 
(Yuri Nakai, Computer Sciences Corporation) who are analyz- 
ing 1985 Landsat-5 orbit solution overlap differences using 
combinations of GSTDN, TDRSS, and Bilateration Ranging 
Transponder System (BRTS) data indicate ephemeris overlap 
differences consistently larger than the average 29 m ob- 
served for Landsat-4. 

2.4.2 ANALYSIS OF LANDSAT-5 LARGER-THAN-EXPECTED OVERLAP 
EPHEMERIS DIFFERENCES 

It was originally thought that these larger-than-expected 
overlap ephemeris differences could be attributed to inade- 
quacies in atmospheric density modeling. Modeling error was 
suspected because the overlap comparison plot of the along- 
track difference, as shown in Figure 2-3, displays sinu- 
soidally increasing differences that are generally typical 
of force modeling mismatch between the orbit solutions in- 
volved. No significant improvement was seen, however, when 
a different atmospheric model (Harris-Pr iester ) was used. 

In an attempt to reduce these ephemeris differences, the 
April 1984 GSTDN data arcs were again processed using the 
standard orbit determination scenario, but with two varia- 
tions. First, orbits were obtained by using the 36-by-36 
GEM-10B geopotential model. Second, a second set of orbits 
was obtained by using the Marsh tracking station position 
coordinates in addition to the GEM-10B model. As was also 
true in the Landsat-4 study, using these improved models did 
not result m substantial improvement in orbit determination 
performance. The results of the 1984 GSTDN orbit determina- 
tion runs are summarized m Table 2-5. 
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Figure 2-3. Landsat-5 GSTDN Orbit Solution Overlap Ephemens Comparison 



Although the previous Landsat-4 experience with unified 
S-band (USB) tracking data revealed systematic errors in the 
range observations/ the May 1985 GSTDN data arcs were re- 
processed using both range and range-rate observations. 

This was done to determine if range errors were caused by 
redundant observation correction preprocessing in R&D GTDS . 
When the preprocessing was deleted, the resulting orbit 
solutions still displayed large range and range-rate resid- 
ual statistics. Probably because of the greater observabil- 
ity afforded by the range observations, the maximum position 
difference was decreased from 110 to 69 m--still large in 
comparison with the Landsat-4 accuracies. The results of 
the 1985 GSTDN orbit determination runs are summarized in 
Table 2-6. 

2.4.3 LANDSAT-5 SOLUTIONS INVOLVING FEWER TRACKING STATIONS 

In comparison to previous Landsat-4 GSTDN data, the 
Landsat-5 solutions involved fewer ground tracking sta- 
tions. In the typical Landsat-4 GSTDN 30-hour arc solu- 
tions, between 6 and 10 ground stations were involved in 
tracking. For the six baseline Landsat-5 GSTDN solutions in 
Tables 2-5 and 2-6, between three and five tracking stations 
were involved. 

Irregularity of the Lanasat-5 tracking passes may also con- 
tribute to the larger-than-expected overlap ephemeris dif- 
ferences. Figure 2-4 displays a plot of the ground tracking 
passes for Landsat-5 GSTDN data, beginning on April 13, 

1984. This plot shows the 30-hour intervals for baseline 
run IDs GDA and GDB and their corresponding overlap interval 
extending from 0 hours to 6 hours on April 14 (840414). 

Because it was suspected that the sparse tracking over the 
GDA-GDB 6-hour overlap contributed to its poor ephemeris 
overlap comparison (from Table 2-5, a maximum position dif- 
ference rss of 106 m) , alternate 30-hour arcs were chosen. 
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NOTES CAPITAL LETTERS INDICATE THE CORRESPONDENCE BETWEEN TRACKING 
PASSES (WITH A TYPICAL DURATION OF 10 MINUTES) AND THE TRACKING 
STATION FROM WHICH THE MEASUREMENT SIGNAL ORIGINATED 

THE TRACKING PASSES ABOVE THE TIME AXIS CONTRIBUTE TO THE OR- 
BIT SOLUTION OVER THE DATA ARC 

THE TRACKING PASSES BELOW THE TIME AXIS ARE PRESENT IN THE 
GSTDN DATA BASE 


Figure 2-4 


Plot ot Ground Tracking Passes and Orbit Solution Arcs Involving 
April 1984 Landsat-5 GSTDN Data 
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These alternate arcs, labeled GDC and GDD, begin at 840413 
18 hours and 840414 18 hours, respectively, and have a 
6-hour overlap beginning at 840414 18 hours. As seen from 
the plot in Figure 2-4, this alternate overlap time includes 
more uniform tracking. This is reflected in the GDC-GDD 
overlap ephemeris comparison RSS of 44 m (from Table 2-7) , 
which is a significant improvement over the GDA-GDB overlap 
(although it is still inferior to the Landsat-4 GSTDN over- 
lap maxima, which have an RSS average of 28.8 m) . Part of 
this improvement may also result because the alternate GDC 
and GDD arcs involve more tracking stations, four and six, 
respectively, compared to three and five for GDA and GDB . 

A second set of GSTDN data orbit solutions were performed, 
using 39-hour arcs. These two arcs, labeled GDE and GDF, 
also have a 6-hour overlap, beginning at 840414 9 hours. 

The maximum position difference rss of this overlap is 72 m 
(Table 2-7) . Although the number of tracking stations in- 
volved remained the same as in the GDC-GDD solutions, the 
larger orbit arcs resulted in a poorer overlap comparison. 
This suggests that the relatively small number of stations 
involved in Landsat-5 GSTDN data cannot be compensated for 
by extending orbit solutions over arcs of greater duration, 
thereby increasing the number of observations and possibly 
the number of tracking stations. Instead, the GDE-GDF over- 
lap, when compared to GDC-GDD errors, indicates that the 
force modeling errors begin to increase at a greater rate 
than the observational benefits. 

2.4.4 ANALYSIS OF LANDSAT-4 GSTDN SOLUTIONS WITH REDUCED 
NUMBER OF TRACKING STATIONS 

To determine if poor Landsat-5 GSTDN solutions result from 
the small number of tracking stations involved, Landsat-4 
GSTDN solutions with a reduced number of tracking stations 
were generated. The results of these solutions are sum- 
marized in Table 2-8. The first entry m the table, the 
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Table 2-7. Additional Landsat-5 GSTDN Data Orbit Solution Characteristics 


RUN 

ID 

ARC START 
TIME 

(date, hour) 

ARC 

LENGTH 

(hours) 

NO OF 
TRACKING 
STATIONS 
INVOLVED 

NO OF 
OBSERVA- 
TIONS 

COMPUTED 
DRAG SCALE 
FACTOR 
Cl = AC D /Cd 

RESIDUAL 

STATISTICS, 

RANGE-RATE 

STANDARD 

DEVIATIONS 

Icm/sec) 

MAXIMUM POSITION DIFFERENCES 
OVER 6-HOUR OVERLAP OF ADJACENT 
ARCS (meters) 

REMARKS 

RADIAL 

CROSS- 

TRACK 

ALONG- 

TRACK 

RSS 

TOTAL 

GDC 

840413 

30 

4 

703 

2 045 

384 





BASELINE MODEL 


18 h 






84 

10 4 

442 

44 3 


GDD 

840414 

30 

6 

723 

0 065 

354 







18 h 











GDE 

840413 

39 

4 

767 

-0 424 

444 







0 h 






14 5 

13 5 

71 4 

71 6 


GDF 

840414 

39 

6 

829 

-0 390 

4 31 







9 h 












NOTES SATELLITE AND TIME PERIOD LANDSAT-5, APRIL 1984 
DATA GSTDN RANGE-RATE DATA 
RSS IS ROOT SUM SQUARED 

SEE FIGURE 2-4 FOR TRACKING PASSES AND ORBIT SOLUTION ARCS 
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Table 2-8 


Summary of Landsat-4 GSTDN Solutions With Peduced Number of Tracking Stations 


RUN 

ID 

ARC START 
TIME 

NO OF 
TRACKING 
STATIONS 
ACCEPTED 

NO OF 
TRACKING 
PASSES 
ACCEPTED 

NO OF 
OBSERVA- 

COMPUTED 
DRAG SCALE 
FACTOR 

RESIDUAL 

STATISTICS, 

RANGE-RATE 

STANDARD 

DEVIATIONS 

(cm/sec) 

MAXIMUM POSITION DIFFERENCES 
OVER 6-HOUR OVERLAP OF ADJACENT 
ARCS (meters) 

REMARKS 

(date, hour) 

TIONS 

ei = AC D /c d 

RADIAL 

CROSS 

TRACK 

ALONG- 

TRACK 

RSS 

TOTAL 

GL4 

821001 
0 h 

10 

16 

145 

-0 236 

2 62 

24 

152 

62 

15 5 

LANDSAT-4 BASELINE 
SOLUTION (PRESENTED 
FOR COMPARISON) 

GL5 

821002 

6 

16 

116 

-0 294 

3 45 







0 h 











GL41 

821001 
0 h 

3 

4 

207 

-0 779 

3 17 

10 3 

23 1 

155 2 

155 4 

LANDSAT-4 SOLUTIONS 
WITH REDUCED NO OF 
TRACKING STATIONS 

GL51 

821002 

4 

9 

468 

-0 347 

294 







0 h 






6 5 

82 

25 6 

25 7 


GL42 

821001 

4 

9 

453 

-0 251 

2 89 





* 


0 h 











GL43 

821001 

4 

6 

304 

-0 286 

2 99 







0 h 






11 7 

13 7 

30 5 

33 3 


GL52 

821002 

2 

7 

367 

-0 349 

2 87 







0 h 






7 0 

54 

28 7 

28 7 


GL44 

821001 

6 

11 

661 

-0 268 

2 60 







0 h 












NOTES SATELLITE AND TIME PERIOD LANDSAT-4, OCTOBER 1982 
DATA GSTDN RANGE-RATE DATA 
RSS IS ROOT SUM SQUARED 

ALL SOLUTIONS HAVE DATA ARC LENGTHS OF 30 HOURS 

SEE FIGURE 2-5 FOR TRACKING PASSES AND ORBIT SOLUTION ARCS 
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overlap comparison between the nonreduced Landsat-4 arcs GL4 
and GL5, is presented as the standard for comparison. The 
remaining entries in the table present comparative ephemeris 
overlap differences for orbits generated over the same time- 
spans as GL4 and GL5 but using different reduced sets of 
tracking stations. 

The Landsat-4 GL4-GL5 set was chosen because it displays the 
smallest ephemeris overlap difference of all the Landsat-4 
GSTDN comparisons. It was reasoned that the GL4-GL5 set 
would therefore produce the greatest range of sensitivity to 
reduced tracking station scenarios. In addition, because of 
the large number of stations involved, it was hoped, that 
appropriate subsets of the GL4-GL5 tracking stations could 
be found to approximately reproduce the tracking geometries 
available in the Landsat-5 GSTDN solutions. Unfortunately, 
the latter purpose was frustrated by the absence of valid 
Landsat-4 observations corresponding to tracking stations in 
the Landsat-5 station set and the geometrical complications 
introduced in choosing substitution stations. 

The reduced station results in Table 2-8 do produce poorer 
overlap comparisons. Specifically, as seen in set GL41-GL51 
in comparison to set GL42-GL51, overlap comparisons dete- 
riorate rapidly after some cutoff number of tracking sta- 
tion. As shown in the table, below this cutoff number the 
number of tracking passes involved begins to play an in- 
creasing importance in the overlap differences. As shown in 
Figure 2-5, which displays the tracking schedules for those 
runs summarized in Table 2-8, the number and duration of 
nontracking intervals also appears to affect the ephemeris 
differences. 

Because of these complications introduced by station geom- 
etry, the number of tracking passes, and the amount and 
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Figure 2-5. Plot of Ground Tracking Passes and Orbit Solution Arcs Involving 
October 1982 Landsat-4 GSTDN Data 


NOTES CAPITAL LETTERS INDICATE THE CORESPONDENCE BETWEEN TRACKING 
PASSES (WITH A TYPICAL DURATION OF 10 MINUTES) AND THE TRACK- 
ING STATION FROM WHICH THE MEASUREMENT SIGNAL ORIGINATED 

THE TRACKING PASSES ABOVE THE TIME AXIS CONTRIBUTE TO THE OR- 
BIT SOLUTION OVER THE DATA ARC 

THE TRACKING PASSES BELOW THE TIME AXIS ARE PRESENT IN THE 
GSTDN DATA BASE 
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duration of nontracking intervals, the reduced station re- 
sults were inconclusive in determining the extent to which 
poor Landsat-5 GSTDN orbit solutions result solely from the 
small number of tracking stations. One very interesting 
statistic is that, for the Landsat-4 GSTDN tracking inter- 
vals (Figure 2-5) , an average of 88 percent of the available 
stations and only 50 percent of the available tracking 
passes were accepted and thus contributed to the solu- 
tion.'*" In contrast, for the Landsat-5 GSTDN tracking in- 
tervals (Figure 2-4) , all the available stations and all the 
available tracking passes were accepted and thus contribute 
to the solution. This may indicate that the Landsat-5 solu- 
tions, because of the lower density of data, incorporated 
observations that would otherwise be edited out. 

The present results and the previous Landsat-4 experiences 
lead to the conclusion that more stations, although each 
contains random errors and possibly even systematic errors, 
result not only in better observation geometry but ‘also in a 
larger number of observations. A large amount of data is 
needed because many observations will be edited out because 
of validity, elevation, or 3a editing. After this edit- 
ing, if a sufficient amount of data remains, the remaining 
error sources can be averaged out in the orbit determination 
process to produce small ephemeris overlap differences. In 
short, more stations provide more measurements, more observ- 
ability, and better statistics in dealing with existing ob- 
servation error sources. 


^Available refers to observat ions--associated with a par- 
ticular tracking pass or station — that are both present m 
the input data and not deleted by user observation accept/ 
reject criteria; accepted data must meet the previous avail- 
able criteria and must also not be edited out of the orbit 
solution due to either validity editing, elevation editing, 
or 3a editing. 
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2.5 ANALYSIS OF POSSIBLE MISMATCHES BETWEEN GLI AND R&D GTDS 


This section presents the results of an analysis to deter- 
mine if possible inconsistencies or inaccuracies introduced 
in the R&D GTDS- orbit determination processing cause the 
large pseudorange measurement residuals seen in both 
Landsat-4 and Landsat-5 GPSPAC data orbits (as described in 
Section 2.2). In the course of this analysis, the R&D GTDS 
code, which processes GPSPAC data, was investigated for cor- 
rectness and consistency with the mathematical specifica- 
tions of the GPSPAC/Landsat-D Interface System (GLI) 
(Reference 2-5). As noted in Section 2.1, the GPSPAC ex- 
tracts GPS data from telemetry that is in turn preprocessed 
to transform it from Earth-fixed to inertial coordinates. 

In addition to other functions, GLI performs this preproc- 
essing of the Landsat GPSPAC data for use in R&D GTDS. 

This investigation verified the correctness of the R&D GTDS 
code that processes the GPSPAC data. However, a coordinate - 
frame mismatch can occur when the coordinate reference frame 
of the R&D GTDS orbit integrator does not match the inertial 
frame of the GPSPAC data's internal state, which is included 
with the observations. For instance, in run ID G4A of 
Table 2-1, the Landsat-5 GPS orbit solution is integrated in 
a true-of-ref erence-date (TOR) coordinate frame (Refer- 
ence 2-6) referenced to April 13, 1984. However, the obser- 
vations' GPS inertial position and velocities are expressed 
in a coordinate frame referenced to April 8, 1984, the date 
EPOCH of the GPSPAC data observation file. Because the R&D 
GTDS does not check for reference frame inconsistency be- 
tween the GPSPAC observations and the integration reference 
frame, it fails to take into account the Earth's precession 
and nutation effects (Reference 2-7) occurring in the 5 days 
between April 8 and April 13. 
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To correct this reference frame inconsistency, each GPS ob- 
servation's position and velocity (the S and SDOT) vector 
(expressed in the GPSPAC EPDATE TOR frame) is rotated by a 
stand-alone routine to the reference date to be used by the 
R&D GTDS orbit integrator. For the G4A run, this rotation 
was from the April 8, 1984, to the April 13, 1984, TOR coor- 
dinate frame. This rotated GPSPAC data were then processed 
in an R&D GTDS DC run. The resultant Landsat-5 GPSPAC orbit 
solution still displays the same large pseudorange residuals 
previously observed. In fact, the converged orbit solution 
results in a state vector of the same magnitude as that gen- 
erated in the nonrotated G4A run, except that it was rotated 
through a small angle. This result was to be expected be- 
cause the reference frame rotation does not affect the 
individual GPSPAC pseudorange and delta-pseudorange measure- 
ments. Instead, it rotates the entire NAVSTAR GPS multi- 
satellite constellation through a fixed angle, which in this 
case corrects for the precession and nutation effects occur- 
ring between April 8 and April 13. 

The comparison plot between the nonrotated G4A orbit and the 
G4A orbit with no reference frame inconsistency is shown in 
Figure 2-6. The comparison displays a constant difference 
in the spacecraft along-track direction. Because Landsat 
has a 98-deg- inclination orbit, this offset is due primarily 
to nutation effects. The difference in the cross-track di- 
rection, due primarily to precession effects, displays a 
sinusoidally varying function of period roughly equal to 
Landsat's own 98-minute period. The zeros in the difference 
function occur approximately over the poles while the ex- 
tremes occur near the equatorial crossings. 

Table 2-9 summarizes the characteristics of comparisons 
between nonrotated and reference rotated GPSPAC solutions 
for Landsat-4 and -5. The solutions were performed for 
Landsat-4 arcs to determine if reference rotation would 
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Figure 2-6. Ephemens Comparison Between GPS Orbit (Run G4A) Solutions With and 
Without Reference Frame Mismatch Correction (3 of 3) 
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Table 2-9. Reference Rotated GPSPAC Data Orbit Solution Characteristics and 
Comparisons With Nonrotated Solutions 


RUN 

ID 

ARC START 
TIME 

(date hour) 

ARC 

LENGTH 

(hours) 

NO OF 
OBSERVA 
TIONS 

SOLVE FOR PARAMETERS 

RESIDUAL 

STATISTICS. 

STANDARD 

DEVIATIONS 

MAXIMUM POSITION DIFFERENCES 
OVER 4 HOUR OVERLAP OF ADJACENT 
ARCS (meters) 

MAXIMUM POSITION DIFFERENCES 
FROM GSTDN ORBITS (meters) 

REMARKS 


B 

(10“ 3 sec) 

B 

(10 -8 sec/ 
sec) 

RANGE 

(meters) 

DELTA 

RANGE 

(cm) 

RADIAL 

C 

ALONG 
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improve the GPSPAC-GSTDN ephemeris overlap compares. In 
most cases, very little difference is seen. Contrary to 
what was expected, in four instances, the reference frame 
rotation gives larger GPSPAC-GSTDN differences; although 
this may be caused by some geometric effect, these results 
are not fully understood and should be studied further. In 
any case, the rotation of the reference frame does not ap- 
pear to cause, or affect, the large pseudorange residuals 
seen in any of the GPSPAC delta-range solutions. 

A possible processing mismatch occurs if the GLI and R&D 
GTDS programs use different Greenwich hour angle (GHA) con- 
stants. To account for error on the order of 100 m, the GHA 
correction is estimated to be a rotation of roughly the same 
order of magnitude as that used to correct the coordinate 

reference frame inconsistency, that is, approximately 2 x 
-4 

10 deg. As was expected, the correction of possible GHA 
errors resulted in comparison plots analogous to the refer- 
ence frame rotation correction. 

This study also considered other conceivable errors that can 
produce range error signatures on the order of 100 m and 
appear GPSPAC range-dependent, but are delta-range insensi- 
tive--in short, the range error signatures match the trend 
of the pseudorange residual errors seen in both the 
Landsat-4 and -5 GPSPAC solutions. Two possible error 
sources, which under certain scenarios meet this error sig- 
nature and can be corrected for in the R&D GTDS processing, 
are time tag errors associated with a measurement or with 
the transformation from Earth-fixed to inertial coordi- 
nates. When these postulated errors were fixed, in isola- 
tion or in combination, the characteristic pseudorange 
residual errors in the R&D GTDS orbit solutions were in- 
creased or were only mildly affected, but they were still 
present. 
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2.6 CONCLUSIONS 


The major conclusions concerning the accuracies of GPSPAC 
data for computing Landsat-5 orbits can be summarized as 
follows : 

1. Landsat-5 orbit determination using GPSPAC data 
Landsat-5 orbit solutions computed from GPSPAC delta- 
pseudorange (Doppler) data are good. Maximum differences 
between GPSPAC and GSTDN solution are generally under 70 m. 
Maximum differences between partially overlapping GPSPAC 
solutions are even smaller. There is a good possibility 
that the GPSPAC solutions are superior to the GSTDN solu- 
tions in accuracy. 

2. For orbit determination instead of real-time navi- . 
gation, simultaneous data from four GPS satellites are not 
necessary. The results indicate that approximately 3 hours 
of data from a single GPS satellite are sufficient to re- 
solve the Landsat-5 orbit and clock. A study of a randomly 
selected sample shows that the Landsat-5 orbit solutions 
based on individual GPS satellites agree to within 80 m, and 
typically less than 40 m. 

3. As with the case of Landsat-4 GPSPAC data, there 
exist some inconsistencies between the Landsat-5 GPSPAC 
pseudorange and delta-pseudorange data. Landsat-5 solutions 
derived from pseudorange data generally differ from the 
GSTDN solutions by maximas over 100 m. Furthermore, large, 
GPS- independent, range observation residuals of over 100 m 
are seen in delta-pseudorange Landsat-5 GPSPAC data solu- 
tions. Based on these, it may be concluded that the pseudo- 
ranges have systematic errors on the order of 100 m m 
addition to the expected clock errors. The cause of these 
systematic errors has not been determined, although bad data 
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from an individual GPS satellite and several easily com- 
mitted preprocessing errors have been eliminated as possible 
reasons . 

4. In connection with the GPSPAC study, the computa- 
tion of Landsat-5 orbits from GSTDN data was undertaken. 
Unfortunately, there are not as many ground tracking sta- 
tions for Landsat-5 as for Landsat-4, and the accuracy of 
the resulting Landsat-5 orbit solutions is not as good as 
that of the corresponding Landsat-4 solutions. Investiga- 
tion shows that this degradation cannot be attributed solely 
to dynamic modeling errors accentuated by the sparsity of 
tracking coverage. 
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SECTION 3 - LOW-ALTITUDE SATELLITE NAVIGATION USING 
SIMULATED TDAS ONE-WAY DOPPLER DATA 


This study investigated the extent of agreement between 
orbit determination results using simulated TDAS one-way 
Doppler data and linear error analysis results from a pre- 
vious study. The simulation and orbit determination results 
were obtained using the R&D GTDS TDAS enhancements described 
in Reference 3-1. These simulation/orbit determination re- 
sults were compared, for an analogous satellite navigation 
scenario, with error analysis results presented in Refer- 
ence 3-2 and produced using the TDAS-enhanced SEA program. 

The scenario used to generate the R&D GTDS simulated data is 
presented in Section 3.1, and the corresponding orbit deter- 
mination scenario, in Section 3.2. From the previous SEA 
error analysis results, the dominant orbit error sources for 
the satellite considered in this study were geopotential 
uncertainty, clock acceleration, and TDAS ephemeris errors. 
The contributions from these error sources, found using the 
R&D GTDS simulation/orbit determination results, and their 
comparison with the error analysis results are discussed in 
Section 3.3. Conclusions are presented in Section 3.4. 

3.1 SIMULATION SCENARIO 

3.1.1 TDAS CONFIGURATION 

The present study considered a TDAS configuration identical 
to that used in the SEA error analysis study in Refer- 
ence 3-2. The configuration consists of three geosyn- 
chronous satellites located at longitude 73 deg, 188 deg, 
and 318 deg east. These satellites are located relative to 
the continental United States as shown in Figure 3-1. Relay 
tracking of the user satellite occurs through either of the 
two frontside TDAS relays. In addition, tracking signals 
from the ground may be relayed through the 73-deg backside 
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TDAS via one of the frontside TDAS , when the target is not 
observable to the frontside TDASs. 

3.1.2 USER SATELLITE 

One user satellite at a 28.8-deg inclination and 600-km al- 
titude was considered. The satellite is considered to have 

2 

an area-to-mass ratio of 0.0027 meters per kilogram 
(m /kg) . This spacecraft orbit scenario* is similar to 
that planned for the Space Telescope mission. A nominal 
drag coefficient of C D = 2.0 is assumed. 

3.1.3 TRACKING DATA 


The R&D GTDS one-way TDAS tracking data simulate tracking 
signals originating from the ground and relayed through the 
TDAS satellites to be received and decoded for range and 
Doppler information by the user satellite. In this simula- 
tion, range-rate Doppler data are modeled as the range dif- 
ference between two consecutive integrated ranges, spaced 
10 seconds (sec) apart. To reproduce the errors present m 
the previous error analysis study, the data simulated are 
corrupted by various systematic and random errors. The 
baseline truth model used in simulating the one-way TDAS 
data is as follows: 


Description 

Geopotential 

Solar flux/atmospheric 
density 

Random measurement error 
Range 

Range-rate equivalent of 
Doppler noise 

Measurement biases 
Range 

Range-rate equivalent of 
Doppler bias 

User clock error 

• • 

B 


Model 

21-by-21 GEM-9 

Harr is-Pr iester F150, C D = 2.0 
5 m 

5 millimeters (mm) /sec 


10 m 

1 mm/sec 


10“1° (sec/sec) /day 
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Descr ipt ion 


Model 


TDAS ephemeris error 

Height 25 m 

Crosstrack 23 m 

Alongtrack 40 m 

Two alternate tracking schedules, similar to those used in 
Reference 3-3, were considered: the broadcast or beacon 

mode, in which continuous tracking is available, and the 
scheduled mode, in which tracking is scheduled for 30 min- 
utes (min) every orbit, cycled through the three TDAS satel- 
lites with a 20- to 25-min data gap in between. For 
convenience, these will be referred to as the forward-link 
beacon tracking (FLBT) and forward-link scheduled tracking 
(FLST) . 

During the tracking periods, observations were simulated at 
a rate of one pair of range and range-rate measurements 
every 20-sec (although, in the orbit solutions, the observa- 
tion weights were adjusted to result in a range-rate solu- 
tion) . The user satellite is sometimes visible to two TDAS 
satellites simultaneously. To be consistent with the error 
analysis study, it was assumed that the user satellite is 
tracked through one TDAS satellite and then switched to the 
next TDAS satellite when the first one leaves the user's 
field of view. In contrast with previous studies (Refer- 
ences 3-1 and 3-4) , no attempt was made to minimize tracking 
through the backside TDAS. The FLBT and FLST tracking 
schedules used are shown in Table 3-1; it is identical to 
that used in the error analysis study in Reference 3-2. 

3.2 ORBIT DETERMINATION SCENARIO 

An extended Kalman filter (EKF) orbit determination solution 
was used to reduce the simulated data discussed above. In 
this solution, the variables estimated were the position and 

velocity vectors, an effective atmospheric drag coefficient 

• 

( P 1 ) , and the bias and drift of the onboard clock (B and B) . 
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Table 3-1. Tracking Schedules for 600-Kilometer-Altitude , 
28-Degree-Inclination User Satellite 


FORWARD-LINK BEACON TRACKING 

FORWARD-LINK SCHEDULED TRACKING 

TDAS1 

TDAS2 

TDAS3 

TDAS1 

TDAS2 

TDAS3 

0-3 

64-95 

3-64 

100-110 

65-75 

5-15 

95-132 

166 - 200 

132-166 

205-215 

170-180 

135-145 

200-236 

269 - 305 

236-269 

305- 315 

270 -280 

240 - 250 

305-340 

372-408 

340-372 

410-420 

375 - 385 

340 - 350 

408-444 

475-510 

444-475 

510-520 

480 -490 

445-455 

510-547 

580-612 

547- 580 

615-625 

580 - 590 

550 - 560 

612-650 

684-715 

650 - 684 

720 -730 

685-695 

655- 665 

715-753 

788-819 

753- 788 

820-830 

790 -800 

755 - 765 

819-856 

890- 924 

856-890 

925 - 935 

895-905 

860 - 870 

924 - 960 

993-1028 

960-993 

1030-1040 

995-1005 

960 - 970 

1028-1064 

1096-1131 

1064-1096 

1135-1145 

1100-1110 

1065-1075 

1131-1168 

1200-1234 

1168-1200 

1235-1245 

1200-1210 

1170-1180 

1234-1271 

1304-1336 

1271-1304 

1340-1350 

1305-1315 

1275-1285 

1336-1375 

1439-1440 

1409-1439 

1375-1409 


1410-1420 

1375-1385 















The following state process noises were introduced: between 

measurements, the variances of the velocity components in- 
crease at a rate of 10 ^ square meters per second cubed 
2 2 

(m /sec )/sec; similarly, the clock drift rate variance m- 

— 6 

creases at a rate of 10 nanoseconds (nsec) per second 

2 

cubed (nsec/sec) /sec between measurements. The complete 
baseline offset force model is summarized as follows: 


Model 

21-by-21 GEM-7 

Harris-Pr iester F200, C D = 2.2 

10“2 sec 
10"® sec/sec 

300 m in x, y, z; (positions) 

30 cm/sec in x, y, z (velocities) 


Description 

Geopotential 

Solar flux/atmospheric 
density 

User clock error 
B (clock bias) 

B (clock drift) 

A priori offset in user 
satellite state (from 
truth model state) 

The baseline a priori statistics 

Description 

A prion state uncertainty 
User H, C, L 

Orbit H, C, L 

User Bias 

Clock Drift 

Drag parameter 

Gravitational constant, GM 

Process noise 

Velocity variance growth rate 
Clock drift rate variance 


are as follows: 

Standard Deviation 


500 m 
1 m/sec 

1 msec 

200 nsec/sec 
1 

0.25 parts per million (ppm) 

10“10 (m^/sec^) /sec 
10 -6 (nsec/sec) Vsec 


3.3 SIMULATION/ORBIT DETERMINATION ERROR RESULTS 
3.3.1 GEOPOTENTIAL ERROR 

To compare the simulation results of this study with the 
error analysis results of the previous study, a common geo- 
potent lal error model was chosen. The common error model, 
which can be accommodated using R&D GTDS, is the difference 
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* -4 

- * 

between the 8-by-8 GEM-1 and the GEM-9 truncated to 15 by 
15. In R&D GTDS, the contribution from the GEM-9/GEM-1 dif- 
ference model was determined using the following method. 

First, data were simulated using the baseline version of the 
truth model except that the geopotential model employed was 
the GEM-9 truncated to 15 by 15. These data were next re- 
duced using two variations of the baseline force model. The 
first variation employed the same GEM-9 15-by-15 geopoten- 
tial model used in simulating the data. In effect, this run 
incorporated no geopotential model mismatch. The second 
orbit determination run employed the 8-by-8 GEM-1, so that 
the run incorporated the GEM- 9 /GEM- 1 model mismatch. Fi- 
nally, the reduced ephemeris files of the two orbit 
solutions were compared, with statistics generated on the 
user's state vector differences over time. To the extent 
that nonlinear effects can be ignored, the differences in 
the ephemeris comparisons should arise solely from the 
GEM- 9 /GEM- 1 model mismatch. 

Very close agreement was seen between the R&D GTDS simula- 
tion and the previous SEA error analysis results for geopo- 
tential errors arising from the GEM-9/GEM-1 difference 
model. In the error analysis study, an FLBT maximum posi- 
tion error for the GEM-9/GEM-1 difference was reported as 
38 m; in this study, 42 m was found using the simulation/ 
orbit determination method discussed previously. As shown 
in Figure 3-2, the geopotential errors display the charac- 
teristic signature that the projections along the spacecraft 
h, C, and L directions are approximately the same, with the 
H projection being the smallest. In addition, the rss total 
geopotential error contribution is observed to be fairly 
constant in time. These observations are in agreement with 
the error signatures for the GEM-l/GEM-9 errors found from 
the SEA error analysis results. 
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Figure 3-2. Geopotential (GEM- 9/GEM- 1 Difference Model) Error Projections in the 
User Orbit Plane as Determined From R&D GTDS FLBT Simulation (3 of 3) 



3.3.2 CLOCK ACCELERATION ERROR 


The clock acceleration error was determined using the fol- 
lowing method. Data were first simulated using the baseline 
truth model and reduced with the baseline force model; a 
converged user ephemeris file was then generated. Next, 
data were again simulated, using the same tracking schedule, 
but the baseline truth model was modified to exclude user 
clock acceleration offset. These simulated data were then 
reduced, using the same baseline force model as before, and 
a converged ephemeris file was generated. As the final step 
in determining the user clock acceleration error, the dif- 
ferences between these two user ephemeris files was calcu- 
lated. 

Very good agreement was found for the errors due to clock 
acceleration between R&D simulation and the SEA error analy- 
sis results of the previous study. In this study, a maximum 
clock acceleration error for FLBT tracking was 26 m; in .the 
previous study, the maximum observed was 30 m. In the error 
analysis study, it was observed that the clock acceleration 
error contribution to the range-rate measurement error in- 
creases with elapsed time and is generally linear and in the 
spacecraft alongtrack direction. This study's simulation 
runs, which are range-rate solutions, are in agreement with 
those observations. 

It should be noted that the present study's R&D solution 
used an EKF , whereas the previous study used a Kalman Filter 
(KF) Even with this difference, the simulation and 
error analysis results show agreement m the interaction 


iThe EKF, a variation to the KF, differs in that the solve- 
for state vector is corrected at each observation instead of 
waiting until the last observation. Descriptions of the EKF 
and KF in terms of their implementation in R&D GTDS are pro- 
vided in Reference 3-5. 
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between clock process noise and clock acceleration error 
observed. The R&D GTDS FLBT results are 26 in compared with 
SEA’s 30 m for clock acceleration error. Figure 3-3 dis- 
plays ephemeris comparison plots of the contribution of ac- 
celeration error resulting from the R&D GTDS FLBT mode. 

In addition, simulation and error analysis runs display con- 
sistent results if clock drift process noise is not used to 
compensate for clock errors in the filter. In SEA, when 
clock drift process noise is absent, clock acceleration 
causes an error reaching 200 m in an FLBT range-rate orbit 
solution after 24 hours. In R&D GTDS, simulations under the 
same scenario give a similarly inflated clock acceleration 
error of 132 m. 

3.3.3 TDAS EPHEMERIS ERROR 

R&D GTDS and SEA use slightly different TDAS ephemeris error 
modeling initialization. As described in Reference 3-1, the 
satellite ephemeris error model present in SEA is used in 
R&D GTDS for TDAS satellites. A description of this ephem- 
eris error model and its input is presented in Refer- 
ence 3-5. In the R&D GTDS implementation, to allow the 
greatest flexibility possible, the input phase angles 9 
and <1> are specified separately for each TDAS relay. In 
contrast, in SEA, only a single seed--input into a random 
number generator-- is used to specify the phase input ini- 
tialization for the up to 18 GPS satellites possible. The 
output of the SEA random number generator has been deter- 
mined, but because these numbers are real and the input 
phases into R&D GTDS are integers, an error of up to 0.5 deg 
can occur. This input mismatch can result in small incon- 
sistencies between R&D GTDS and SEA ephemeris error contri- 
butions. In the case of small phase angles, this mismatch 
is enlarged because of the trigonometric functions in- 
volved. Nevertheless, close agreement is seen between the 
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simulation and error analysis results for the contribution 
of TDAS ephemeris error. 

In R&D GTDS , the method of determining the contribution to 
TDAS ephemeris error is similar to that used m determining 
clock acceleration error. Specifically, instead of using 
the baseline truth model without clock acceleration, data 
are simulated without TDAS ephemeris error. The data are 
reduced using the baseline force model, and the converged 
ephemeris is then compared with an ephemeris generated using 
the baseline truth and force models. 

From SEA, the contribution from TDAS ephemeris error is 
23 m; from R&D GTDS simulation results for the FLBT mode, 
this error is 19 m. The ephemeris error, as shown in Fig- 
ure 3-4, has its largest contribution in the alongtrack di- 
rection where it increases rapidly at first but then appears 
to saturate or display very slow growth. Similar signatures 
are seen in the SEA error analysis results. 

3.3.4 ERROR RESULTS WITH LESS FREQUENTLY SAMPLED DATA 

Previous error analysis results show that systematic error 
is the dominant error source. Little performance degrada- 
tion would therefore be expected if the tracking data are 
reduced by sampling. Error analysis results also indicate 
that there is very little difference in navigation perform- 
ance if tracking data are sampled and processed at 3-min 
intervals instead of the 20-sec intervals used in the base- 
line scenario. This observation is similarly verified in 
this study's simulation results. In varying the sample 
rate, the observation weights must, however, be properly 
scaled to offset the increase in the navigation error covar- 
iance resulting from the rarefied tracking data. Alter- 
nately, as employed in the SEA error analysis study, the 
change in sample rate may be compensated for through scaling 
the process noise. 
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3.3.5 COMPARISON BETWEEN FLBT AND FLST TRACKING MODES 

In agreement with the SEA error analysis results, this 
study's FLST simulation results displayed navigation per- 
formance inferior to that in the FLBT mode. Table 3-2 con- 
trasts the error contributions in the two modes and their 
comparison with the results obtained from the SEA error 
analysis. The error analysis and simulation results display 
the same increasing trend in geopotential error contribution 
from FLBT to FLST mode. R&D and SEA results also show that 
clock acceleration error is reduced in the FLST mode. 

As discussed in the error analysis study, caution must be 
exercised in comparing FLST and FLBT results directly be- 
cause the two modes have different amounts of tracking but 
nonetheless use the same baseline process noise. In addi- 
tion, in comparing R&D GTDS simulation results with the pre- 
vious SEA error analysis results, a certain amount of caution 
must be exercised because R&D GTDS employs an EKF whereas 
SEA uses a KF. The difference in processing between the two 
types of filters may cause some of the differences between 
the error analysis and simulation results, specifically the 
errors seen for TDAS ephemeris. 

3.4 CONCLUSIONS 

Surprising agreement is seen between R&D GTDS simulation re- 
sults and SEA error analysis results for a 600-km altitude 
and 28-deg-inclination spacecraft for both the FLST and FLBT 
modes. The agreement is a confirmation of the error analy- 
sis results and provides confidence in the methods used in 
the SEA error analysis program and the meaningfulness of its 
results. Because it directly parallels the actual orbit 
determination process, there is no direct substitution for 
the method of simulation and orbit determination used m 
this study. In comparison to SEA, simulation systems such 
as R&D GTDS are much larger and more environmentally 
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Table 3-2. Navigation Error Comparison Between R&D GTDS 
Simulation and Previous SEA Error Analysis 
Results 


ERROR SOURCE 

MAXIMUM POSITION ERROR (m) 
DURING 24 HOURS OF NAVIGATION 
FOR TWO TRACKING MODES 

FLST 

FLBT 

SEA a 

R&D GTDS 

SEA a 

R&D GTDS 

20 

SEC 

180 

SEC 

20 

SEC 

20 

SEC 

180 

SEC 

20 

SEC 

GEOPOTENTIAL (GEM-9/GEM-1 







DIFFERENCE MODEL) 

63 

872 

93 5 

38 

42 2 

0 

CLOCK ACCELERATION 

17 

14 8 

14 8 

30 

24 3 

25 6 

TDAS EPHEMERIS 

23 

20 8 

21 4 

23 

15 5 

18 6 


a SEA RESULTS WERE REPORTED FOR ONLY THE 20-SEC OBSERVATION SAMPLING FREQUENCY 
b RUN WAS NOT PERFORMED 
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dependent, and so they are inherently more difficult to use 
and are computationally more expensive. The consistency of 
this study's simulation results with previous error analysis 
would therefore suggest that, after suitable calibration 
between simulation and error analysis has been ensured, ef- 
ficiency of effort may be obtained in performing the major- 
ity of the navigation analysis with an error analysis 
program such as SEA. 
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SECTION 4 - LUMPED GEOPOTENTIAL MODEL FOR GEM- 9 


The following spherical harmonic representation of the geo- 
potential is generally recognized as the most convenient for 
orbital analysis: 


V = 


GM 


°° Z l 

-EE t* 

Z— 2 m— 0 


JLm (sin ♦> (C ta cos mX 


+ sin mX) 


(4-1) 


where GM = geocentric gravitational constant 

a g = mean equatorial radius of the Earth 

P$, m = normalized associated Legendre function of 
degree & and order m 

r, <t>, X = distance to the center of the Earth, lati- 
tude, and longitude 

C & , = normalized spherical harmonic coefficients 

The quantity GM provides a scale for the whole field; a 


provides a scale for the altitude dependence, r; and C 


&m 


and 


characterize the geographic, or nonspher ical , varia- 
tions. A geopotential model is defined by these param- 
eters. Of necessity, the field is truncated at some finite 
degree and order. Over the years, NASA/GSFC has developed a 
series of progressively more accurate GEMs based on satel- 
lite tracking data, satellite-borne altimetry, and surface 
gravimetry. GEM-10B (Reference 4-1) , complete to degree and 
order 36 by 36, is currently the most accurate model. The 
Goddard Trajectory Determination System (Reference 4-2) , 
which is used for operational orbit determination at GSFC, 
has provisions for GEM-9 (Reference 4-3) further truncated 
to degree and order 21 by 21. 
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Even with accurate geopotential models such as GEM-9 and 
GEM-10B, the uncertainty in these models remains a major 
error source for low-altitude Earth satellite orbit determi- 
nation. Geopotential error arises because of truncation 
(omission error) and because of errors in the spherical har- 
monic coefficients (commission error) . Generally, with the 
derivation of a GEM, an error covariance matrix describing 
the uncertainties and correlations of the spherical harmonic 
coefficients is a standard byproduct. In principle, it 
would be straightforward to use this covariance matrix in 
linear error analysis to compute the effect of commission 
error on satellite orbit determination accuracy. However, 
the size 'of the covariance matrix makes this approach pro- 
hibitively expensive in computations. 

As a simpler alternative, a so-called lumped geopotential 
error model was suggested (Reference 4-4) and implemented in 
orbit determination error analysis programs. Briefly, this 
approach takes the weighted differences of the geopotential 
coefficients of two independent geopotential models and com- 
putes the orbit determination error resulting from the 
lumped effect of these differences. Subsequently, standard 
deviations, which are the scaled-up formal uncertainties of 
GEM spherical harmonic coefficients, were used in place of 
the coefficient differences. The theoretical objection to 
the use of standard deviations was first raised in Refer- 
ence 4-5. Its practical inadequacy was brought out in vivid 
graphical displays in Reference 4-6. 

Section 4.1 describes the rationale of the lumped geopoten- 
tial error model. Section 4.2 discusses the global distri- 
bution of gravity error for different geopotential error 
models. Section 4.3 discusses navigation and orbit predic- 
tion errors introduced by geopotential uncertainty. Sec- 
tion 4.4 compares results obtained using TDRS data with 
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those from the lumped geopotential error model for GEM-9. 
Section 4.5 presents conclusions. 

4.1 THE LUMPED ERROR MODEL 

The derivation of each GEM is always accompanied by a care- 
ful error analysis to assess its accuracy. Reference 4-7 on 
GEM- 7 and GEM- 8 is an outstanding example of the conscien- 
tiousness and thoroughness associated with such assess- 
ments. Analyses of GEM-9 and GEM-10 accuracies are available 
in References 4-3 and 4-1, respectively. Reference 4-8 con- 
tains a discussion of gravity model improvement and its im- 
plications for operational orbit determination. 

As discussed before, each GEM is accompanied by an error co- 
variance matrix associated with the geopotential coeffi- 
cients. The covariance matrix is computed based on assumed 
precisions of tracking and other data from which the geo- 
potential is derived. Because not all error sources in 
tracking data and in spherical harmonic coefficient estima- 
tion methods can be accounted for, the computed error covar- 
iance is not a true indication of the accuracy of the 
spherical harmonic coefficients. Part of the objective of 
geopotential model accuracy analysis is to derive a calibra- 
tion factor, which is used to scale up the standard devia- 
tions from the error covariances to more realistic levels. 

The calibration factor is typically around 3.3. 

One of the obvious methods of assessing a geopotential model 
is to compare it with other independently derived geopoten- 
tial models. The lumped geopotential error model, first 
proposed by Martin and Roy (Reference 4-4) , takes one-half 
the differences of two uncorrelated models of comparable 
accuracy as a measure of the accuracies of the individual 
models. The rationale can be explained as follows. Let 
(C a ) and (C 0 ) be the geopotential coeffi- 

cients of the two models. If the two models are uncorrelated 
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but of comparable accuracy, their average is closer to the 
truth, i . e . , 


( C !m) 


truth 


1 

~ 2 


( C|lm )model A + ( Cim ) 


model B 


(4-2) 


and an estimate of the individual errors can be obtained as 


( £m )model A ( 5,111 Imodel A ( lm )trut 

.Nlmodel A ' (M 


1 

2 


model BJ 


(4-3) 


Of course, an estimate using the sample mean based on onl^ 
two samples, as in Equation (4-3), is not very reliable. 

The lumped model introduces additional samples or statistics 
by considering not a single but the whole set of C^'s 
and S £ m 's in the geopotential model. In other words, it is 
not saying that Ac^ m as given in Equation (4-3) is a good 
indicator of the accuracy of individual coefficients, but 
that the aggregate effect of all AC^'s on the computed 
orbit is a reasonable representation of the effect of the 
geopotential uncertainty on orbit determination. 

The key assumption of the lumped model is the requirement 
that the two models be uncorrelated. This generally means 
that the models be derived by two different organizations 
based on different data sets. It also implies that this 
lumped model does not adequately represent omission errors. 
The assumption that the two models are of equal accuracy is 
not essential, as the simple average in Equation (4-3) can 
be replaced by a weighted average. In the spirit of minimum 
variance estimators, the weights can be chosen to be in- 
versely proportional to the variances of the individual 


\ 
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models. In practice, an educated guess or external calibra- 
tion can be used to assign the weights. 

The lumped model has been used in many orbit determination 
error analysis studies, perhaps not so much because it is a 
good model but because nothing better is available. Part of 
the difficulty of this approach is the availability of inde- 
pendent models. In a study (Reference 4-9) of onboard orbit 
determination using GEM-9 truncated to various degrees, it 
was suggested that the lumped model be used to represent the 
uncertainty in the truncated GEM-9 as follows: GEM-9 stand- 

ard deviations would be used to represent the uncertainties 
of those spherical harmonic coefficients retained in the 
truncated model, and truncation errors would be represented 
by the truncated high degree and order GEM-9 spherical har- 
monic coefficients themselves. There are two theoretical 
objections (Reference 4-5) to the use of standard deviations 
instead of coefficient differences in the lumped model: 

• The errors in the geopotential coefficients are 
likely to have both positive and negative signs, 
whereas the standard deviations are all positive. 

• The standard deviations do not contain information 
about correlations among the coefficients. 

It may be argued that, because the primary interest is in 
orbit determination errors, rather than in the errors in the 
geopotential coefficients, the use of the positive standard 
deviations may not be objectionable if the orbit error sen- 
sitivities to these coefficients are randomly distributed to 
serve to randomize their aggregate effect. Elrod, however, 
showed that this is not the case (Reference 4-6) . Fig- 
ure 4-1, reproduced from Reference 4-6, shows that the use 
of GEM-9 standard deviations results in a nonuniform global 
distribution of gravitational acceleration errors that is 
greater in the Northern Hemisphere with a singularity near 
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0 deg longitude and 60 deg latitude. Such anomalous distri- 
bution is not supported by other evidence. Furthermore, 
Elrod showed that the anomaly does not occur if the rss re- 
sult is obtained rather than the algebraic sum of the lumped 
error model. There is no question that the use of GEM-9 
standard deviations in the lumped error model is faulty. 

The effect of false singularity is perhaps not serious for 
error analysis of the batch orbit determination method, but 
it can exaggerate local navigation errors in sequential 
navigation methods. 

Because the standard deviations do not convey information 
about correlations, an obvious alternative is to assume that 
all errors in the spherical harmonic coefficients of a geo- 
potential are uncorrelated. The total error of the geopo- 
tential model can thus be obtained as the rss of the error 
contribution of individual coefficients. Although it is 
known that certain geopotential coefficients are highly cor- 
related, this model is perhaps as reasonable as any, short 
of having to consider the actual correlations. The major 
drawback of this model is the excessive amount of computa- 
tion required to compute separately the error contributions 
of individual error coefficients, versus the spirit of the 
lumped error model. At best, only the level of orbit deter- 
mination error resulting from geopotential uncertainty can 
be expected; therefore, the computationally expensive rss 
approach may not be justifiable, and simpler alternatives 
would be desirable. The following alternatives are proposed 
as candidates for study; 

• Random-Sign Method--The lumped model will use the 
GEM-9 standard deviations, with positive and negative signs 
randomly assigned to them. 
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• Random-Phase Method--The terms involving the spher- 
ical harmonic coefficients in the geopotential representa- 
tion, Equation (4-1) , may be written as 


c £m cos + s fi.m sin mX = Y C L + S lm cos (mX ' V (4_4) 
In this method, the amplitudes 



are computed from GEM-9 standard deviations, but the phase 
angles, <j> , are selected from a uniform random distribution 
between 0 and 2n (for zonals, the phase is randomly chosen as 
0 or tt, i.e., random sign). 

In addition to these models dependent on internal accuracy 
estimates, two geopotential models uncorrelated with GEM-9 
are used in conjunction with GEM-9 to form geopotential dif- 
ference error models. One of these is the 1969 Smithsonian 
Astrophysical Observatory (SAO) Standard Earth (Refer- 
ence 4-10) , referred to as the SAO model below. The other 
will be referred to as the MD1 model. One-half of the dif- 
ference of this model and GEM-5 has been considered repre- 
sentative of the accuracy of GEM-5 (Reference 4-11) . Both 
the SAO and MDl models are inferior to GEM-9 in accuracy. 

4.2 GLOBAL DISTRIBUTION OF GRAVITY ERROR FOR DIFFERENT GEO- 
POTENTIAL ERROR MODELS 

The gravitational acceleration as specified by a geopoten- 
tial model will have errors that vary with the geographical 
location (longitude and latitude) and decrease with the al- 
titude. A gravity error map such as that illustrated in 
Figure 4-1 shows the magnitude (but not the direction) of 
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the gravitational acceleration error at a given altitude as 
a function of longitude and latitude, as predicted by a geo- 
potential error model. Different error models give rise to 
different error characteristics. Of interest are orders of' 
magnitude and geographical fluctuations of the errors. Fig- 
ures 4-2 through 4-22 are gravity error maps at 200-, 400-, 
and 600-km altitudes for the following geopotential error 
models: 


• GEM-9/MD1 one-half difference model 

• GEM- 9 /S AO one-half difference model 

• GEM- 9 /GEM- 5 difference model 

• GEM-5/MD-1 one-half difference model 

• GEM-9 standard deviation model 

• GEM-9 random sign model 

• GEM-9 random-phase model 

• GEM-9 uncorrelated model 

The difference models, discussed in Section 4.1, are self- 
explanatory. The GEM-9 standard deviation model uses posi- 
tive GEM-9 standard deviations as spherical harmonic error 
coefficients. The gravity error map for this model at 
200-km altitude is shown in Figure 4-1, so only the maps at 
400- and 600-km altitudes are shown here. The random-sign 
and random-phase models are also explained in Section 4.1. 

A number of variants of these models were also considered, 
but they did not show significant difference from those de- 
scribed here. The uncorrelated model is based on the rss of 
the contributions of individual coefficients; computing these 
errors is expensive, and the gravity error maps were taken 
from Reference 4-6. To provide a reference for gauging the 
magnitudes of gravity errors, Figure 4-23 shows a map repre- 
senting the gravitational acceleration (not the error) re- 
sulting from the nonaxisymmetr ic (i.e., the nonzonal) portion 
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Figure 4-2. GEM-9/MD1 One-Half Difference Model, 


200-Kilometer Altitude 
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Figure 4-3. GEM-9/MD1 One-Half Difference Model, 400-Kilometer Altitude 
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Figure 4-4. GEM-9/MD1 One-Half Difference Model, 600-Kilometer Altitude 
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Figure 4-5. GEM-9/SA0 One-Half Difference Model, 200-Kilometer Altitude 
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Figure 4-6. GEM-9/SAO One-Half Difference Model, 400-Kilometer Altitude 
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Figure 4-7. GEM-9/SA0 One-Half Difference Model, 600-Kilometer Altitude 



<L A T> 


UNCERTAINTY IN NONSPHERICAL GRAVITATIONAL ACCFIERAtlON IMGAL) 


90 

♦ 

1 

1 

i 

i 

1 

i 

• 

i 

i 

i 

i 

i 

i 

i 

i 

2 

2 

2 

2 

2 

2 

2 

i 

1 

i 

t 

i 

1 

i 

i 

i 

i 

i 

1 

1 

1 

i 

85 

+ 

1 

1 

i 

1 

i 

i 

1 

1 

0 

i 

1 

i 

2 

2 

3 

3 

3 

3 

3 

3 

3 

3 

3 

2 

2 

2 

1 

1 

i 

t 

1 

1 

1 

1 

1 

1 

1 

80 

♦ 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

3 

4 

4 

4 

4 

4 

4 

4 

4 

4 

3 

3 

2 

2 

2 

2 

2 

2 

2 

2 

2 

75 

♦ 

2 

2 

3 

3 

3 

3 

3 

3 

2 

2 

2 

2 

2 

2 

2 

1 

2 

4 

5 

5 

5 

4 

4 

4 

5 

6 

6 

6 

6 

6 

7 

7 

6 

5 

3 

2 

2 

70 

+ 

3 

4 

4 

3 

4 

3 

3 

3 

3 

2 

2 

2 

2 

3 

4 

4 

5 

6 

6 

7 

7 

8 

8 

7 

6 

7 

7 

7 

8 

9 

10 

1 0 

8 

6 

4 

3 

3 

65 

4 

5 

6 

5 

3 

3 

4 

4 

5 

6 

5 

5 

4 

4 

4 

5 

6 

8 

9 

1 0 

10 

1 0 

13 

13 

1 0 

7 

8 

9 

8 

7 

9 

1 0 

9 

7 

6 

4 

4 

5 

60 

♦ 

4 

4 

3 

2 

4 

7 

8 

6 

5 

6 

6 

5 

5 

6 

7 

9 

(0 

t t 

1 2 

1 1 

1 1 

15 

16 

1 0 

5 

6 

8 

6 

3 

6 

8 

8 

7 

4 

3 

3 

4 

55 

♦ 

5 

4 

2 

4 

1 0 

1 2 

1 1 

6 

3 

8 

8 

6 

6 

8 

1 0 

1 3 

1 2 

1 1 

1 1 

8 

8 

1 3 

1 4 

9 

5 

4 

7 

6 

2 

4 

5 

7 

1 0 

6 

3 

3 

5 

50 

1- 

6 

5 

2 

7 

16 

19 

1 1 

7 

4 

8 

8 

6 

8 

1 2 

1 4 

16 

1 1 

8 

9 

6 

4 

9 

8 

8 

6 

4 

6 

8 

8 

7 

4 

6 

1 2 

9 

8 

4 

6 

45 

♦ 

5 

6 

3 

1 0 

1 4 

19 

1 0 

7 

8 

1 t 

7 

4 

8 

1 3 

1 2 

1 4 

1 0 

6 

8 

7 

8 

1 0 

4 

6 

7 

5 

5 

9 

1 2 

8 

4 

5 

9 

9 

9 

5 

5 

40 

♦ 

1 

5 

5 

9 

9 

t 1 

7 

6 

7 

1 3 

6 

3 

5 

10 

6 

7 

8 

6 

8 

6 

9 

1 1 

6 

7 

9 

6 

4 

7 

1 1 

7 

4 

3 

5 

4 

6 

3 

1 

35 

» 

6 

5 

5 

7 

5 

7 

1 0 

6 

5 

9 

4 

3 

2 

6 

3 

3 

4 

5 

8 

6 

6 

8 

8 

7 

8 

6 

4 

6 

8 

5 

3 

7 

7 

6 

9 

7 

6 

30 

4 

7 

4 

7 

8 

4 

7 

t 1 

1 0 

7 

7 

2 

3 

4 

5 

5 

3 

3 

9 

7 

4 

4 

4 

7 

5 

6 

5 

4 

8 

5 

3 

3 

9 

7 

9 

1 1 

7 

7 

25 

4 

5 

4 

7 

7 

4 

6 

7 

8 

6 

6 

3 

3 

5 

4 

5 

5 

2 

10 

6 

6 

2 

2 

4 

4 

3 

4 

5 

1 0 

8 

5 

2 

6 

4 

6 

9 

3 

5 

20 

4 

2 

3 

9 

4 

2 

4 

5 

6 

2 

4 

4 

1 

3 

5 

5 

7 

2 

8 

5 

7 

2 

2 

5 

5 

2 

4 

6 

8 

1 1 

6 

2 

2 

4 

5 

1 0 

4 

2 

1 5 

4 

3 

5 

9 

4 

1 

3 

6 

1 0 

3 

5 

7 

5 

4 

a 

3 

5 

2 

6 

2 

7 

5 

2 

6 

4 

2 

2 

4 

4 

7 

5 

3 

9 

6 

5 

9 

5 

3 

1 0 

4 

3 

4 

5 

4 

0 

1 

8 

1 0 

4 

6 

9 

6 

7 

6 

3 

4 

2 

5 

3 

5 

4 

3 

4 

5 

5 

| 

4 

4 

4 

5 

5 

1 1 

6 

4 

6 

4 

3 

5 

4 

2 

2 

5 

4 

1 

3 

6 

8 

6 

6 

1 0 

2 

7 

2 

2 

5 

5 

6 

5 

4 

2 

1 

4 

6 

7 

2 

5 

8 

4 

4 

3 

8 

6 

4 

6 

2 

2 

0 

* 

1 

3 

7 

5 

3 

5 

4 

7 

7 

4 

8 

3 

5 

4 

2 

6 

5 

6 

5 

3 

3 

3 

7 

G 

6 

4 

7 

8 

4 

4 

2 

4 

6 

4 

7 

4 

1 

-5 

4 

5 

5 

6 

6 

5 

4 

5 

7 

4 

2 

6 

4 

3 

4 

3 

5 

3 

6 

6 

2 

1 

4 

10 

4 

5 

4 

6 

4 

3 

3 

3 

3 

5 

6 

6 

6 

5 

- 1 0 

4 

7 

4 

3 

5 

6 

2 

6 

6 

4 

3 

5 

4 

3 

2 

3 

7 

3 

6 

6 

3 

3 

4 

7 

6 

5 

3 

4 

1 

2 

6 

2 

1 

5 

7 

5 

6 

7 

- 1 5 

4 

4 

2 

3 

5 

7 

1 

5 

5 

4 

4 

4 

3 

5 

3 

2 

6 

1 

5 

5 

5 

2 

4 

6 

6 

4 

2 

8 

1 

3 

7 

1 

3 

4 

7 

6 

5 

4 

-20 

4 

2 

2 

2 

5 

5 

3 

3 

6 

3 

4 

4 

2 

4 

3 

3 

3 

3 

5 

3 

5 

2 

2 

5 

5 

5 

6 

9 

2 

4 

6 

3 

4 

4 

6 

6 

6 

2 

-25 

• 

3 

2 

2 

5 

6 

6 

4 

9 

3 

2 

3 

2 

2 

3 

5 

4 

4 

5 

3 

4 

3 

2 

5 

4 

5 

5 

6 

3 

4 

3 

3 

3 

4 

7 

3 

7 

3 

-30 

4 

4 

2 

2 

4 

5 

5 

3 

8 

4 

0 

2 

2 

1 

3 

9 

9 

5 

4 

4 

5 

2 

3 

4 

5 

4 

4 

2 

3 

2 

2 

2 

3 

3 

7 

5 

8 

4 

- 35 

4 

5 

3 

4 

5 

6 

4 

2 

5 

5 

3 

2 

3 

3 

3 

1 0 

1 0 

7 

5 

6 

7 

7 

5 

2 

3 

2 

3 

4 

3 

3 

2 

2 

3 

4 

7 

8 

8 

5 

- 40 

4 

7 

6 

8 

6 

6 

6 

4 

4 

5 

5 

4 

5 

4 

4 

5 

7 

8 

9 

10 

9 

9 

4 

1 

4 

4 

3 

3 

3 

4 

5 

3 

4 

3 

7 

1 0 

8 

7 

-45 

4 

1 0 

8 

7 

5 

5 

7 

6 

5 

4 

3 

6 

7 

3 

4 

3 

5 

6 

8 

10 

9 

7 

2 

3 , 

5 

4 

3 

3 

1 

3 

5 

4 

5 

3 

5 

9 

8 

1 0 

-50 

4 

9 

9 

8 

5 

4 

5 

6 

6 

5 

1 

5 

8 

6 

6 

3 

5 

6 

4 

5 

4 

5 

4 

4 

5 

3 

3 

2 

2 

3 

2 

3 

5 

4 

2 

4 

5 

9 

-55 

4 

5 

9 

1 0 

5 

3 

5 

3 

5 

5 

2 

3 

5 

6 

7 

4 

4 

5 

2 

3 

4 

7 

7 

4 

3 

2 

4 

3 

2 

5 

4 

5 

5 

8 

5 

4 

3 

5 

-60 

4 

6 

8 

6 

4 

2 

3 

4 

4 

4 

3 

2 

4 

6 

6 

6 

6 

6 

5 

4 

4 

6 

8 

5 

4 

4 

6 

4 

3 

6 

6 

6 

6 

9 

8 

7 

4 

6 

-65 

4 

7 

7 

5 

5 

5 

4 

4 

1 

2 

2 

t 

3 

4 

5 

.6 

9 

9 

8 

6 

3 

3 

7 

9 

9 

7 

7 

6 

6 

7 

6 

5 

5 

7 

8 

8 

7 

7 

-70 

4 

7 

1 0 

t 1 

1 1 

1 0 

7 

5 

3 

2 

2 

2 

4 

5 

6 

7 

8 

9 

8 

7 

5 

6 

8 

10 

1 0 

9 

9 

9 

1 0 

9 

8 

7 

6 

7 

7 

6 

6 

7 

-75 

4 

8 

1 0 

t t 

1 1 

9 

8 

6 

4 

4 

4 

4 

6 

7 

8 

8 

7 

6 

4 

3 

2 

4 

6 

8 

9 

9 

9 

9 

1 0 

9 

8 

8 

8 

9 

9 

8 

8 

8 

-80 

4 

a 

8 

7 

6 

5 

4 

3 

3 

3 

3 

4 

5 

5 

6 

6 

6 

5 

4 

3 

3 

4 

6 

7 

7 

8 

8 

8 

8 

8 

8 

9 

9 

1 0 

1 0 

1 0 

9 

8 

-85 

4 

6 

5 

4 

4 

3 

2 

1 

1 

1 

1 

2 

2 

3 

3 

3 

4 

4 

4 

5 

5 

5 

6 

6 

6 

6 

; 

-» 

7 

7 

7 

7 

7 

7 

7 

7 

6 

6 

90 

4 

4 

4 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 


I - 

- I - 

- I - 

• - I - 

- I - 

• - 1 - 

- I - 

- I - 

- I - 

■ - 1 - 

- I 

- - r - 

- 1 

- - 1 - 

- 1 

- * I ' 

• I 

- - I - 

■ - I 

- - 1 - 

- I 

- - 1 

I 

1 - 

I 

i 

I 

- - I - 

- 1 

- - I - 

- 1 

- I - 

- 1 

- 1 

- I 

- I - 

- I 

- - I 

LONG> 


0 


20 


40 


60 


80 


1 00 


120 


140 


160 


180 


200 


??n 

7 

40 


260 


280 


300 


320 


340 


360 


Figure 4-8. GEM-9/GEM-5 Difference Model, 200-Kilometer Altitude 



<LAT> 


UNCERTAINTY IN NONSPHE R I C A t. GRAV I TA I I ONAI ACCFLFRA I l ON IMGALI 



Figure 4-9. GEM- 9 /GEM- 5 Different Model, 400-Kilometer Altitude 
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Figure 4-10. GEM-9/GEM-5 Difference Model, 600-Kilometer Altitude 
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Figure 4-11. GEM-5/MD1 One-Half Difference Model, 200-Kilometer Altitude 
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Figure 4-12. GEM-5/MD1 One-Half Difference Model, 400-Kilometer Altitude 
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Figure 4-13. GEM-5/MD1 One-Half Difference Model, 600-Kilometer Altitude 
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Figure 4-14. GEM-9 Standard Deviation Model, 400-Kilometer Altitude 
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Figure 4-16. GEM-9 Formal Uncertainties With Random Signs, 200-Kilometer Altitude 
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Figure 4-17. GEM-9 Formal Uncertainties With Random Signs, 400-Kilometer Altitude 
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Figure 4-18. GEM-9 Formal Uncertainties With Random Signs, 600-Kilometer Altitude 
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Figure 4-19. 


GEM-9 Formal Uncertainties With Random Phase, 200-Kilometer 
Altitude 
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UNCERTAINTY IN NONSPHERICAL GRAVITATIONAL ACCELERATION (MGALt 



Figure 4-20. GEM-9 Formal Uncertainties With Random Phase, 400-Kilometer 
Altitude 
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Figure 4-21. GEM-9 Formal Uncertainties With Random Phase, 600-Kilometer 
Altitude 
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Figure 4-23. 


Gravity Map at a 200-Kilometer Altitude for the Nonaxisymmetric 
Portion of GEM-9 



of the GEM-9 spherical harmonic coefficients. The maps 
depict the following: 

• As expected, gravity errors for every model de- 
crease with altitude. 

• In addition to the anomalous distribution of grav- 
ity errors discussed for the GEM-9 standard deviation model, 
difference models involving the MD1 and a GEM model show a 
band of large differences at latitudes near the Equator. 

• The GEM-9 uncorrelated model predicts a very uni- 
form distribution of errors with a magnitude around 

-5 2 

4 milligals (mgals) (10 m/sec ) at 200 km. The GEM-9 
standard deviation model predicts errors, apart from the 
anomalous region, at -a level considerably smaller than those 
given by the other models. The other models generally give 
rise to errors from 1 to 10 mgals and show much greater 
fluctuations from one location to another. 

4.3 NAVIGATION AND ORBIT PREDICTION ERRORS INTRODUCED BY 
GEOPOTENTIAL UNCERTAINTY 

The interest in geopotential error models is not so much in 
the gravity errors themselves but in their effect on satel- 
lite navigation accuracy. Satellite navigation error 
analyses were thus conducted using the SEA program to eval- 
uate, under postulated navigation scenarios, the expected 
orbit determination and prediction errors resulting from the 
geopotential errors given by different error models. The 
baseline scenario considered was for low-altitude Earth 
satellite navigation based on a sequential filter using con- 
tinuous TDAS range and Doppler tracking data sampled at 
3-min intervals for 1 day, followed by 5 days of propagation 
without tracking. Navigation at 200 and 600 km and at 28- 
ana 57-deg inclinations was considered. No process noise 
was introduced in the filter; thus, the 5-day propagation 
results should agree with those based on a batch orbit 
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determination process. The expected navigation errors re- 
sulting from geopotential uncertainties are summarized in 
Table 4-1 for the different geopotential error models. 

4.4 COMPARISON WITH LANDSAT-5 ORBIT DETERMINATION RESULTS 

Some results on Landsat-5 orbit determination using GSTDN 
and single-TDRS tracking data are presented in Refer- 
ence 4-12. Measures of orbit accuracies are provided by 

• Overlap ephemeris comparisons of definitive orbit 
solutions (34-hour data arcs with 10 hours of over- 
lap) 

• Ephemeris comparisons of definitive and predictive 
orbit solutions (1- or 2-day predictions based on 
34-hour data arcs) 

The overlap comparisons of definitive solutions derived from 
TDRS data have maximum differences of 82, 71, 35, 74, and 
35 m. One-day predictions have maximum differences of 186, 
131, 40, 172, and 79 m from definitive solutions. Two-day 
predictions have maximum differences of 412 and 172 m. Al- 
though these differences, or errors, result from a combina- 
tion of geopotential and atmospheric uncertainties and 
tracking errors, it is expected that geopotential uncertain- 
ties play a major role. These numbers can thus serve as 
bases for the calibration of geopotential error models. For 
this purpose, error analyses were set up to evaluate, under 
the same tracking and batch orbit determination scenario, 
orbit errors resulting from the GEM-9/SAO one-half differ- 
ence error model and the GEM-9/MD1 one-half difference 
model. The GEM-9/SAO error model predicts a maximum posi- 
tion error of 57 m during the tracking period and maximum 
errors of 184 and 236 m during 1-day and 2-day propaga- 
tions. These numbers appear quite reasonable compared with 
actual orbit determination results quoted above. The GEM-9/ 
MD1 error model predicts a maximum position error of 77 m 
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Table 4-1. Navigation Errors Resulting From Geopotentxal 
Uncertainties 


ORBIT 

CHARACTERISTICS 

GEOPOTENTIAL ERROR 
MODEL 

ORBIT POSITION ERROR 
(meters) 

MAXIMUM DURING 
24 HOURS OF 
TRACKING 

MAXIMUM DURING 
5-DAY 
PREDICTION 

200-km ALTITUDE, 

GEM-9 STANDARD DEVIATION 

21 

205 

28 -deg INCLINATION 




'A GEM-9/SAO DIFFERENCE 

149 

481 


A GEM-9/MD-1 DIFFERENCE 

168 

3209 


'A GEM-5/MD-1 DIFFERENCE 

186 

3495 


RANDOM SIGN 

148 

496 


RANDOM PHASE 

108 

575 

200-km ALTITUDE, 

GEM-9 STANDARD DEVIATION 

56 

1171 

57-deg INCLINATION 




A GEM-9/SAO DIFFERENCE 

126 

1271 


'A GEM-9/MD-1 DIFFERENCE 

253 

44886 


'A GEM-5/MD-1 DIFFERENCE 

227 

42745 


RANDOM SIGN 

105 

4562 


RANDOM PHASE 

238 

2511 

600-km ALTITUDE, 

GEM-9 STANDARD DEVIATION 

9 

71 

28-deg INCLINATION 




'A GEM-9/SAO DIFFERENCE 

60 

246 


'A GEM-9/MD-1 DIFFERENCE 

73 

597 


'A GEM-5/MD-1 DIFFERENCE 

71 

576 


RANDOM SIGN 

68 

154 


RANDOM PHASE 

44 

118 

600-km ALTITUDE, 

GEM-9 STANDARD DEVIATIONS 

14 

235 

57-deg INCLINATION 




'A GEM-9/SAO DIFFERENCE 

58 

686 


'A GEM-9/MD-1 DIFFERENCE 

81 

1537 


'A GEM-5/MD-1 DIFFERENCE 

77 

2330 


RANDOM SIGN 

35 

395 


RANDOM PHASE 

44 

118 


4-34 


S8/UU) ZllO 


























during the tracking period but unrealistic, smaller maximum 
errors of 58 and 64 m during propagations. 

4.5 CONCLUSIONS 

The results described in the preceding sections show that 
the choice of an appropriate model is very difficult. The 
models studied may be classified into two categories: those 

dependent on internal accuracy estimates and those dependent 
on comparisons with other models. 

Error models based on internal accuracy estimates run the 
risk of being unduly optimistic. This weakness can, how- 
ever, be remedied by a calibration scale factor. The real 
difficulty is the need to consider correlations between 
errors in the geopotential coefficients. This is imprac- 
tical because of the large computational requirements. 

Short of actually considering the correlations, the GEM-9 
uncorrelated standard deviation model seems to be the most 
appropriate, although correlations do exist. However, even 
this model is computationally prohibitive. Comparison of 
the error map of this model with those of the geopotential 
difference models shows that the geographical fluctuation of 
the error as predicted by this model may be too mild. On 
the other hand, based on the same comparisons, it seems that 
the predicted gravity error magnitude for GEM-9, approxi- 
mately 5 mgals at a 200-km altitude, is reasonable. 

The random-sign and random-phase models, which are simple to 
implement, give rise to gravity error maps that seem reason- 
able, with local error fluctuations possibly somewhat exces- 
sive. The orbit determination errors based on these models 
also appear reasonable, although the error growths during 
prediction do not clearly exhibit the distinct characteris- 
tic of a sinusoidal variation superimposed on a linear 
growth. The main objection to these models is that they are 
empirical models without sound theoretical justifications. 
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The GEM-9 standard deviation model, when compared with other 
models, shows gravitational errors that are too low in the 
Southern Hemisphere and near the Equator. Therefore, its 
prediction of orbit errors for low-inclination satellites is 
too optimistic. On the other hand, this model has an un- 
realistic concentration of large errors near 0-deg longitude 
and 65-deg latitude that may give rise to anomalous orbit 
error spikes near these regions. This model is definitely 
faulty and should not be used. 

A comparison between different geopotential models provides 
information about the accuracies of the individual models. 
The geopotential difference error models use the weighted 
differences of the spherical harmonic coefficients of two 
geopotential models as error coefficients. Tf one of the 
geopotential models is considerably more accurate than the 
other, the straight difference of the two models can be con- 
sidered, with good confidence, as representative of the 
errors of the less-accurate model. If the two models are 
derived from independent sources, there is some justifica- 
tion in using their weighted differences as characterizing 
the accuracies of individual models, although the justifica- 
tion is weak and cannot be rigorously proven. 

SAO and MD1 are models independent of GEM-9. GEM-9 was de- 
rived several years after SAO and MD1 , and there should be 
no disagreement that GEM-9 is a considerably improved 
model. Thus, the GEM- 9 /SAO and GEM-9/MD1 difference models 
serve as good error models for the SAO and MDl geopotential 
models. Doubling the numbers shown in Figure 4-5, which are 
for one-half the model differences, implies that the gravi- 
tational errors for the SAO model at a 200-km altitude vary 
from 2 to 26 mgals. The error map shows local fluctuations 
but not any particularly recognizable features. Similarly, 
it might be said that Figure 4-2, with the errors doubled, 
describes the accuracy of the MDl model. It should be noted 

4-36 


0112 



that, in Figure 4-2, there is a band of large errors at 
equatorial latitudes reaching, at a 200-km altitude, as 
large as 60 mgals (in accuracy), which is comparable in mag- 
nitude to the total contribution of the nonaxisymmetr ic 
portion of GEM-9. This band and its large magnitude is 
bothersome, because no ready explanations exist. Away from 
this band, the error distribution appears reasonable. 

Of course, it is not the errors for these older models but 
the accuracy of GEM-9 that is of interest here. Since, at 
best, only the level of orbit determination error resulting 
from geopotential uncertainty can be expected, it may not be 
unreasonable to take the scaled difference of GEM-9 and a 
less accurate model as an error model for GEM-9, with the 
scale factor determined from some means of calibration. For 
instance, it might be speculated that GEM-9 is twice as ac- 
curate as MD1 and that therefore one-half of the GEM-9/MD1 
difference would serve as an error model for GEM-9. This 
is, of course, not completely satisfactory because gravity 
errors are expected to become more uniformly distributed for 
higher order and more accurate models. The band of large 
errors, shown in Figure 4-2, is probably reflective of the 
nonuniform inaccuracies of MD1 , and a simple scaling of the 
error model that preserves this feature of concentrated 
error would not be a good error model for GEM-9. The same 
consideration also applies to the GEM-9/SAO difference 
model. However, as shown in Figure 4-5, this difference 
model has a more random error distribution devoid of partic- 
ular unexplainable features. Thus the GEM- 9/S AO difference 
model, with a proper scaling factor, appears to be a reason- 
able error model for GEM-9. 

The primary interest is not in the global gravity error dis- 
tribution itself but in the effect of geopotential errors on 
orbit determination accuracies. Table 4-1 shows the orbit 
determination and prediction errors according to the 
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different geopotential error models. The following may be 
observed from the results presented in the table: 

• With the exception of the GEM-9 standard deviation 
model, the orbit determination errors predicted by the dif- 
ferent models do not differ by a factor greater than two 
during the tracking period. The error of the GEM-9 standard 
deviation model is too low. Greater differences occur for 
orbit propagation beyond the tracking period. 

• The GEM-9/MD1 and GEM-5/MD1 one-half difference 
models are similar. Perhaps as a result of the large band 
of anomalous errors discussed earlier, the orbit determina- 
tion errors predicted by these models are erratic (The ex- 
cessive errors for the 200-km-altitude , 57-deg-inclination 
orbit and the small propagation errors for Landsat-5 as dis- 
cussed in Section 4.4.) 

It may be concluded from this study that the GEM-9 standard 
deviation model and GEM-9/MD1 and GEM- 5/MDl difference model 
are all faulty and should not be used. The random-sign and 
random-phase models are difficult to justify theoretically 
and do not seem to offer any advantage over the GEM-9/SA0 
difference model. Although the GEM-9 uncorrelated model has 
not been extensively studied, it is not expected to be much 
more accurate than the GEM-9/SA0 difference model, yet the 
computational burden is much greater. It appears, from this 
limited study, that the scaled GEM-9/SA0 difference model is 
the best error model for GEM-9. A scaling factor of one- 
half is tentatively suggested before the availability of 
additional calibration. 
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SECTION 5 - ORBIT SMOOTHER PROCESS NOISE STUDY 


Because the batch method of orbit determination is a subset 
of an orbit smoother, an improved orbit determination accu- 
racy would be expected with a smoother. As discussed in 
Appendix A, the SEA program now has the capability of eval- 
uating the performance of orbit smoothers. Limited investi- 
gations using this new capability have been performed as an 
exploratory study of orbit smoothers and also to validate 
the fading memory process noise option of the SEA smoother 
capability. Two different orbit determination scenarios are 
considered: 

• The baseline TDAS navigation scenario for a 600-km- 
altitude, 28-deg inclination satellite considered 
in Section 3 

• The 34-hour arc Landsat-5 orbit determination 
scenario referred to in Section 4.4 (This investi- 
gation is undertaken in part to determine if the 
degradation of Landsat-5 orbit accuracy discussed 
in Section 2 results from geopotential errors com- 
pounded by the reduction in tracking coverage.) 

The TDAS navigation smoother results are summarized in 
Table 5-1. Two process noise options are considered: the 

linear growth option and the fading memory option. The 
linear growth option assumes that, m addition to the 
modeled evolution of the covariance matrix, there is an ad- 
ditive increase (noise) in the variances of the Cartesian 
velocity components proportional to the time elapsed since 
the last measurement. The fading memory option assumes a 
multiplicative increase of the whole covariance matrix. The 
orbit smoother does indeed improve the performance over that 
of a batch orbit determination process for this particular 
TDAS scenario. The major error sources are geopotential 
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errors and clock acceleration; the TDAS ephemens errors are 
secondary, and other error sources are negligible. Although 
the fading memory option has not been investigated m as 
much detail as the linear growth option, the former does not 
appear to offer any advantage over the latter. 

Table 5-1. Orbit Smoother Performance for the TDAS 


Navigation 

Process Noise Option 

Linear Growth Rate 
10 (m/sec) 2 /sec 
10 ^ (m/sec) 2 /sec 
10 12 . (m/sec) 2 /sec 

10 - ^ (m/sec) 2 /sec 

— 8 2 
10 (m/sec) /sec 

Fading Memory Multi- 
plicative Factor 

1.1 

1.25 

1.5 


Scenario 


Orbit 

Position Errors 

(m) 

Maxima 

Root Mean 

Square 

104 

46 


78 

41 


67 

33 


40 

24 


42 

23 



61 

24 

63 

36 

98 

38 


The Landsat-5 results are presented in Table 5-2. Tracking 
data consist of single TDRS and ground station data. The 
only error source considered is the one-half the GEM-9/SAO 
geopotential difference discussed m Section 4. Generally, 
increasing the process noise tends to decrease dynamic error 
at the expense of measurement error. Although measurement 
noise is excluded. Table 5-2 shows that an optimum noise 
level exists above which the dynamic error also increases. 
Another interesting result shown in Table 5-2 is that, with 
the addition of ground tracking data, orbit errors actually 
increase if process noise is not introduced. This is not a 
generally valid conclusion but shows that the behavior of 
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dynamic error is difficult to predict. On the other hand. 
Table 5-2 shows that, with the use of process noise, orbit 
errors are decreased when more tracking data are available. 

Table 5-2. Landsat-5 Orbit Smoother Error Caused by 
One-Half GEM-9/SAO Differences 


Trackinq Data 

Process Noise 
Linear Growth 
Rate 
2 

(m/sec) /sec 

Orbit Position 
Maximum Root 

Error (m) 
Mean Square 

Single TDRS and 

10-8 

69 

24 

Lanasat-5 ground 

10-1° 

56 

20 

tracking stations 

10" 12 

47 

24 


0 

57 

28 

As above with ad- 

10-12 

42 

17 

dition of several 

0 

69 

30 


Lanasat-4 ground 
tracking stations 

Other results not shown in Tables 5-1 and 5-2 include the 
following: 

• Maximum errors generally occur near two ends of the 
data arc where smoothing has smaller effects. 

• Orbit solutions propagated beyond the data arc are 
generally less accurate for smoothed solutions than 
for batch solutions. This is not unexpected be- 
cause the process noise introduced in an orbit 
smoother is artificial rather than physical. 
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SECTION 6 - CONCLUSIONS 


Several topics related to low-altitude satellite orbit de- 
termination have been investigated. The topics studied and 
major conclusions reached are summarized below. 

• Landsat-5 Orbit Determination using GPSPAC data-- 
Landsat-5 orbit solutions computed from GPSPAC delta- 
pseudorange (Doppler) data are good. Maximum differences 
between GPSPAC and GSTDN solutions are generally under 
70 m. Maximum differences between partially overlapping 
GPSPAC solutions are even smaller. There is a good possi- 
bility that the GPSPAC solutions are superior to the GSTDN 
solutions in accuracy. 

For orbit determination instead of real-time navigation, 
simultaneous data from four GPS satellites are not neces- 
sary. The study results indicate that approximately 3 hours 
of data from a single GPS satellite- are sufficient to re- 
solve Landsat-5 orbit and clock. A study of a randomly se- 
lected sample shows that the Landsat-5 orbit solutions based 
on individual GPS satellites agree to within 80 m, and 
typically less than 40 m. 

As with Landsat-4 GPSPAC data, some inconsistencies exist 
between the Landsat-5 GPSPAC pseudorange and delta- 
pseudorange data. Landsat-5 solutions derived from pseudo- 
range data generally differ from the GSTDN solutions by 
maximas over 100 m. Furthermore, large, GPS- independent 
range observation residuals of over 100 m are seen in 
delta-pseudorange Landsat-5 GPSPAC data solutions. Based on 
these, it may be concluded that the pseudoranges have system- 
atic errors on the order of 100 m in addition to the ex- 
pected clock errors. The causes of these systematic errors 
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have not been determined, although bad data from an individ- 
ual GPS satellite and several easily committed preprocessing 
errors have been eliminated as possible reasons. 

In connection with the GPSPAC study, the computation of 
Landsat-5 orbits from GSTDN data was undertaken. Unfor- 
tunately, there are not as many ground tracking stations for 
Landsat-5 as for Landsat-4, and the accuracy of the result- 
ing Landsat-5 orbit solutions is not as good as the corre- 
sponding Landsat-4 solutions. Investigations show that this 
degradation cannot be attributed solely to dynamic modeling 
errors accentuated by the sparsity of tracking coverage. 

• TDAS Simulation — Simulated TDAS one-way Doppler 
data were generated. The data, with errors added, were in- 
put to an extended Kalman filter in R&D GTDS for the naviga- 
tion of low-altitude Earth satellites. Comparisons with 
truth models yielded information about the navigation per- 
formance. The results showed very good agreement with those 
obtained earlier from error analysis using the SEA program. 
The agreement provides confidence in the TDAS capabilities 
of both programs and indicates that TDAS navigation perform- 
ance can be studied economically using SEA instead of the 
computationally expensive simulations. 

• Error Model for GEM-9 — Several candidate error 
models for GEM-9 were studied. Although the specification 
of an appropriate model was difficult, the use of the 
GEM-9/SAO difference model with a proper scaling factor is 
recommended. A scaling factor of one-half is tentatively 
suggested before the availability of additional calibration. 

• Orbit Smoothers--Some exploratory study of orbit 
smoothers was conducted using SEA. The fading memory proc- 
ess noise option of SEA was validated. However, this option 
does not seem to offer particular advantages over the linear 
growth option. This limited study showed that an orbit 
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smoother can indeed outperform a batch orbit determination 
process in accuracy. The following two observations may be 
made concerning orbit smoothers: (1) If accuracy is of pri- 

mary concern, a smoothed orbit near the ends of the data arc 
can be discarded because smoothing action is less at the 
ends and orbit errors are generally greater. (2) Because 
process noises are generally artificial rather than physical, 
orbit propagation beyond the data arc from a smoothed orbit 
' is generally not as accurate as the propagation from a batch 
solution. 

• Validation of the SEA Smoother/TDAS Capability--The 
SEA Smoother/TDAS capability was fully validated by test- 
ing. An inconvenience of the SEA smoother is that it will 
not output smoother error analysis results at times not co- 
incident with a measurement. It is possible to trick the 
program into outputting results by introducing dummy meas- 
urements with negligible data weights, i.e., large measure- 
ment variances. 
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APPENDIX A ~ VALIDATION AND VERIFICATION OF THE SEA SMOOTHER/ 

TDAS CAPABILITY 


The Sequential Orbit Determination Error Analysis (SEA) Pro- 
gram has been enhanced to incorporate error analysis capa- 
bilities for the following orbit determination scenarios: 

• Orbit determination using the smoother method 

• Orbit determination using the proposed Tracking and 
Data Acquisition System (TDAS) data 

These enhancements to SEA were designed and implemented in 
1984 under Task 42100. 

This appendix outlines the results of the validation and 
verification testing of these enhancements. For the con- 
venience of discussion, the enhanced SEA program will be 
referred to as SEA Version 4.1 to distinguish it from the 
unenhanced SEA Version 3.1. 

The testing to validate and verify the SEA enhancements was 
divided into three phases: 

• Phase 1 tests ensured that the enhancements have 
not corrupted the integrity of the original capabilities of 
SEA Version 3.1. 

• Phase 2 tests compared the SEA Version 4.1 error 
analysis results with those of the independent Orbit Analy- 
sis (ORAN) Program. (ORAN is an error analysis program for 
the batch orbit determination method, which can be con- 
sidered a subset of an orbit smoother. Thus, ORAN can be 
used to verify some, but not all, of the SEA smoother capa- 
bilities. ) 

• Phase 3 tests focused on the smoother algorithm's 
numerical stability and expected behavior under controlled 
operating conditions. 
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These three test phases and their results are described in 
more detail in Sections A.l, A. 2, and A. 3, respectively. 

During the course of testing several software modifications 
or corrections were required. After they were implemented, 
the testing sequence was reinitialized. These modifications 
or corrections are described in the sections describing the 
testing phases during which they occurred. 

Finally, all keyword deck setups and associated 30 b control 
language (JCL) used in testing have been archived so that 
they may be used as prototype benchmarks for future modifi- 
cations to SEA. For reference, the listings of all deck 
setups are reproduced in Section A. 5. References to setups 
are made in Sections A.l, A. 2, and A. 3. The convention for 
labeling a particular setup is as follows: PiSjRk, where i, 

j, and k are integers associated with the phase, a given 
scenario, and a given run number, respectively, for which 
the test was performed. For brevity, the corresponding SEA 
Version 3.1 control runs used in Phase 1 of the testing are 
not listed. 

A.l PARALLEL FUNCTIONAL VERIFICATION 
A. 1.1 TESTING DESCRIPTION 

This phase of testing verified the functional integrity be- 
tween SEA Version 3.1 and SEA Version 4.1 for nonsmoother 
options. It was undertaken to ensure that the smoother and 
TDAS enhancements had not corrupted the existing SEA options 
and processing. Parallel functional verification was per- 
formed under four different tracking scenarios: 

1. Direct tracking of a low-altitude satellite from 
10 ground stations — The satellite orbit is at a 200-km alti- 
tude and 28-deg inclination. An effective drag coefficient 
is alternately a solve-for or a consider parameter m the 
orbit determination process. Tracking data are available 
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whenever the satellite is visible to a station (deck setup 
PlSlTl) . 

2. Tracking of a user satellite by six Global Posi- 
tioning System (GPS) satellites. The user orbit is at a 
400-km altitude and 57-deg inclination. Tracking data are 
available during periods of visibility. Several different 
geopotential error models are considered (deck setup P1S2T1) . 

3. Scheduled relay tracking of a user satellite 
through two Tracking and Data Relay Satellite System (TDRSS) 
relays. The user satellite is at a 600-km altitude and a 
98-deg inclination (deck setup P1S3T1) . 

4. Continuous beacon tracking of a user satellite 
through three TDAS satellites (one backside, two front- 
side) — The user satellite is at a 600-km altitude and a 
28-deg inclination (deck setup P1S4T1) . 

Each of these four tracking scenarios constituted more than 
one test run. In the test runs performed for this phase, 
most nonsmoother options were tested, although not in all 
possible combinations. 

SEA Version 3.1 does not include either the smoothing or 
TDAS options; thus, the basis for comparison in this phase 
of testing was the forward-filtered results and those re- 
ports that are smoother/TDAS independent. In testing asso- 
ciated with the TDAS scenario, this comparison was slightly 
modified: the forward-filtered results of SEA Version 4.1 

were compared with a TDAS-only modification of SEA Ver- 
sion 3.1. This was to test for any interference that might 
be present when the TDAS and smoother updates were combined 
in SEA Version 4.1. In addition, the SEA Version 4.1 tests 
were run with and without requesting the smoother option to 
test for any smoother-introduced inconsistency in the 
forward-filtered results. 
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A. 1.2 TESTING RESULTS 


All tests performed using SEA Version 3.1 agreed with their 
corresponding SEA Version 4.1 runs. It may be concluded 
that the smoother/TDAS updates have not corrupted any of the 
many options preexisting in SEA. 

This agreement is demonstrated by the output listings shown 
in Figures A-l and A-2, both of wmch are associated with 
the GPS tracking scenario. Figure A-l displays the forward- 
filtered error budget at 1440 min past epoch resulting from 
the nonupdated SEA Version 3.1. Figure A-2 displays the 
corresponding output resulting from SEA Version 4.1 (run 
P1S3T1) . Although only three decimal places of accuracy are 
output, the length of propagation is long enough (1 day) so 
that any divergence would be displayed. 

A. 2 INDEPENDENT VERIFICATION 

A. 2.1 TESTING DESCRIPTION 

For the independent verification phase, the smoother error 
analysis results of SEA Version 4.1 were compared and con- 
trasted with corresponding results from the batch orbit de- 
termination error analysis program ORAN. SEA and ORAN are 
two independent programs with different capabilities; com- 
parisons were thus made under restricted conditions in which 
correspondence or near correspondence between the two pro- 
grams exists. In particular, this means that no TDAS track- 
ing and no process noise were introduced into this phase of 
testing. Also, clock errors m SEA Version 4.1 were emu- 
lated by measurement biases in ORAN. 

Under these restrictions, the same geopotential error models, 
measurement errors, and tracking systems were considered. 
However, the atmospheric drag models for the two programs 
are somewhat different. SEA uses the Harr is-Pr lester model, 
and ORAN uses the Jacchia model. 
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Figure A-2. SEA Version 4.1 Forward-Filtered Error Budget 
at 1440 Minutes. User Tracking Trhough Six 
GPS Satellites 



Because the two programs use different error analysis algo- 
rithms, some differences in the results produced were ex- 
pected even if the initial conditions and input values were 
the same. Special attention was therefore accorded the de- 
gree of disagreement between SEA and ORAN results to ensure 
that discrepancies were not larger than bounds established 
by physical considerations. 

Independent verification was performed under two different 
tracking scenarios using corresponding SEA and ORAN setups: 

1. Direct tracking of two TDRS-type spacecraft through 
a single ground tracking station--Dynamic errors considered 
are lumped geopotential modeling errors (GEM-9 15x15 - GEM-1 
8x8) and gravitational constant. Measurement bias errors 
are considered on both range and range-rate data (deck set- 
ups P2S1T1 and P2S1T10 (the 0 indicates the ORAN run) ) . 

2. Relay tracking of a user satellite through two 
TDRSS-type satellites — The user satellite is at a 400-km 
altitude and a 28-deg inclination. Dynamic and measurement 
errors are of the same type as those in the previous 
scenario. In addition, an effective drag uncertainty is 
considered (deck setups P2S2T1 through P2S2T20) . 

In this pnase of testing, the smoother results of SEA Ver- 
sion 4.1 were compared with the batch results from ORAN. In 
addition, corresponding observation measurements and propa- 
gated ephemerides between SEA Version 4.1 and ORAN served as 
additional bases for comparison. 

A. 2.2 TESTING RESULTS 

In the course of testing, it was found that the two programs 
usually agree quite well: two to three decimal digits were 

quite common for most numeric results. Close comparison be- 
tween the two programs was achieved in the direct tracking 
scenario. Figures A-3 and A-4 display the error budget 
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Figure A-3. SEA Version 4.1 Smoothed Error Budget at 
1080 Minutes. Direct Tracking on Two 
TDRS-Type Spacecraft 
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output from SEA Version 4.1 and ORAN (setups P2S1T1 and 
P2S1T10) at 1080 min after epoch. 

In the relay tracking scenario, a limitation in ORAN compli- 
cated the process of comparison between the two programs. 
ORAN does not model the TDRSs directly as relay satellites 
with given orbits; an artifice has to be used to make ORAN 
simulate the relay role of the TDRSs. ^ Naturally, greater 
disparities should be expected between the results. 

An example of the relay tracking scenario results is shown 
m Figures A-5 and A-6, in which the error budget outputs 
for 180 min past epoch are shown for SEA and ORAN, respec- 
tively (deck setups P2S2T1 and P2S2T10) . All error budgets 
agree quite well except for the errors caused by the gravity 
coefficient uncertainty and the drag uncertainty. The ORAN 
result shows considerably greater effect for the gravity 
coefficient uncertainty. This was expected because, in the 
artifice of ORAN simulation, this uncertainty also accounts 
for part of the TDRSS ephemeris errors. In other words, the 
ORAN result includes both the dynamic effect of gravity co- 
efficient uncertainty and the geometric effect of TDRSS 
ephemeris error that can be attributed to gravity coeffi- 
cient uncertainty in deriving the TDRS orbits. The drag 
uncertainty effects from the two programs are different 
because, as discussed above, two different atmospheric drag 
models were used. 

A comparison between SEA and ORAN was also performed for the 
scenario in which all three satellite orbits are solved for 
simultaneously. Figures A-7 and A-8 demonstrate the results 
of setups in which both relay and direct tracking of the 


Ieg&G Wasnington Analytical Services Center, Inc., TDRSS 
Era Orbit Determination System Review Study , B. T. Fang and 
B. P. Gibbs, December 1975. 
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Figure A-6. ORAN Batch Error Budget at 180 Minutes. Relay 
Tracking of 400-Kilometer-Altitude Spacecraft 
Through Two TDRSs 
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Figure A-7. SEA Version 4.1 Smoothed Error Budqet at Epoch. Relay 
Tracking Scenarios in Which User and Relay Orbits Are 
Solved for Simultaneously 
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Figure A-8. ORAN Batch Error Budget at Epoch. Relay Tracking 
Scenarios in V7hich User and Relay Orbits Are 
Solved for Simultaneously 
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relays are present (deck setups P2S2T2 and P2S2T20) . These 
figures display the error budgets for all three satellites 
at the epoch or starting time and show very good agreement. 

A. 3 FUNCTIONAL PERFORMANCE VALIDATION 

A. 3.1 TESTING DESCRIPTION 

As mentioned in Section A. 2.1, the SEA smoother capability 
cannot be verified in the presence of process noise by 
examining the agreements between SEA and ORAN. Thus, the 
final phase of testing concentrated on the effects of proc- 
ess noise on the smoother behavior. (See Section 5 for a 
related study.) 

If the smoother was implemented correctly, the following 
results should be observed in the tests: 

• The smoothed state variances should be smaller than 
the corresponding filtered state variances. This should be 
true irrespective of the magnitude of the process noise. 

• The smoothed state variances should increase with 
the increase of the process noise. 

• The process noise should introduce a fading memory 
feature into an estimator. Thus, the effect of dynamic 
errors is attenuated with the use of process noise whereas 
the effect of measurement noise is amplified. Furthermore, 
the filter is a one-sided fading memory estimator whereas 
the smoother is a two-sided fading memory filter. The fluc- 
tuation of errors with time is thus considerably smaller for 
a smoother than for a filter. 

This phase of testing was performed using the following 
tracking scenarios with controlled process noise parameters: 

1. Scheduled relay tracking of a user satellite through 
two TDRSs--The user satellite was at a 600-km altitude and a 
28-deg inclination. A nominal velocity variance process 
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-9 2 3 

noise of 1 x 10 m /sec and several variations were used 
(deck setup P3S1T1) . 

2. Continuous tracking of a user satellite through 
three TDAS satellites (one backside, two frontside) --The user 

satellite is at a 600-km altitude and a 28-deg inclination. 

-9 2 

A nominal velocity variance process noise of 1 x 10 m / 

3 

sec and several variations were used (deck setup P3S2T1) . 

A. 3.2 TESTING RESULTS 

During the course of this phase of testing, the following 
errors relating to the process noise in backward filtering/ 
smoothing were uncovered and corrected: 

• The process noise variance had an incorrect sign 
because the elapsed time was computed based on 
backward time differences. 

• The process noise variance was not propagated using 
the backward transition matrix. This was a subtle 
error that was intimately related to the method 
chosen to model the process noise in forward 
filtering . 

• The backward f liter/smoother executive did not in- 
clude the backward process noise added since the 
last measurement, when it computed the smoothed 
error results. 

After these corrections were made, tests performed confirmed 
that all process noise effects listed m Section A. 3.1 were 
indeed observed, thus providing confidence that the SEA 
smoother capability was working correctly. Figures A-9 
through A-12 show some sample error analysis results at 
1140 min of a 1440-min TDAS tracking arc (run corresponding 
to Scenario 2 of Section A. 3.1). The reduction of error 
from filtering (Figure A-9) to smoothing (Figure A-10) and 
the decrease in dynamic error (lumped geopotent lal) and 
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Figure A-9. SEA Version 4.1 Forward-Filtered Error Budget 
at 1140 Minutes. Nominal Velocity Variance 
Rate Process Noise of 1 x 10 -9 Meter 2 / 

Second3 
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Figure A-ll. SEA Version 4.1 Smoothed Error Budget at 

1140 Minutes. Velocity Variance Rate Process 
Noise of 1 x 10“ 6 Meter2/Second3 
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Figure A-12. SEA Version 4.1 Smoothed Error Budget at 

1140 Minutes. Velocity Variance Rate Process 
Noise of 1 x lO - ^ Me ter^/Second^ 



increase in measurement error (TDAS ephemeris error and 
others) as the process noise increases (Figure A-10, A-ll, 
and A-12) , should be noted. 

A. 4 CONCLUSIONS 

The following conclusions can be drawn from the series of 
tests performed: 

• The SEA smoother /TDAS enhancements were implemented 
successfully. The only major errors uncovered were those 
process noise errors described in Section A. 3. 2. 

• The SEA smoother /TDAS capabilility is now fully 
validated and ready for use. 

• The close agreement between the process-noise-free 
SEA smoother results and the batch ORAN results verifies 
that SEA now also has the error analysis capability for a 
batch orbit determination process. Because SEA is a well- 
structured, well-documented, and easily enhanced program, it 
can be easily modified to analyze new tracking systems and 
data types, whether these are used in a sequential or batch 
orbit determination mode. 

One of the limitations of a smoother (or a smoother error 
analysis program) is that, for discrete time measurements, 
smoothed states (or smoother errors) are only available at 
those measurement times peculiar to the tracking scenario. 
This means that the smoother results will not be output dur- 
ing data gaps between tracking passes. Smoother output from 
SEA during data gaps can be obtained by entering fictitious 
measurements at the desired output times. These arti- 
ficially introduced, fictitious measurements are given large 
standard deviations so that they do not affect the error 
analysis accuracy. 
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A. 5 LISTING OF TEST DECK SETUPS 


The following pages contain listings of representative deck 
setups used in the validation and verification of the SEA 
smoother/TDAS capabilities. 
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GRETAS. F.^QE 5S 
of poor quality 


//ZBPAP1T 1 JOB ( GC002 , 31 1H, FFF ) , 'SEA P1S1T1 ' , MSGLEVEL* (1,1), 

// MSGCLASS = A, TIME=30, CLASS =C, NOTIFY *ZBPAP 

/* JOBPARM LINES=30 
/•ROUTE PRINT PRTSS 

//•MEMBER P1S1T1 DATA IN TESTPAN LIBRARY 

//•DIRECT TRACKING OF A 200-KM 28 DEG INC S/C, W/DRAG A SOLVE-FOR 
//•RANGE AND RANGE RATE MEASUREMENTS EVERY 180 0 SECONDS 
//•OETERMIN TRACKING SCHEDULE. FOR 13 TRACKING STATIONS, FOR 24 HRS 
//•USING SMTHDU5M SEA SMOOTHER/TOAS UPDATES, SEA VER 4 1 KEYWORDS 
//•WITH SMOOTHING 

//STEP1 EXEC PGM=PAN#1 ,REGI0N=3OOK,C0ND=( 1 ,LE) 

//PANDD1 DD DSN=GCDEV . MVT SEA. PANLIB DATA .UNIT* DISK , DISP=SHR 
//PANDD2 DD DSN=&&SEAUPD . UNIT^VIO , SPACE* ( CYL ,(2,1), RLSE ) , 

// DISP* (NEW, PASS ) , DCB* ( RECFM*FB , LRECL*80 , BLKSIZE=3B00 ) 

//SYSPRINT DD SYSOUT»» 

//SYSPUNCH DD DUMMY 

//SYSIN DD DSN=ZBPAP SMTHDU5M DATA, UNIT=DISK, DISP=SHR 
// EXEC FORTRANH , PARM= ' XREF ' 

//SYSLIN DD SPACE* (CYL. (2.1)) 

// SOURCE. SYSTERM DD DUMMY 
//SOURCE . SYSPRINT DD DUMMY 

//SYSIN DD DSN*&&SEAUPD,UNIT=VIO, DISP* (OLD, DELETE) 

// EXEC LINK, PARM*'LET, LIST, MAP, SIZE=(200K, 20K) ' ,REGI0N=25OK 
//SYSLIB DD 0SN*SYS2 FORTLIB , DISP*SHR 
// DD DSN*GCOEV SEAMVS LOAD , DISP*SHR 

//SYSPRINT DD SYSOUT** 

//SYSUT1 DD SPACE * ( TRK ,(55,1,1)) 

//SYSLIN DD 
// DD * 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM=GSFC, REGI0N*4OOK 

//STEPLIB DD OSN*&&LODMOD. DISP* (OLD, DELETE) 

/ /FT01F001 DO UNIT *DISK , DCB* ( RECFM* VBS , LRECL*44 , BLKSIZE*4404 ) , 

// SPACE* ( CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

/ /FT02F001 DD UNIT *DISK , DCB* ( RECFM* VBS , LRECL = 44 , BLKSIZE = 4404 ) , 

// SPACE* (CYL, (1,1)) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT =DISK , DCB* ( RECFM* VBS , LRECL = 44 , BLKSIZE *4404 , 

// BUFNO* 1 ) , DISP* ( NEW, PASS ) , SPACE* ( CYL . ( 1 , 1 ) ) , 

// OSN*&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD DDNAME * D ATA5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYSOUT** , DCB* (RECFM*VBA, LRECL=137 , BLKSIZE=3990) , 

/ /FT06F001 DD SYSOUT** , DCB* ( RECFM«VBA, LRECL* 137 , BLKSIZE* 141 ) . 

// SPACE«(CYL, (3,1), RLSE). 

// DSN*&&PRTOUT SEA PRINTER OUTPUT 

//•FT08F001 DD DDNAME * I NPCRD SEA KEYWORD CARD INPUT 

//* INPCRD DD * 

//FT08F001 DD * SEA KEYWORD CARD INPUT 


EPOCHTIM 

0 


99 

861101 

000000 0000 




USERSATO 

1 

11 

1 

6578140. 

0 0001 

28.8 



USERSATO 

1 

12 

1 

0 

0 

0 



SATELITE 

2 

11 

0 

42166750. 

0.0004 

5 0 



SATELITE 

2 

12 

0 

358 

0 

0 



SATELITE 

3 

11 

0 

42163592 42 

0 0004 

5.0 



SATELITE 

3 

12 

0 

228 

0 

0 



STATIONS 

1 

0 

0 

323003.857 

2532329 162 

1441 37 



STATIONS 

2 

0 

0 

352029.642 

2430737 214 


0 

9 133 10D+03 

STATIONS 

3 

0 

0 

352029 642 

2430735 042 


0 

9 187 10D+03 

STATIONS 

5 

0 

0 

283029.905 

2791826.183 


-0 

558200D+02 

STATIONS 

6 

0 

0 

-075717.371 

3454022.570 


0 

528370D+03 

STATIONS 

8 

0 

0 

322105 001 

2952031 937 


-0 

337500D+02 

STATIONS 

11 

0 

0 

220734.461 

2002005 439 


0 

1139750+04 

STATIONS 

12 

0 

0 

283029.774 

2791823 853 


-0 

544500D+02 

STATIONS 

13 

0 

0 

352031.937 

2430735 561 


0 

9 127 10D+03 

STATIONS 

4 

0 

o 

402719.656 

3554953.596 


0 

808990D+03 

STATIONS 

7 

0 

0 

-330903.596 

2892001 065 


0 

7066 10D+03 

STATIONS 

9 

0 

0 

385954 656 

2830926 161 


-0 

250000D+01 

STATIONS 

10 

0 

0 

131838.265 

1444412 524 


0 

1 15950D+03 

EARTH 

15 

15 




150 0 



SPCPARAM 

1 

0 

0 

0.00272 


2 0 



EBOUTPUT 

- 

99 

2 

1440 0 

30 0 

30 0 



CLKBIAS 

1 


1 

1000000 





CLKDRIFT 

1 


1 

200 0 

0 000001 
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USERDRAG 1 101 

GRAVCOEF -1 0.25 

/* 

// DD UNIT =DISK , DSN=ZBBTF ZBEXS.GEM DATA(REC0EF2 ) , DISP = SHR 
// DD * 


SATSOLPR 

1 

-1 

1 


MEASBIAS 

2 

-1 

2000000010002. 

10.0 

MEASBIAS 

3 

-1 

13000000010013. 

001 

MEASBIAS 

2 

-1 

10000000010000. 

10 0 

MEASBIAS 

3 

-1 

10000000010000 

.001 

MEASBIAS 

2 

-1 

5000000010000. 

10 0 

MEASBIAS 

3 

-1 

12000000010000 

.001 

EPHEMERR 

1 


25.0 

23 0 40 0 

EPHERROR 



99 


CLKACCEL 

1 

-1 

.11574 


COVARANC 

1 




2 . 50+05 

2 

. 5D+05 2 5D+05 

1.0D+00 1 OD+OO 1 OD+OO 

NOISECOV 

1 

1 0 



0 



0 . 0 

1 OOOOOD-IO 1 . OOOOOD- 10 1 00000D-10 


/* 

//*FT09F001 DD SYSOUT®* , DCB= ( RECFM= VBA, LRECL® 137 , BLKSIZE=3990 ) . 

/ /FT09F001 DD SYSQUT = * . OCB= ( RECFM = VBA . LRECL® 137 , BLKSIZE® 141 ) , 

// SPACE® (CYL .(2,1), RLSE ) , 

// DSN=&&N0M0UT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//•FT10F001 DD SYSOUT®* , DCB® ( RECFM=VBA, LRECL® 137 , BLKSIZE=3990 ) , 
//FT10F001 DD SYSOUT®* , DCB=(RECFM®VBA, LRECL® 137 , BLKSIZE® 141 ) , 

// SPACE=(CYL, (2,1), RLSE), 

// DSN®&&MEAS0UT PROCESSED MEASUREMENTS PRINTER OUTPUT 

/ /FT20F001 DD DISP=(NEW, PASS ) , SPACE® (CYL, (1,1)), 

// DCB® ( RECFM=VBS , LRECL® 144 , BLKSIZE=2884 ) , UNIT=DISK . 

// DSN=&&FESUM SEQUENTIAL’ ERROR BUDGET SUMMARY 

//FT21F001 DD DISP® (NEW, PASS ), SPACE® ( CYL .( 1 , 1 )) , 

// DCB® ( RECFM»VBS . LRECL® 144 , BLKSIZE=2884 ) , UNIT=DISK , 

// DSN®&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

/ /FT27F001 DD UNIT=DISK, DCB® ( RECFM«VBS , LRECL® 196 , BLKSIZE® 1964 , 

// BUFNO® 1 ) , DISP® (NEW, PASS ) , SPACE® ( CYL , (1,1)), 

• // DSN*&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

//FT40F001 DD UNIT=DISK, DCB® ( RECFM=FB , LRECL=1688 , BLKSIZE® 18568 ) , 

// DISP® ( NEW , DELETE ) , SPACE® ( CYL , ( 2 , 1 ) , RLSE ) , 

// DSN®&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO SYSUDUMP DD DUMMY 

//*SYSUDUMP DD SYSOUT®* , SPACE® (TRK , 1 ) , 

//* DCB® (RECFM=VBA, LRECL® 137, BLKSIZE® 1922) 

//GO DATA5 DD * 

3 861101000000 0000 -1 


100020001 

3 

1800.0 

10000.0 



0 00 

1440.0 





100030001 

3 

1800.0 

10000 0 



0 00 

1440.0 





100000002 

2 

1800.0 

1.5 



0 00 

1440 0 





200000001 

2 

20 0 

1 5 



638 0 

640 0 


732.0 

733 

0 

200000001 

3 

20 0 

0 002 



638 0 

640 0 


732 0 

733. 

.0 

300000001 

2 

20.0 

1.5 



638 0 

640 0 


732.0 

733 

0 

300000001 

3 

20 O 

0 002 



638.0 

640 0 


732 0 

733 

0 

500000001 

2 

20 0 

1.5 



365 0 

369 0 


459.0 

463 

0 

552.0 

556 0 


646.0 

649 

0 

500000001 

3 

20.0 

0 002 



365 0 

369.0 


459 0 

463 

0 

552 0 

556 0 


646 0 

649 

0 

600000001 

2 

20 0 

1 5 



665 0 

668 0 
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600000001 

3 

20 

0 


0 002 



665 0 

668 0 







700000001 

2 

600 0 


1 5 



1200 0 

1440.0 







700000001 

3 

600 0 


0.002 



0.00 

240 0 







800000001 

2 

20 . 

.0 


1 5 



369 0 

373.0 



462 

0 

466 

0 

556 0 

559 0 







800000001 

3 

20 . 

.0 


0 002 



369 0 

373.0 



462 

.0 

466 

0 

556 0 

559 0 







1100000001 

2 

20 

0 


1.5 



627 0 

631.0 



721 

.0 

724 

0 

815.0 

817 0 







1100000001 

3 

20 

0 


0 002 



627.0 

631 0 



721 

0 

724 

0 

815 0 

817 0 







1200000001 

2 

20 

0 


1 5 



365 0 

369 0 



459 

0 

463 

0 

552 0 

556 0 



646 

0 

649 

0 

1200000001 

3 

20 

0 


0 002 



365 0 

369 0 



459 

0 

463 

0 

552.0 

556 0 



646 

0 

649 

0 

1300000001 

2 

20 

0 


1.5 



638.0 

640 0 



732 

0 

733 

0 

1300000001 

3 

20 

0 


0 002 



638 0 

640 0 



732 

0 

733 

0 


0 

/* 

//NTSO EXEC PGH*NOTI FY , COND= EVEN 

// 
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.//ZBPAP2T1 JOB ( GC002 , 31 1H, FFF ) , 'SEA P1S2T1 ' .TIME =30, 

// MSGCLASS=X , MSGLEVEL 5 ( 1,1) . NOTIFY=ZBPAP . CLASS=C 

/* JOBPARM LINES=30 
//•ROUTE PRINT RMT6 

//•MEMBER P 1 S2T 1 . DAT A IN TESTPAN LIBRARY 
//*6 GPS'S TRACKING A 400-KM 57 DEG INC S/C, 

//•RANGE AND RANGE RATE MEASUREMENTS EVERY 300 0 SECONDS 
//•USING SMTHDU5M SEA SMOOTHER/TDAS UPDATES, SEA VER 4 1 KEYWORDS 
//•WITH SMOOTHING 

//STEP1 EXEC PGM=PAN#1 , REGION=300K , COND= ( 1 , LE ) 

//PANDD1 DD DSN=GCDEV . MVT . SEA PANLIB DATA , UNIT= DISK , DISP=SHR 
//PANDD2 DD DSN=&&SEAUPD , UNIT=VIO , SPACE 5 ( CYL , (2,1) ,RLSE) , 

// DISP=( NEW. PASS) ,DCB=(RECFM=FB,LRECL=80,BLKSIZE=3600) 

//SYSPRINT DD SYSOUT 55 
//SYSPUNCH DD DUMMY 

//SYSIN DD DSN=ZBPAP . SMTHDU5M DAT A , UNIT =OISK , DISP=SHR 
// EXEC FORTRANH. PARM='XREF' 

//SYSLIN DD SPACE 5 ( CYL ,(2,1)) 

//SOURCE. SYSTERM DD DUMMY 
//SOURCE. SYSPRINT DD DUMMY 

//SYSIN DD DSN 5 &&SEAUPO , UNIT =VIO , DISP= ( OLD , DELETE ) 

// EXEC LINK. PARM 5 ' LET, LIST. MAP, SIZE 5 (200K.20K) ' ,REGI0N=250K 

//SYSLIB DD DSN=SYS2 . FORTLIB , OISP=SHR 

// DD DSN 5 GCDEV.SEAMVS.LOAD.DISP 5 SHR 

//SYSPRINT DD DUMMY 

//SYSUT1 DD SPACE=(TRK, (55, 1 . 1 ) ) 

//SYSLIN DD 
// DD • 

INCLUDE SYSLIB(SEA) 

ENTRY MAIN 

//GO EXEC PGM 5 GSFC,REGI0N 5 400K 

//STEPLIB DD DSN 5 &&LODMOD,DISP 5 ( OLD, DELETE) 

//FT01F001 DD UNIT s DISK , DCS 5 ( RECFM=VBS , LRECL S 44 , BLKSIZE=4404 ) , 

// SPACE 5 (CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT 5 DISK,DCB=(RECFM 5 VBS,LRECL 5 44,BLKSIZE=4404) , 

// SPACE 5 (CYL, (1,1)) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT 5 DISK,DCB 5 (RECFM 5 VBS,LRECL 5 44,BLKSIZE 5 4404, 

// BUFNO 5 1 ) , DISP 5 ( NEW , PASS ) , SPACE 5 ( CYL , ( 1 , 1 ) ) , 

// DSN 5 &&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DO DDNAME 5 DAT A5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYS0UT 5 »,DCB 5 (RECFM 5 VBA,LRECL 5 137,BLKSIZE=3990) , 
//FT0SF001 DD SYSOUT 55 ,DCB=(RECFM 5 VBA, LRECL 5 137 , BLKSIZE 5 141 ) , 

// SPACE 5 (CYL,(3,1),RLSE). 

// DSN=&&PRTOUT SEA PRINTER OUTPUT 

//•FT08F001 DD DDNAME 5 INPCRD SEA KEYWORD CARD INPUT 

//•INPCRD DD • 

//FT08F001 DD • SEA KEYWORD CARD INPUT 


EPOCHTIM 

2 

1 

-9 

880101 . 

000000 . 0000 


USERSATO 

1 

11 

1 

6978140. 

0 . 0005 

57.0 

USERSATO 

1 

12 

1 

0 . 

0 . 

0 

SATELITE 

2 

11 


26560123.0 

.001 

63.0 

SATELITE 

2 

12 


120.0 

0.0 

100 0 

SATELITE 

3 

11 


26560123.0 

.001 

63.0 

SATELITE 

3 

12 


120.0 

0 0 

140.0 

SATELITE 

4 

11 


26560123.0 

.001 

63.0 

SATELITE 

4 

12 


120.0 

0 0 

180 0 

SATELITE 

5 

11 


26560123.0 

.001 

63 0 

SATELITE 

5 

12 


240.0 

0 0 

60 0 

SATELITE 

6 

11 


26560123 0 

001 

63.0 

SATELITE 

6 

12 


240.0 

0 0 

100 0 

SATELITE 

7 

11 


26560123.0 

001 

63 0 

SATELITE 

7 

12 


240.0 

0.0 

140.0 

EARTH 

4 

4 

0 

.0 



SPCPARAM 

1 

0 

0 

.0117618 

.0 

2 0 

GRAVCOEF 

0 

0 

-1 

1.0 

.0 

0 

/» 







// DD OSN=ZBBTF . ZBEXS . GEM . DATA( REC0EF4 ) , DISP=SHR , UNIT=DISK 


// DD » 







USERORAG 

1 

0 

-1 

0.25 



CLKBIAS 

1 

0 

1 

10000. 



CLKDRIFT 

1 

0 

1 

10000 



CLK ACCEL 

1 

0 

-1 

.0001 



MEASBIAS 

2 


-1 

200010000. 

5.0 
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MEASBIAS 

2 


-1 

300010000. 

5 0 


MEASBIAS 

2 


-1 

400010000. 

5 0 


MEASBIAS 

2 


-1 

500010000 

5 0 


MEASBIAS 

2 


- 1 

600010000 

5.0 


MEASBIAS 

2 


-1 

700010000 

5 0 


MEASBIAS 

3 


-1 

200010000 0 

.001 


MEASBIAS 

3 


- 1 

300010000 0 

.001 


MEASBIAS 

3 


- 1 

400010000 0 

001 


MEASBIAS 

3 


- 1 

500010000.0 

.001 


MEASBIAS 

3 


- 1 

600010000.0 

001 


MEASBIAS 

3 


-1 

700010000.0 

.001 


EPHEMERR 

1 

0 

0 

10 0 

10.0 

20 0 

EPHERROR 

0 

0 

0 

99.0 

.0 

0 

EBOUTPUT 

0 

99 

2 

1440.0 

30 0 

30 0 

COVARANC 

1 

0 

0 




100000.0 


100000 

0 100000 0 

1.0 

10 10 

NOISECOV 

1 

1 

0 




0 




0. 0 

1 00000D-10 

1 . OOOOOD- 10 1 . OOOOOD- 10 

/* 







/* 








//*FT09F001 DD SYSOUT = * , DCB= ( RECFM=VBA, LRECL=137 , BLKSIZE=3990 ) , 

/ /FT09F001 OD SYSOUT=* , DCB*(RECFM=VBA, LRECL=137 , BLKSIZE* 141 ) , 

// SPACE* (CYL, (2,1) , RLSE) , 

// DSN=&&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//‘FTIOFOOI OD SYSOUT** , DCB= ( RECFM=VBA, LRECL=137, BLKSIZE =3990 ) , 

/ /FTIOFOOI DD SYSOUT**, DCB*(RECFM*VBA,LRECL=137,BLKSIZE=141 ) , 

// SPACE =(CYL, (2,1) ,RLSE) , 

// DSN=&&MEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

//FT20F001 DD DISP* ( NEW, PASS ), SPACE* ( CYL ,( 1 , 1 )) , 

// DCB* ( RECFM*VBS , LRECL* 144 , BLKSIZE=2884 ) . UNIT*DISK , 

// DSN*&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP* (NEW . PASS ). SPACE* ( CYL ,( 1 , 1 )) , 

// DCB* ( RECFM=VBS , LRECL* 144 , BLKSIZE=2884 ) ,UNIT=DISK . 

// DSN*&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DD UNIT=0ISK,DCB*(RECFM=VBS,LRECL*196,BLKSIZE=1964, 

// BUFNO* 1 ) , DISP* ( NEW . PASS ) . SPACE* ( CYL , ( 1 , 1 ) ) , 

// . DSN*&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

//FT40F001 DD UNIT *DISK , DCB* ( RECFM*FB , LRECL* 1688 , BLKSIZE = 1 8568 ) , 

// DISP* (NEW, DELETE), SPACE* (CYL, (2,1 ) , RLSE ) , 

// DSN*&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO SYSUDUMP DD DUMMY 

//•SYSUDUMP DO SYSOUT**, SPACE* ( TRK, 1 ) , 

//* DCB* ( RECFM=VBA, LRECL*137 , BLKSIZE* 1922 ) 

//GO. DAT AS DO * 

3 860101000000.0000 -1 


0 00 

20001 

2 

1440.0 

300.0 

1.0 

0.00 

20001 

3 

1440 0 

300 0 

0.001 

0.00 

30001 

2 

1440.0 

300.0 

1.0 

0.00 

30001 

3 

1440.0 

300.0 

0.001 

0.60 

40001 

2 

1440.0 

300.0 

1.0 

0.00 

40001 

3 

1440 0 

300 0 

0.001 

0.00 

50001 

2 

1440.0 

300.0 

1 0 

0.00 

50001 

3 

1440.0 

300.0 

0.001 

0.00 

60001 

2 

1440.0 

300 O 

1 0 

0 00 

60001 

3 

1440.0 

300.0 

0 001 

0 00 

70001 

2 

1440 0 

300.0 

1.0 

0.00 

70001 

3 

1440.0 

300.0 

0.001 


/* 
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//NTSO EXEC PGM=NOTIFY,COND=EVEN 

// 



//ZBPAP3T 1 JOB (GC002 3 1 1H . FFF ) , ' SEA P 1 S3T 1 ' , MSGLEVEL® (1,1). 

// MSGCLASS=X , TIME®30 , CLASS =C.NOTIFY=ZBPAP 

/ * JOBPARM LINES=40 
//•ROUTE PRINT PRTSS 

//•MEMBER P1S3T1 .DATA IN TESTPAN LIBRARY 

//•TDRSS TRACKING OF A 600-KM 98 DEG INC S/C 

//•RANGE AND RANGE RATE MEASUREMENTS EVERY 180 0 SECONDS 

//•USING SMTHDU5M SEA SMOOTHER/TDAS UPDATES. SEA VER 4.1 KEYWORDS 

//•WITH SMOOTHING 

//STEP1 EXEC PGM=PAN#1 , REGI0N=300K, C0ND=( 1 , LE) 

//PANDD1 DD OSN=GCDEV MVT SEA PANLIB DATA, UNIT=DISK, DISP=SHR 
//PANDD2 DD DSN=&&SEAUPD . UNIT=VIO . SPACE® ( CYL .( 2 , 1 ), RLSE ) . 

// DISP= ( NEW, PASS ) , DCB= ( RECFM=FB , LRECL=80 . BLKSIZE=3600 ) 

//SYSPRINT DD SYSOUT®* 

//SYSPUNCH DD DUMMY 

//SYSIN DD OSN*ZBPAP SMTHDU5M DATA,UNIT=DISK,DISP=SHR 
// EXEC FORTRANH , P ARM® ' XREF ' 

//SYSLIN DD SPACE= (CYL, (2,1)) 

//SOURCE. SYSTERM DD DUMMY 
//SOURCE. SYSPRINT DD DUMMY 

//SYSIN DD DSN*&&SEAUPD , UNIT =VIO , DISP= ( OLD , DELETE ) 

// EXEC LINK, PARM®' LET, LIST. MAP, SIZE® (200K,20K)',REGION=250K 
//SYSLIB DD DSN=SYS2 FORTLIB . DISP=SHR 
// DD DSNsGCDEV.SEAMVS.LOAD,DISP=SHR 

//SYSPRINT DD SYSOUT=* 

//SYSUT1 DD SPACEs(TRK, (55, 1 , 1 ) ) 

//SYSLIN DD 
// DD • 

INCLUDE SYSLIBt SEA) 

ENTRY MAIN 

//GO EXEC PGM=GSFC , REGI0N=400K 

//STEPLIB DD DSN=&&LODMOD,DISP= (OLD, DELETE) 

/ /FT01F001 DD UNIT* DISK , DCB= ( RECFM=VBS , LRECL S 44 . BLKSIZE s 4404 ) , 

// SPACE® ( CYL ,(5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT=DISK,DCB=(RECFM®VBS,LRECL=44-,BLKSIZE®4404) , 

// SPACE® (CYL, (1,1)) MERGING FOR TRACKING DATA SCHEDULE 

/ /FT03F001 DD UNIT=DISK,DCB«(RECFM®VBS, LRECL®44,BLKSIZEs4404, 

// BUFN0®1 ),DISP=( NEW, PASS), SPACE® (CYL, (1,1)), 

// DSNs&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD DDNAME®0ATA5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYSOUT®* , DCB* ( RECFM® VBA , LRECL® 137 , BLKSIZE® 3990 ) , 
//FT06F001 DD SYSOUT®* , DCB® ( RECFM®VBA, LRECL® 137 , BLKSIZE® 141 ) , 

// SPACE® ( CYL ,(3,1), RLSE ) , 

// DSN=&&PRTOUT SEA PRINTER OUTPUT 

//•FT08F001 DD DDNAME ® I NPCRD SEA KEYWORD CARD INPUT 

//•INPCRD DD * 

//FT08F001 DD * SEA KEYWORD CARD INPUT 


EPOCHTIM 

1 


99 

851101. 


000000 . 0000 




USERSATO 

1 

11 

1 

6978140. 


0.0001 

98 8 



USERSATO 

1 

12 

1 

0 . 


0 . 

0 



SATELITE 

2 

11 

1 

6778140. 


0.0017 

28 



SATELITE 

2 

12 

1 

160 0 


0 

0 



SATELITE 

3 

11 

0 

42166750. 

0 

0004 

5.0 



SATELITE 

3 

12 

0 

358 


0 . 

0 



SATELITE 

4 

11 

0 

42163592.42 

0 

0004 

5 0 



SATELITE 

4 

12 

0 

228 


0 

0 



SATELITE 

5 

1 1 

0 

42163592 42 

0 

0004 

5 0 



SATELITE 

5 

12 

0 

113 


0 . 

0 



STATIONS 

1 

0 

0 

323003 857 

2532329.162 

1441 37 



STATIONS 

3 

0 

0 

385953.936 


2830929 141 


0 

400000D+01 

STATIONS 

4 

0 

0 

402719 656 


3554953.596 


0 

808990D+03 

STATIONS 

6 

0 

0 

-075717 371 


3454022.570 


O 

528370D+03 

STATIONS 

7 

0 

0 

-330903.596 


2892001.065 


0 

7066 10D+03 

STATIONS 

8 

0 

0 

322105 001 


2952031 937 


-0. 

, 337500D+02 

STATIONS 

10 

0 

0 

131838 265 


1444412 524 


0 

1 15950D+03 

STATIONS 

11 

0 

0 

220734 461 


2002005 439 


0 

1 1 3975D+04 

STATIONS 

12 

0 

0 

283029.774 


2791823 853 


-0 

544500D+02 

STATIONS 

EARTH 

13 

0 

0 

352031 937 


2430735 561 

150 0 

0. 

- 9 127 10D+03 

SPCPARAM 

1 

0 

0 

0 0011765 


2 0 

2.0 



EBOUTPUT 

- 

99 

2 

1020 0 


30 0 

30.0 



MEASRMNT 

CLKBIAS 

1 


1 

1000000 


1020' 
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CLKDRIFT 1 1 200 0 0 000001 

USERDRAG 1 101 

GRAVCOEF -1 0 25 

/* 

// DO UNIT* DISK , DSN=ZBBTF . ZBEXS . GEM DATA( REC0EF2 ) , DISP=SHR 
// DD * 


SATSOLPR 

1 

-1 

. 1 







MEASBIAS 

2 

-1 

1000300010000. 

0 00001 





MEASBIAS 

2 

-1 

6000000010000 

10.0 





MEASBIAS 

3 

1 

10000000010000. 

001 





MEASBIAS 

2 

-1 

7000000020000 . 

10 0 





MEASBIAS 

3 

-1 

3000000020000 . 

001 





EPHEMERR 

1 


25.0 


23.0 


40 

.0 


EPHERROR 



99. 







CLKACCEL 

1 

-1 

.11574 







COVARANC 

1 









1 . 0D+05 

1 

0D+05 

1 . 0D+05 

1 OD+OO 

1 

OD+OO 

1 

OD+OO 

COVARANC 

2 

1 








1 0D+06 


0 

0 

0 . 


0 

1 

OD-06 


0 . 

1 

00+06 

0 . 

0 . 

5 

0D-01 

1 

OD-06 


0 . 


0 

1 .0D+06 

1. 0D-06 


0 

1 

OD-06 


0 


0 

1 0D-06 

5 . OD+OO 


0 

5 

0D-09 


0 . 

5 

0D-01 

0 . 

0 . 

5 

OD+OO 

5 

0D-09 

1 OD 

-06 

1 

OD-06 

1 00-06 

5 0D-09 

5 

0D-09 

5 

OD+OO 

NOISECOV 

1 

1 0 








0 


1 

0 

0 1 

OOOOOD- 10 1 

OOOOOD- 10 1. 

OOOOOD- 10 

NOISECOV 

2 

1 0 








1 OOOOOD-12 

1. 00000D-12 1, 

,000000-12 

0 


0 


0 


/* 

//*FT09F001 DD SYSOUT" , DCB 9 ( RECFM=VBA, LRECL= 137 , BLKSIZE=3990 ) , 
//FT09F001 DD SYSOUT 99 , DCB= ( RECFM=VBA, LRECL 9 137 , BLKSI2E= 141 ) , 

// SPACE 9 ( CYL ,(2,1), RLSE ) , 

// - DSN 9 &&N0M0UT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//*FT10F001 DD SYSOUT". DCB 9 (RECFM=VBA, LRECL 9 1 37, BLKSIZE=3990) , 
//FT10F001 DD SYSOUT" , DCB 9 ( RECFM=VBA, LRECL 9 137 , BLKSIZE 9 141 ) , 

// SPACE 9 (CYL, (2,1) , RLSE) , 

// DSN 9 &&MEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

/ /FT20F001 DD DISP 9 ( NEW, PASS ), SPACE 9 (CYL ,( 1 . 1 )) , 

// DCB 9 ( RECFM 9 VBS , LRECL 9 144 , BLKSIZE=2884 ) . UNIT=DISK , 

// DSN»&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP 9 ( NEW, PASS ), SPACE 9 ( CYL ,( 1 , 1 )) . 

// DCB 9 ( RECFM 9 VBS , LRECL 9 144 . BLKSIZE = 2884 ) , UNIT 9 DISK , 

// DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DD UNIT 9 DISK , DCB 9 ( RECFM 9 VBS , LRECL 9 1 96 , BLKSIZE 9 1964 , 

// BUFNO* 1), DISP 9 (NEW, PASS), SPACE 9 (CYL, (1, 1)), 

// DSN=&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

//FT40F001 DD UNIT=DISK, DCB 9 ( RECFM=FB . LRECL 9 1688 , BLKSIZE 9 18568 ) , 

// DISP 9 (NEW, DELETE) ,SPACE=(CYL, (2, 1),RLSE). 

// DSN=&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO SYSUDUMP DD DUMMY 

//*SYSUDUMP DD SYSOUT". SPACE 9 ( TRK, 1), 

//* DCB 9 ( RECFM 9 VBA, LRECL 9 1 37 , BLKSIZE 9 1 922 ) 

//G0.DATA5 DD 9 


3 861101000000.0000 -1 


100030001 

2 

180.0 

30.0 





0 00 

1440 0 







100030001 

3 

180 0 

0.01 





0 00 

1440 0 







100040001 

2 

180 0 

30 0 





0.00 

1440 0 





- 


100040001 

3 

180 0 

0.01 





0 00 

1440 0 






" 

1000400050002 2 

180.0 

30 0 





0.00 

120.0 







100000002 

2 

60.0 

1.5 





30 0 

39.0 


1122.0 

1 128 

0 


- 

1218.0 

1227.0 


1317 0 

1323 

0 



100000002 

3 

60.0 

0 002 





30 0 

39.0 


1122 0 

1128 

0 



1218 0 

1227 0 


1317 0 

1323 

0 




300000002 2 60 0 15 
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1128 0 

1134.0 


1032.0 

1038 

0 

300000002 

3 

60 0 

0 002 



1032.0 

1038.0 


1128.0 

1134 

0 

400000002 

2 

60 0 

1.5 



852.0 

861.0 





400000002 

3 

60.0 

0 002 



852 0 

861 0 





600000002 

2 

60 0 

1 5 



354.0 

363.0 


453 0 

462 

0 

552.0 

558 0 


1245 0 

1254 . 

.0 

1344.0 

1350 0 





600000002 

3 

60.0 

0 002 



354 0 

363 0 


453.0 

462 

0 

552 0 

558.0 


1245 0 

1254 

0 

1344 0 

1350.0 





700000002 

2 

60.0 

1.5 



237.0 

255 0 


339 0 

348 . 

,0 

435.0 

447 0 


534 0 

543 

0 

700000002 

3 

60 0 

0 002 



237.0 

255 0 


339 0 

348 

0 

435.0 

447 0 


534 0 

543 

0 

800000002 

2 

60 0 

1 5 



936.0 

942 0 


1032 0 

1041 

0 

1131 0 

1140.0 


1230 0 

1236 

0 

800000002 

3 

60.0 

0.002 



936 0 

942 0 


1032 0 

1041 

0 

1131 .0 

1140.0 


1230 0 

1236 

0 

900000002 

2 

180.0 

1 5 



60.0 

30.0 





1000000002 

2 

60.0 

1.5 



3 00 

12 0 


102 0 

111 

0 

201 0 

207 0 


300.0 

306 

0 

399.0 

405 0 


495.0 

504 

0 

594 0 

603.0 


1389 0 

1395 

0 

1000000002 

3 

60.0 

0 002 



3.00 

12 0 


102 0 

1 1 1 

0 

201.0 

207.0 


300.0 

306 . 

.0 

399 0 

405.0 


495 0 

504 

0 

594.0 

603.0 


1389 0 

1395 

.0 

'l 100000002 

2 

60.0 

1.5 



18 0 

27.0 


114 0 

126 

0 

213 0 

225.0 


312 0 

321 

0 

1206.0 

1212.0 


1305 0 

1311 

0 

1401 0 

1410 0 





1100000002 

3 

60.0 

0.002 



18.0 

27 0 


114 0 

126 

0 

213 0 

225 0 


312 0 

321 

0 

1206 0 

1212 0 


1305 0 

1311 

0 

1401 0 

1410.0 





1200000002 

2 

60.0 

1 5 



933 0 

939 0 


1029 0 

1038 

0 

1028.0 

1034.0 


1224 0 

1233 

0 

1323 0 

1329 O 





1200000002 

3 

60.0 

0.002 



933 0 

939.0 


1029.0 

1038 

0 

1028.0 

1034.0 


1224 0 

1233 

0 

1323 0 

1329.0 





1300000002 

2 

60 0 

1 5 



27 0 

36 0 


1218.0 

1224 

0 

1314 0 

1320.0 


1413 0 

1419 

0 

1 300000002 

3 

60 0 

0 002 



27 0 

36 0 


1218 0 

1224 

0 

1314 0 

1320.0 


1413 0 

1419 

0 
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/* 

//NTSO EXEC PGM=NOTIFY,COND=EVEN 

// 



/ /ZBPAP4T 1 JOB ( GC002 , 3 1 1H , FFF ) , ' SEA P1S4T1', 

// MSGCLASS=X , TIME =30 , MSGLEVEL= ( 1 , 1 ) ,NOTIFY‘ZBPAP,CLASS=C 

/* JOBPARM LINES=30 
//‘ROUTE PRINT RMTB 

//‘MEMBER P1S4T1 DATA IN TESTPAN LIBRARY 

//‘USING SMTHDU5M SMOOTHER/TDAS UPDATES, SEA VER 4 1 KEYWORDS 
//‘TDAS TRACKING SCENARIO, CONTINUOUS BEACON TRACKING 
//STEP1 EXEC PGM=PAN*1 , REGI0N=300K , COND= ( 1 , LE ) 

//PANDD1 DO DSN=GCDEV MVT SEA PANLIB DATA , UNIT‘DISK , DISP‘SHR 
/ /PANDD2 DD DSN=&.&SEAUPD . UNIT=VI0 . SPACE= ( CYL ,(2,1), RLSE ) , 

// DISP= ( NEW , PASS ) , DCB= ( RECFM=FB , LRECL=80 . BLKSIZE‘3600 ) 

//SYSPRINT DD SYSOUT“ 

//SYSPUNCH DD DUMMY 

//SYSIN DD DSN=ZBPAP SMTHDU5M DATA , UNIT=DISK, DISP=SHR 
// EXEC FORTRANH, PARM = 'XREF ' 

//SYSLIN DD SPACE=(CYL, (2, 1 ) ) 

//SOURCE SYSTERM DD DUMMY 
//SOURCE SYSPRINT DD DUMMY 

//SYSIN DD DSN=&&SEAUPD,UNIT=VIO , DISP= ( OLD , DELETE ) 

// EXEC LINK , PARM= ' LET, LIST , MAP , SIZE= ( 200K, 20K ) ' , REGION® 2S0K 
//SYSLIB DD DSN=SYS2 FDRTLIB , DISP=SHR 
// DD DSN=GCDEV SEAMVS LOAD , DISP=SHR 

//SYSPRINT DD SYSOUT“ 

//SYSUT1 DD SPACE=(TRK, (55, 1 , 1 ) ) 

//SYSLIN DD 
// DD ‘ 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM*GSFC , REGI0N=4OOK 

//STEPLIB DD DSN=&aLODMOD , DISP= ( OLD , DELETE ) 

/ /FT01F00 1 DD UNIT ‘DISK , DCB= ( RECFM*VBS , LRECL=44 , BLKSIZE =4404 ) , 

// SPACE=(CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT=DISK , DCB= ( RECFM=VBS , LRECL=44 , BLKSIZE=4404 ) , 

// SPACE=(CYL, ( 1 , 1 ) ) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT ‘DISK , DCB= ( RECFM=VBS , LRECL=44 , BLKSIZE = 4404 , 

// BUFNO= 1 ) ,DISP=( NEW. PASS) ,SPACE=( CYL. (1,1)), 

// DSN=&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD DDNAME=DATA5 TRACKING SCHEDULE CARD INPUT 

//‘FT06F001 DD SYSOUT“ , DCB= ( RECFM=VBA, LRECL= 137 , BLKSIZE = 3990) , 
//FT06F001 DD SYSOUT“ , DCB= ( RECFM=VBA , LRECL=137 , BLKSIZE = 141 ) , 

// SPACE=( CYL, (3,1), RLSE), 

// DSN=8l&PRT0UT SEA PRINTER OUTPUT 

//‘FT08F001 DD DDNAME=INPCRD SEA KEYWORD CARD INPUT 

//* INPCRD DD * 

//FT08F001 DD » SEA KEYWORD CARD INPUT 


EPOCHTIM 

3 


1 

840805 

120000 0000 


USERSATO 

1 

11 

1 

6681144 

0 00078765 

57 01940796 

USERSATO 

1 

12 

1 

96 4218 

360 

0 

SATELITE 

2 

1 1 

0 

42166750 

0 0004 

5 0 

SATELITE 

2 

12 

0 

360 

160 

292 

SATELITE 

3 

11 

0 

42163592 42 

0 0004 

5 0 

SATELITE 

3 

12 

0 

240 

60 

22 

SATELITE 

4 

11 

0 

42163592 42 

0 0004 

5 0 

SATELITE 

4 

12 

0 

240 

60 

292 

EARTH 

IS 

15 




150 0 

SPCPARAM 

1 

0 

0 

0 00214 


2 2 

STATIONS 

1 

0 

0 

323003 857 

2532329 162 

1441 37 

GRAVCOEF 



-1 

1 




// DD UNIT=DISK,DSN=ZBBTF ZBEXS GEM DATA( RECOEF ) . DISP=SHR 

// DD * 

USERDRAG 1 1 1. 

CLKDRIFT 1 1 1000 

CLKACCEL 1 -1 0 12 

SATSOLPR 1 -1 1. 

EPHEMERR 1 10 0 20 0 50 0 

EPHERROR 99 

EBOUTPUT 1 2 1320 10 10 

COVARANC 1 

1 0D+06 1 0D+06 1 0D+06 1 OD+OO 1 OD+OO 1 OD+OO 

NOISECOV 1 1 0 

0 0 01 OOOOOD-IO 1 OOOOOD- 10 1 OOOOOD-IO 

/* 

//‘FT09F001 DD SYSOUT“ , DCB = ( RECFM‘VBA, LRECL= 137 , BLKSIZE ‘39 90) , 
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//FT09F001 DD SYSOUT = * , DCB= ( RECFM = VBA, LRECL= 137 , BLKSIZE= 141), 

// SPACE= ( CYL ,(2,1), RLSE ) , 

// DSN=&8N0M0UT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//*FT 10F001 DD SYSOUT=» , DCB= (RECFM=VBA,LRECL=137,BLKSIZE=3990), 

//FT 10F001 DD SYSOUT=* , DCB= ( RECFM=VBA ,LRECL=137, BLKSIZE= 141), 

// SPACE=(CYL, (2, 1 ) .RLSE) , 

// DSN=&SMEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

//FT20F001 DD DISP= ( NEW , PASS ), SPACE= ( CYL ,( 1 , 1 )) , 

// DCB= ( RECFM=VBS , LRECL= 144 , BLKSIZE=2884 ) , UNIT=DISK, 

// DSN=&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP= ( NEW, PASS ). SPACE= ( CYL ,( 1 , 1 )) , 

// DCB= ( RECFM=VBS , LRECL= 144 . BLKSIZE=2884 ) . UNIT=DISK . 

// DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

/ /FT27F001 DD UNIT = DISK , DCB= ( RECFM=VBS . LRECL= 196 , BLKSIZE= 19B4 , 

// BUFNO=1) ,DISP=(NEW.PASS) .SPACE* (CYL, ( 1,1)), 

// DSN=&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

/ /FT40F001 DD UNIT=DISK . DCB= ( RECFM=FB , LRECL= 1688 , BLKSIZE= 18568 ) , 

// DISP= ( NEW, DELETE ) , SPACE = ( CYL , ( 2 , 1 ) , RLSE ) , 

// DSN=&8SMST0R SMOOTHER INFORMATION STORAGE FILE 

//GO SYSUDUMP DD DUMMY 
/ /*SYSUDUMP DD SYSOUT=* , SPACE = (TRK, 1 ) , 

//* DCB= (RECFM = VBA, LRECL= 137 , BLKSIZE= 1922 ) 

//GO DATA5 DD * 

3 840805120000 0000 -1 


000100020001 

3 180 


0 

01 


0 

25 

80 



120 

185 

215 

285 



310 

380 

415 

475 



515 

570 

610. 

665. 



700 

760 

795 

860. 



890 

960. 

985 

1055 



1090 

1160 

1190 

1285, 



1285 

000100030001 

3 180 


0 

01 


30 

80 

130 



185 

230 

285. 

335. 



380 

430. 

. 475 

525 



565. 

625. 

660 

720. 



755 

820 

860 

925 



960 

1030 

1055 

1125, 



1150 

1225 

1245. 





0001000300040001 3 180 


0 

01 


25 

30 

120. 



130 

215 

230 

310 



335 

415. 

430. 

515. 



525 

610 

625 

705 



720 

805 

820 

900 



925 

995 

1030 

1090 



1120 

1190 

1220. 

1300 



1320 


/* 

//NTSO EXEC PGM=NOTIFY , COND=EVEN 

// 
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®?I3L7ftL 

POOR 

//2BPAPS0R JOB ( GC002, 31 1H, FFF ) , 'SEA P2S1T1 ' , MSGLEVEL 2 (1,1), 

// MSGCLASS=A,TIME 2 30, CLASS 2 C, NOTIFY =ZBPAP 

/•JOBPARM LINES=30 

/•ROUTE PRINT RMT6 

//* MEMBER P2S1T1 IN TESTPAN 

II* SOLVE FOR ONLY TWO SATELLITES, PROPOGATE FOR 1 DAY 

//* WITH 15X15 GEM9 - 8X8 GEM1 MISMATCH, GM (DATA CUT-OFF AT 1440MIN) 

II* USE ZBPAP SMTHDU5M . DATA TO UPDATE THE SEA SOURCE 
//STEP1 EXEC PGM=PAN#1 , REGI0N=30OK,COND=( 1 , LE) 

//PANDD1 DD DSN=GCDEV . MVT . SEA. PANLIB . DATA, UNIT=DISK , DISP=SHR 
//PANDD2 DD DSN=&&SEAUPD , UNIT 2 VIO , SPACE 2 ( CYL, (2,1) , RLSE ) , 

// DISP= ( NEW , PASS ) , DCB= ( RECFM=FB , LRECL=80 , BLKSIZ£=3B00) 

//SYSPRINT DD SYSOUT=* 

//SYSPUNCH DD DUMMY 

//SYSIN DD DSN=ZBPAP.SMTHDU5M.DATA,UNIT 2 OISK,DISP 2 SHR 
// EXEC FORTRANH , PARM= ' XREF ' .TERM 2 '*' 

//SYSLIN DD SPACE 2 (CYL, (2,1)) 

//SYSPRINT DD DUMMY 

//SYSIN DD DSN 2 &&SEAUPD . UNIT=VIO , DI5P 2 ( OLD , DELETE ) 

// EXEC LINK, PARM 2 'LET , LIST, MAP, SIZE 2 ( 200K, 20K) ' , REGI0N=25OK, 

// NBLK 2 100 

//SYSLIB DD DSN 2 SYS2.F0RTLIB,DISP=SHR 
// DD DSN 2 GCDEV.SEAMVS.LOAD,DISP 2 SHR 

//SYSPRINT DD SYSOUT 2 * 

//SYSUT1 DD SPACE 2 (TRK, (55,1,1)) 

//SYSLIN DD 
// DD * 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM 2 GSFC,REGI0N 2 500K 

//GO. STEPLIB DD 0SN 2 &&L0DM0D,DISP=( OLD, DELETE) 

//FT01F001 DD UNIT s DISK , DCB 2 ( RECFM 2 VBS , LRECL=44 , BLKSIZE=4404 ) , 

// SPACE 2 (CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

/ /FT02F001 DD UNIT 2 DISK , DCB 2 ( RECFM s VBS , LRECL=44 , BLKSIZE s 4404 ) , 

// SPACE 2 ( CYL, ( 1 , 1 ) ) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT 2 DISK, DCB 2 ( RECFM»VBS , LRECL»44 , BLKSIZE«4404 , 

// BUFNO 2 1 ) , DISP 2 ( NEW , PASS ) , SPACE 2 ( CYL , ( 1 , 1 ) ) , 

// DSN 2 &&SORTRK SORTED TRACKING DATA SCHEDULE 

/ /FT05F001 DD 0DNAME 2 DATA5 TRACKING SCHEDULE CARD INPUT 

//»FT08F001 DD SYSOUT 22 , DCB 2 ( RECFM 2 VBA , LRECL 3 1 37 , BLKSIZE 2 1 922 ) , 
//FT08F001 OD SYSOUT 22 , DCB 2 ( RECFM=VBA, LRECL=137 , BLKSIZE 2 14 1 ) , 

// SPACE«(CYL,( 1,1), RLSE), 

// DSN 2 &&PRTOUT SEA PRINTER OUTPUT 

//•INPCRD DD 2 

//*FT08F001 OD DDNAME 2 INPCRD SEA KEYWORD CARD INPUT 

//FT08F001 DD 2 SEA KEYWORD CARD INPUT 


EPOCHTIM 

1 


1 

800301 . 

000000 0000 


USERSATO 

1 

11 

1 

42166663.0 

0.2499 

D-7 7.0 

USERSATO 

1 

12 

1 

319. 

0 . 

158.92521261 

SATELITE 

2 

11 

1 

42166663.0 

0. 17865865 

D-7 7. 

SATELITE 

2 

12 

1 

189. 

0 . 

158.9521261 


EARTH 88 40.0 

GRAVCOEF -1 0.25 

/* 

// DD UNIT=DISK,DSN 2 ZBBTF.ZBEXS.GEM.DATA(RECOEF) ,disp=shr 

// DD 2 

SPCPARAM 1 0 0 0.036 1 5 

SPCPARAM 200 0.036 1.5 

MEASBIAS 2 -1 0001000100000001. 5.0 

MEASBIAS 2 -1 0001000200000001 5 0 

MEASBIAS 3 -1 0001000100000001. 0.0005 

MEASBIAS 3 -1 0001000200000001 0.0005 

STATIONS 1 0 0 323003.857 2532329.162 1441 37 

EBOUTPUT 1 2 1440. 180. 30.0 

COVARANC 1 

1.0D+06 1.0D+06 1 . 0D+06 1 OD+OO 1.0D+00 1.0D+00 

/* 

//»FT09F001 DD SYSOUT 22 , DCB 2 ( RECFM=VBA, LRECL 2 137 , BLKSIZE 2 1922 ) , 

/ /FT09F001 DD SYSOUT 22 , DCB 2 ( RECFM=VBA, LRECL 2 137 , BLKSIZE 2 141 ) , 

// SPACE=(CYL, ( 1, 1 ) .RLSE) , 

// OSN 2 &&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//•FTIOFOOI DD SYSOUT 22 , DCB 2 ( RECFM 2 VBA, LRECL 2 137, BLKSIZE 2 1922 ) , 
//FT10F001 DD SYSOUT 22 , DCB 2 ( RECFM 2 VBA, LRECL 2 137 , BLKSIZE 2 141 ) , 


PAGE is 
® !! *Lrry 
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// SPACE =( CYL, (1,1), RLSE ) , 

// DSN = &&ME ASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

//FT20F001 DD DISP = ( NEW, PASS ), SPACE* ( CYL ,( 1 , 1 )) , 

// DCB= ( RECFM-VBS , LRECL* 144 , BLKSIZE=2884 ) , UNIT=OISK . 

// DSN=&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP=( NEW, PASS) ,SPACE=(CYL, ( 1 , 1 ) ) , 

// DCB=( RECFM=VBS , LRECL* 144 , BLKSIZE=2884 ) , UNIT=DISK , 

// DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DD UNIT=DISK,DCB=(RECFM=VBS,LRECL=196,BLKSIZE=1964, 

// BUFNO=1 ) ,DISP=( NEW, PASS) ,SPACE=(CYL, (1.1)). 

// DSN = && AKKHI N VISIBLE TRACKING DATA SCHEDULE 

/ /FT40F001 DD UNIT=DISK. DCB= ( RECFM=FB , LRECL= 1688 . BLKSIZE* 18568 ) , 

// DISP* (NEW, DELETE) .SPACE = ( CYL. (2, 1 ) . RLSE ) . 

// DSN=&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO SYSUDUMP DD DUMMY 

//•SYSUDUMP DD SYSOUT - * , SPACE = ( TRK , 1 ) , 

/ / * DCB= ( RECFM= VBA , LRECL =137,BLKSIZE=1922) 

//GO DATA5 DD « 

3 800301000000.0000 0 


00010001 

0001 2 

10 

30.0 


0 0 

1 0 


120 0 

121.0 

240.0 

241.0 


360.0 

361.0 

480.0 

481 0 


600 0 

601.0 

720.0 

721 .0 


840.0 

841.0 

960.0 

961.0 


1080 

1081 . 

1200. 

1201 


1320. 

1321 

1440. 

1440. 




00010001 

0001 3 

10. 

0.003 


0,0 

1 .0 


120 0 

121.0 

240 0 

241.0 


360 0 

361.0 

480.0 

481.0 


600.0 

601.0 

720.0 

721 0 


840.0 

841.0 

960.0 

961 0 


1080. 

1081 . 

1200 

1201 . 


1320. 

1321. 

1440 

1440. 




00010002 

'0001 2 

10. 

30.0 


1.0 

2.0 


121.0 

122.0 

241 0 

242 0 


361.0 

362.0 

481.0 

482.0 


601.0 

602.0 

721.0 

722.0 


841.0 

842.0 

961 0 

962.0 


1081. 

1082. 

1201. 

1202 


1321. 

1322. 

00010002 

0001 3 

10 

0.003 


1 0 

2.0 


121.0 

122.0 

241.0 

242.0 


361.0 

362.0 

481.0 

482.0 


601 0 

602.0 

721.0 

722.0 


841 .0 

842.0 

961.0 

962.0 


1081 . 

1082. 

1201 . 

1202. 


1321. 

1322 


/* 

//NTSO EXEC PGM =NOT I FY , COND= EVEN 

// 
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//ZBPAP2S2 JOB (GC002,311H,FFF) , 'SEA P2S2T 1 ' , MSGLEVEL® (1,1), 

// HSGCLASS=X , TIME*30 , CLASS=C , NOTIFY*ZBPAP 

/•JOBPARM LINES=30 

/•ROUTE PRINT RMTB 

//* MEMBER P2S2T1 IN TESTPAN 

//* SOLVE FOR ONLY ONE OF 3 SATELLITES, PROPOGATE FOR 1 DAY 
//« USE ZBPAP . SMTHDU5M . DATA TO UPDATE THE SEA SOURCE 
//STEP1 EXEC PGM=PAN#1,REGI0N=3OOK,C0ND=(1,LE) 

//PANDD1 DD DSN=GCDEV MVT . SEA. PANLIB . DATA , UNIT “DISK , DISP*SHR 
//PANOD2 DD DSN*&&SEAUPO,UNIT*VIO, SPACE* f CYL, (2,1) ,RLSE) , 

// DISP= ( NEW , PASS ) , DCS 5 ( RECFM=FB , LRECL=80 , BLKSIZE=3800 ) 

//SYSPRINT DD SYSOUT** 

//SYSPUNCH DD DUMMY 

//SYSIN DD DSN=ZBPAP.SMTHDU5M.DATA,UNIT=DISK,DISP*SHR 
// EXEC FORTRANH,PARM='XREF' .TERM 5 '*' 

//SYSLIN DD SPACE 5 CCYL, (2,1)) 

//SYSPRINT DD DUMMY 

//SYSIN DO DSN=&&SEAUPD , UNIT =VIO , DISP* ( OLD , DELETE ) 

// EXEC LINK, PARM='LET. LIST, MAP, SIZE=(200K, 20K) ' ,REGI0N=25OK. 

// NBLK=100 

//SYSLIB DD DSN®SYS2 . FORTLIB , DISP=SHR 
// DO DSN*GCDEV.SEAMVS.LOAD,DISP=SHR 

//SYSPRINT DD SYSOUT** 

//SYSUT1 DD SPACE=(TRK. (55, 1, 1 ) ) 

//SYSLIN DD 
// DD * 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM*GSFC.REGI0N*500K 

//GO . STEPLIB DD DSN*&&LODMOD , DISP= ( OLD , DELETE ) 

/ /FTOIFOOI DD UNIT=DISK, DCB= ( RECFM=VBS , LRECL*44 , BLKSIZE=4404 ) , 

// SPACE® (CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT =OISK , DCS® ( RECFM=VBS , LRECL=44 , BLKSIZE =4404 ) , 

// SPACE = ( CYL , ( 1 , 1 ) ) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT®DISK,DCB>(RECFM>VBS,LRECL®44,BLKSIZE®4404, 

// BUFNO® 1 ) , DISP=(NEW, PASS) , SPACE® (CYL, (1,1)), 

// OSN®&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD D0NAME>DATA5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYSOUT®® , DCB®(RECFM«VBA, LRECL® 137 , BLKSIZE* 1922 ) , 
//FT06F001 DO SYSOUT®* , DCB=(RECFM=VBA, LRECL* 137 .BLKSIZE® 141 ) , 

// SPACE® ( CYL ,(1,1), RLSE ) , 


// 


DSN* 

&&PRTOUT SEA PRINTER OUTPUT 


//*FT08F001 

DD 

DDNAME * INPCRD 

SEA KEYWORD CARD INPUT 


//*INPCRD 

DC 

1 * 





/ /FT08F001 

1 DO * 


SEA KEYWORD CARD INPUT 


EPOCHTIM 

1 


1 

800301 . 

000000.0000 


USERSATO 

1 

ii 

1 

6778140. 

0.0017 

28. 

USERSATO 

1 

12 

1 

0 . 

0 . 

0 

SATELITE 

2 

11 


42166883.0 

0.2499 D- 

7 7.0 

SATELITE 

2 

12 


319 

0 . 

158.92521261 

SATELITE 

3 

11 


42166683.0 

0.17885865 D- 

7 7. 

SATELITE 

3 

12 


189. 

0 . 

158 9521261 

EARTH 

8 

8 





SPCPARAM 

1 

0 

0 

0.0011765 

1.5 

2 0 

SPCPARAM 

2 

0 

0 

0.036 

0 . OOOOOOOOOOOOOOOO 1 

0 . 00000000000000 1 

SPCPARAM 

3 

0 

0 

0.036 

0 . 00000000000000001 

0 . 00000000000000 1 

STATIONS 

1 

0 

0 

323003 . 857 

2532329.162 

1441.37 

USERDRAG 



-1 

0.25 



GRAVCOEF 

/* 



-1 

0.25 




// DD UNIT *DISK , DSN*ZBBTF . ZBEXS . GEM . DATA( RECOEF ) , DISP*SHR 
// DD » 


MEASBIAS 

2 

-1 

0001000200010000. 

5 0 


MEASBIAS 

2 

-1 

OOO 1 0003000 1 oooo . 

5.0 


MEASBIAS 

3 

-1 

000 1 0002000 1 oooo . 

0.0005 


MEASBIAS 

3 

-1 

0001000300010000. 

0.0005 


EBOUTPUT 

COVARANC 

1 

1 

2 

1440.0 

180.0 

30.0 

1.0D+05 

1 

.00+05 1.0D+05 

1.0D+00 

1.0D+00 1.0D+00 


/* 

//*FT09F001 DD SYSOUT** , DCB* (RECFM=VBA, LRECL® 137 , BLKSIZE* 1922 ) , 
//FT09F001 DD SYSOUT*® , DCB*( RECFM=VBA, LRECL* 137 , BLKSIZE* 141 ) , 

// SPACE=(CYL,( 1,1), RLSE), 
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// OSN-&&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//•FTIOFOOI DO SYSOUT 3 * , DCB 3 ( RECFM=VBA , LRECL 3 137 , BLKSIZE 3 1922 ) , 
//FTIOFOOI DD SYSOUT 3 * , DCB=( RECFM*VBA, LRECL 3 137 , BLKSIZE= 141 ) , 

// SPACE 3 (CYL,(1,1),RLSE). 

// DSN=a&MEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

/ /FT20F001 DD DISP 3 ( NEW , PASS) , SPACE =■( CYL. ( 1 . 1 )) , 

// DCB= ( RECFM 3 VBS , LRECL 3 144 , BLKSIZE=2884 ) . UNIT=DISK , 

// DSN=&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP 3 (NEW, PASS) .SPACE 3 (CYL. ( 1 , 1 )) , 

// DCB 3 ( RECFM 3 VBS , LRECL 3 144, BLKSIZE=2884 ) , UNIT=DISK, 

// DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DD UNIT 3 DISK,DCB 3 (RECFM 3 VBS,LRECL 3 19B,BLKSIZE 3 18B4. 

// BUFNO 3 1 ) . DI SP 3 ( NEW , PASS ) . SPACE 3 ( CYL , ( 1 , 1 ) ) , 

// DSN 3 &&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

//FT40F001 DD UNIT 3 DISK, DCB 3 ( RECFM=FB , LRECL 3 1688 , BLKSIZE 3 1 85B8 ) , 

// DISP 3 ( NEW , DELETE ) . SPACE 3 ( CYL . ( 2 . 1 ) . RLSE ) . 

// DSN=&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO . SYSUDUMP DD DUMMY 

// 3 SYSUDUMP DD SYSOUT 33 ,SPACE 3 (TRK,1), 

//* DCB=(RECFM 3 VBA, LRECL 3 137, BLKSIZE 3 1922 ) 

//G0.DATA5 DD 3 

3 800301000000.0000 0 


000100020001 

10. 

2 

360. 

600. 

30. 

000100020001 

10 

3 

360. 

600. 

0.3 

000100030001 

10. 

2 

360. 

600. 

30. 

000100030001 

10. 

3 

360. 

800. 

0.3 


/» 

//NTSO EXEC PGM=NOTIFY , C0ND=EVEN 

// 


A- 3 8 



//ZBCPYSOT JOB ( GJ002 , 31 1H, FFF ) , 'SEA P2S2T2 ' , MSGLEVEL® (1,1), 

// MSGCLASS=X, TIME=30 , CLASS=C , NOTIFY=ZBCPY 

/* JOBPARM LINES=30 

/•ROUTE PRINT PRTSS 

//* MEMBER P2S2T2 IN TESTPAN 

//* SOLVE FOR ALL THREE SATELLITES. PROPOGATE FOR 1 DAY 
//* *««* USING EQUAL WEITHTS FOR A DATA TYPE- REGARDLESS OF S/C 
//* USE ZBCPY . SMTHDU5M DATA TO UPDATE THE SEA SOURCE 
//STEP1 EXEC PGM*PAN#1 ,REGI0N*3OOK,C0ND=( 1 ,LE) 

//PANDD1 OD DSN=GCDEV . MVT . SEA. PANLIB . DATA . UNIT=DISK , DISP=SHR 
//PANDD2 DD DSN*&&SEAUPD, UNIT*VIO , SPACE® (CYL , (2,1) ,RLSE) . 

// DISP= ( NEW , PASS ) , DCB= ( RECFM=FB , LRECL=80 , BLKSIZE=3600 ) 

//SYSPRINT DD SYSOUT®* 

//SYSPUNCH DD DUMMY 

//SYSIN DD OSN»ZBCPY SMTHDU5M.DATA.UNIT=DISK,DISP=SHR 
// EXEC FORTRANH , PARM= ' XREF ' , TERM® ' * ' 

//SYSLIN DO SPACE=(CYL, (2,1)) 

//SYSPRINT DD SYSOUT®* 

//SYSIN DD DSN*&8SEAUPD,UNIT=VI0, DISP® (OLD, DELETE) 

// EXEC LINK , PARM= ' LET , LIST , MAP , SIZE = ( 200K , 20K ) ' , REGI0N=250K , 

// NBLK= 100 

//SYSLIB DD DSN=SYS2 . FORTLIB , DISP=SHR 
// DD DSN*GCDEV.SEAMVS.LOAD,DISP=SHR 

//SYSPRINT DD SYSOUT®* 

//SYSUT1 DD SPACE® (TRK, (55,1,1)) 

//SYSLIN DD 
// DD * 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM*GSFC , REGI0N=500K 

//GO . STEPLIB DD DSN®&&LODMOD , DISP® ( OLD , DELETE ) 

//FT01F001 OD UNIT*DISK,DCB»(RECFM*VBS,LRECL=44.BLKSIZE*4404) , 

// SPACE® ( CYL ,(5,1)) SORTING FOR TRACKING DATA SCHEDULE 

/ /FT02F001 DD UNIT ®DISK , DCB® ( RECFM® VBS , LRECL=44 , BLKSIZE=4404 ) , 

// SPACE® ( CYL ,(1,1)) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT ®DISK , DCB® ( RECFM® VBS , LRECL®44 , BLKSIZE®4404 , 

// BUFNO® 1 ) , DISP® ( NEW . PASS ) . SPACE® ( CYL , ( 1 . 1 ) ) , 

// DSN>&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD DDNAME®DATA5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYSOUT®* ,DCB«(RECFM«VBA.LRECL=137,BLKSIZE=1922) , 
//FT06F001 DD SYSOUT®* , DCB® ( RECFM-VBA, LRECL® 137 , BLKSIZE® 141 ) , 

// SPACE® ( CYL ,(1,1), RLSE ) , 

// DSN=&&PRTOUT SEA PRINTER OUTPUT 

//FT08F001 DD DDNAME * INPCRD SEA KEYWORD CARO INPUT 

//*FT09F001 DD SYSOUT®* , DCB® ( RECFM=VBA, LRECL® 137 , BLKSIZE® 1922 ) , 
//FT09F001 DD SYSOUT®* , DCB® ( RECFM*VBA, LRECL® 137 , BLKSIZE* 141 ) , 

// SPACE® (CYL, (1,1), RLSE), 

// DSN*&&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//*FT10F001 OD SYSOUT®* , DCB® (RECFM=VBA. LRECL® 137, BLKSIZE* 1922) , 
//FT10F001 DD SYSOUT®*, DCB® (RECFM»VBA, LRECL® 137, BLKSIZE® 141), 

// SPACE»(CYL,< 1,1), RLSE), 

// DSNs&SlMEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

/ /FT20F001 DD DISP* (NEW, PASS) .SPACE® (CYL, ( 1 . 1 )) , 

// DCB* ( RECFM*VBS , LRECL® 144 . BLKSIZE=2884 ) . UNIT*DISK , 

// DSN*&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

/ /FT21F001 OD DISP* ( NEW, PASS ) , SPACE* ( CYL .(1,1)), 

// DCB* ( RECFM=VBS , LRECL* 144 , BLKSIZE® 2884 ) , UNIT = DISK , 

// DSN*&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DO UNIT=DISK,DCB=(RECFM=VBS,LRECL=196.BLKSIZE=1964, 

// BUFNO® 1), DISP® (NEW. PASS). SPACE* (CYL, (1,1)). 

// OSN=&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

/ /FT40F001 OD UNIT*DISK , DCB® ( RECFM=FB , LRECL® 1688 , BLKSIZE® 18568 ) . 

// DISP® ( NEW , DELETE ) , SPACE® ( CYL , ( 2 . 1 ) , RLSE ) . 

// DSN*&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO . SYSUDUMP DD DUMMY 
/ /*SYSUDUMP DD SYSOUT®* , SPACE® (TRK, 1 ) , 


//* 

//INPCRD 

DD 

* 

DCB 

=( RECFM® VBA, LRECL® 

137, BLKSIZE* 1922) 


EPOCHTIM 

1 


1 

800301 . 

000000.0000 


USERSATO 

1 

11 

1 

6778140. 

0.0017 

28. 

USERSATO 

1 

12 

1 

0 

0 

0. 

SATELITE 

2 

11 

1 

42166663.0 

0 2499 

D-7 7 0 

SATELITE 

2 

12 

1 

319, 

0. 

158.92521261 
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SATELITE 

3 

11 

1 

42166663.0 

0 17865865 

D-7 7 

SATELITE 

3 

12 

1 

189. 

0 . 

158 

EARTH 

8 

8 




40 

SPCPARAM 

1 

0 

0 

0.0011765 

2.0 

2.0 

SPCPARAM 

2 

0 

0 

0.036 

1.5 


SPCPARAM 

3 

0 

0 

0.036 

1.5 


STATIONS 

1 

0 

0 

323003 857 

2532329.162 

1441 

MEASBIAS 

2 


-1 

000 1 0002000 1 0000 . 

5.0 


MEASB1AS 

2 


-1 

0001000300010000. 

5 0 


MEASBIAS 

3 


-1 

0001000200010000 

0.005 


MEASBIAS 

3 


-1 

0001000300010000 

0.005 


EBOUTPUT 


0 

2 

1440. 

180. 

30.0 

COVARANC 

1 






1 . 0D+05 


1 

. 0D+05 1.0D+05 

1 OD+OO 

1. OD+OO 


/* 

//G0.DATA5 DO * 


3 80030 1 000000 . 0000 


000100020001 

2 

600 

30 


10. 

360. 




000100020001 

3 

600. 

0 3 


10. 

360. 




000100030001 

2 

600. 

30. 


10. 

360. 




000100030001 

3 

600. 

0.3 


10. 

360. 




00010002 

0001 2 

10. 

30.0 


0.0 

1.0 


120 0 

121 0 

240.0 

241.0 


360.0 

361 0 

480.0 

481.0 


600.0 

601.0 

720.0 

721.0 


840.0 

841.0 

960.0 

961 0 


1080. 

1081 . 

1200. 

1201. 


1320. 

1321 . 

1440. 

1440. 




00010002 

0001 3 

10. 

0.3 


0.0 

1.0 


120.0 

121.0 

240.0 

241.0 


360.0 

361 .0 

480.0 

481.0 


600.0 

601.0 

720.0 

721.0 


840.0 

841.0 

960 0 

961.0 


1080. 

1081. 

1200. 

1201. 


1320. 

1321 . 

1440. 

1440. 




00010003 

0001 2 

10. 

30.0 


1.0 

2.0 


121.0 

122.0 

241.0 

242.0 


361.0 

362.0 

481.0 

482.0 


601.0 

602 0 

721.0 

722.0 


841.0 

842 0 

961.0 

962.0 


1081 . 

1082. 

1201. 

1202. 


1321 . 

1322. 

00010003 

0001 3 

10. 

0.3 


1.0 

2.0 


121.0 

122.0 

241 .0 

242.0 


361.0 

362.0 

481.0 

482.0 


601.0 

602.0 

721.0 

722.0 


841.0 

842 0 

961 0 

962 0 


1081 

1082. 

1201 

1202 


1321 . 

1322. 


/* 

//NTSO EXEC PGM=NOTIFY,CDND=EVEN 

// 


9521261 

0 

37 

I 

1 OD+OO 
1 
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//ZBPAPSWN JOB ( GC002 , 31 1H, FFF ) , 'SEA P3S1T1 ' , MSGLEV£L=( 1,1). 

// MSGCLASS=A, TIME=30 , CLASS=C , NOTIFY=ZBPAP 

/* JOBPARM LINES=30 

/•ROUTE PRINT PRTSS 

//* MEMBER P3S1T1 IN TESTPAN 

//* SOLVE FOR ONLY ONE OF 3 SATELLITES, PROPOGATE FOR 1 DAY 
//* WITH VELOCITY PROCESS NOISE OF 1 0D-09 ON STATE VELOCITIES 
//* USE ZBP AP . SMTHDU5M . DAT A TO UPDATE THE SEA SOURCE 
//* SMTHDUP3 MODIFIED TO SET VAR DELTAT WITH THE OPPOSITE SIGN 
//STEP1 EXEC PGM=fAN#1,REGI0N=300K,C0ND=(1,LE) 

//PANDD1 DD DSN=GCDEV . MVT . SEA . PANLIB . DATA, UNIT=DISK , DISP=SHR 
//PANDD2 DD DSN*&&SEAUPD , UNIT = VIO , SPACE = ( CYL ,(2,1), RLSE ) , 

// DI5P= ( NEW , PASS ) , DCB= ( RECFM=FB , LRECL=80 , BLKSIZE=3600 ) 

//SYSPRINT DD SYSOUT=* 

//SYSPUNCH OD DUMMY 

//SYSIN DD OSN=ZBPAP . SMTHDU5M DATA , UNIT=DISK, DISP=SHR 
// EXEC FORTRANH , PARM= ' XREF ' , TERM® ' * ' 

//SYSLIN DD SPACE = ( CYL , ( 2 , 1 ) ) 

//SYSPRINT OD DUMMY 

//SYSIN DD DSN*&&SE AUPD , UNIT *VIO , DISP= ( OLD , DELETE ) 

// EXEC LINK, PARM*' LET, LIST, MAP, SIZE* (200K.20K) ' , REGI0N=250K, 

// NBLK=100 

//SYSLIB DD DSN=SYS2 . FORTLIB , DISP*SHR 
// DD DSN s GCDEV . SEAMVS . LOAD , DISP=SHR 

//SYSPRINT DD SYSOUT** 

//SYSUT1 DD SP ACE = ( TRK ,(55,1,1)) 

//SYSLIN DD 
// DD * 

INCLUDE SYSLIB (SEA) 

ENTRY MAIN 

//GO EXEC PGM*GSFC,REGI0N=500K 

//GO . STEPLIB DD DSN*&&LODMOD , DISP* ( OLD , DELETE ) 

//FTOIFOOI DD UNIT*OISK , DCB= ( RECFM* VBS , LRECL=44 , BLKSIZE=4404 ) , 

// SPACE* (CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT*DISK,DCB*(RECFM*VBS.LRECL=44,BLKSIZE=4404) , 

// SPACE* ( CYL, ( 1 , 1 ) ) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT*DISK,DCB*(RECFM*VBS,LRECL=44',BLKSIZE=4404, 

/ / BUFNO* 1 ) , DISP* ( NEW , PASS ) , SPACE* ( CYL , ( 1 , 1 ) ) . 

// DSN-&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DD D0NAME*DATA5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYS0UT**,DCB*(RECFM*VBA,LRECL=137,BLKSIZE=1922) , 
//FT06F001 DD SYS0UT**,DCB*(RECFM*VBA,LRECL*137,BLKSIZE=141 ) , 

// SPACE*(CYL,( 1,1), RLSE), 

// DSN*&&PRTOUT SEA PRINTER OUTPUT 

//•FT08F001 DD DDNAME * I NPCRD SEA KEYWORD CARD INPUT 

//•INPCRD OD * 

//FT08F001 DD * SEA KEYWORD CARD INPUT 


EPOCHTIM 

1 


1 

800301 . 

000000.0000 


USERSATO 

1 

11 

1 

6778140. 

0.0017 

28 

USERSATO 

1 

12 

1 

0 . 

0 . 

0 

SATELITE 

2 

11 


42166663 . 0 

0.2499 

D-7 7 0 

SATELITE 

2 

12 


319. 

0 . 

158.92521261 

SATELITE 

3 

11 


42166663.0 

0. 17865865 

D-7 7 

SATELITE 

3 

12 


189 

0 

158 9521261 

EARTH 

8 

8 





SPCPARAM 

1 

0 

0 

0 0011765 


0 00000000001 

SPCPARAM 

2 

0 

0 

0.036 

1 5 


SPCPARAM 

3 

0 

0 

0.036 

1 .5 


STATIONS 

1 

0 

0 

323003.857 

2532329 162 

1441 37 

GRAVCOEF 

/• 



-1 

0.25 



// DD UNIT 

•DISK, 

OSN*ZBBTF . ZBEXS . GEM 

. DATA( RECQEF ) , DISP 

= SHR 

// DD * 







MEASBIAS 

2 


-1 

0001000200010000. 

5.0 


MEASBIAS 

2 


-1 

000 1 0003000 1 0000 . 

5.0 


MEASBIAS 

3 


-1 

0001000200010000. 

0.0005 


MEASBIAS 

3 


-1 

000 1 0003000 1 0000 . 

0 . 0005 


EBOUTPUT 


1 

2 

360. 

30. 

30.0 

COVARANC 

1 







1 0D+05 1.0D+05 1 0D+05 I.OD+OO 1 OD+OO I.OD+OO 

NOISECOV 1 1 0 

0. OD+OO 0. OD+OO 0. OD+OO 1.0D-09 1 0D-09 1 0D-09 

/* 
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//•FT09F001 DD SYSOUT=* , DCB= ( RECFM=VBA . LRECL 2 1 37 , BLKSIZE 2 1922 ) , 
//FT09F001 DD SYSOUT = * , DCB= ( RECFM=VBA, LRECL=137 , BLKSIZE= 141 ) , 

// SPACE= ( CYL ,(1,1), RLSE ) , 

// DSN=&&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//*FT10F001 DD SYSOUT 2 * , DCB= ( RECFM=VBA, LRECL= 137 , BLKSIZE= 1922 ) , 
//FTIOFOOI DD SYSOUT 2 * , DCB= ( RECFM=VBA, LRECL 2 137 , BLKSIZE 2 141 ) , 

SPACE= ( CYL , ( 1 , 1 ) , RLSE ) , 

DSN-&&MEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

DISP= ( NEW , PASS ) . SPACE 2 ( CYL . ( 1 , 1 ) ) , 

DCB= ( RECFM= VBS , LRECL= 144 , BLKSIZE =2884 ) , UNIT=DISK , 
DSN=&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

DISP= ( NEW , PASS ) , SPACE 2 ( CYL . ( 1 . 1 ) ) , 

DCB= ( RECFM=VBS . LRECL 2 144 , BLKSIZE=2884 ) . UNIT=DISK , 
DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

UNIT =DISK . DCB= ( RECFM=VBS . LRECL 2 196 , BLKSIZE 2 1964 , 

BUFNO 2 1 ) . DISP= ( NEW , PASS ) , SPACE 2 ( CYL , ( 1 , 1 ) ) , 

DSN=&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

UNIT=DISK , DCB= ( RECFM=FB , LRECL 2 1688 , BLKSIZE 2 1 8568 ) , 

DISP* ( NEW , DELETE ) . SPACE 2 ( CYL , ( 2 , 1 ) , RLSE ) , 

OSN=&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

DD DUMMY 

DD SYSOUT 2 *,SPACE=(TRK,1), 

DCB=(RECFM=VBA, LRECL 2 137, BLKSIZE 2 1922) 


DO 


DD 


DD 


DD 


// 

// 

//FT20F001 

// 

// 

/ /FT21F001 

// 

// 

//FT27F001 

// 

// 

//FT40F001 

// 

// 

//GO . SYSUDUMP 
//*SYSUDUMP 
//* 

//G0.DATA5 DD 
3 

000100020001 

10 . 

000100020001 

10 . 

000100030001 

10 

000100030001 

10 


80030 1 000000 . 0000 
2 600 


360 


360 


600 


600. 


360 


3 600. 


30. 


0.3 


30 


0.3 


360 


/* 

//NTSO EXEC PGM 2 NOTIFY,COND 2 EVEN 

// 
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//ZBCPYSTS JOB ( GJ002, 31 1H, FFF ) , 'SEA P3S2T1 ' ,TIME=30, 

// MSGCLASS=X,MSGLEVEL=( 1,1) ,NOTIFY=ZBCPY ,CLASS=C 

/* JOBPARM LINES=100 
/•ROUTE PRINT PRTSS 
//•MEMBER P3S2T1 IN TESTPAN 

//•BEAKON TRACKING ( FLBT ) FOR A SPACE TELESCOPE MISSION MODEL 
//•RANGE AND RANGE RATE MEASUREMENTS EVERY 180 0 SECONDS 
//•USING SMTHDU5M SEA SMOOTHER/TDAS UPDATES, SEA VER 4 1 KEYWORDS 
//•WITH SMOOTHING 

//•WITH LOW ACCCURACY TRACKING EVERY 30 MINUTES TO INCREASE THE 

//•NUMBER OF SMOOTHED ERROR BUDGETS 

//STEP1 EXEC PGM=PAN#1 , REGI0N=300K , COND= ( 1 , LE ) 

//PANDD1 DD DSN=GCDEV . MVT . SEA . PANLIB . DATA, UNIT s OISK , DISP=SHR 
//PANDD2 DD DSN 3 8&SEAUPD,UNIT 3 VIQ, SPACE 3 (CYL, (2,1) ,RLSE) , 

// DISP 3 ( NEW , PASS ) , DCB= ( RECFM=FB , LRECL=80 , BLKSIZE=3B00 ) 

//SYSPRINT DD SYSOUT=* 

//SYSPUNCH DD DUMMY 

//SYSIN DD DSN s ZBCPY . SMTHDU5M. DATA, UNIT 3 DISK, DISP=SHR 
// EXEC FORTRANH , PARM 3 'XREF' .TERM 3 ' • ' ,OUT= '* ' 

//SYSLIN DD SPACE 3 ( CYL, (2,1)) 

//SOURCE. SYSTERM DD SYSOUT 3 * 

//SOURCE. SYSPRINT DD DUMMY 

//SYSIN DD DSN=&&SEAUPO , UNIT =VIO , DISP 3 ( OLO , DELETE ) 

// EXEC LINK, PARM='LET, LIST, MAP, SIZE 3 ( 200K , 20K) ' , REGI0N=250K 
//SYSLIB DD DSN 3 SYS2.F0RTLIB,DISP 3 SHR 
// DD DSN=GCDEV . SEAMVS LOAD , DISP 3 SHR 

//SYSPRINT DD SYSOUT 3 * 

//SYSUT1 DD SPACE 3 (TRK, (55,1,1)) 

//SYSLIN DD 
// OD • 

INCLUDE SYSLIB(SEA) 

ENTRY MAIN 

// EXEC PGM 3 GSFC,REGI0N 3 400K 

//STEPLIB DO DSN 3 &&L0DM0D, DISP 3 (OLD, DELETE) 

//FT01F001 DD UNIT 3 OISK, DCB 3 ( RECFM=VBS , LRECL 3 44 , BLKSIZE 3 4404 ) , 

// SPACE 3 (CYL, (5,1)) SORTING FOR TRACKING DATA SCHEDULE 

//FT02F001 DD UNIT*DISK,DCB 3 (RECFM 3 VBS, LRECL«44 , BLKSIZE*4404 ) , 

// SPACE 3 (CYL, (1,1)) MERGING FOR TRACKING DATA SCHEDULE 

//FT03F001 DD UNIT -DISK, DCB 3 ( RECFM 3 VBS , LRECL»44 , BLKSIZE *4404 , 

// BUFN0»1 ) .DISP 3 (NEW, PASS) , SPACE 3 (CYL, (1,1)), 

// DSN-&&SORTRK SORTED TRACKING DATA SCHEDULE 

//FT05F001 DO DDNAME 3 DAT A5 TRACKING SCHEDULE CARD INPUT 

//•FT06F001 DD SYSOUT 33 , DCB 3 ( RECFM 3 VBA, LRECL 3 137 , BLKSIZE=3990) . 

//• DSN 3 &&PRTOUT SEA PRINTER OUTPUT 

//FT06F001 DD SYSOUT 3 *, 

// DCB- ( RECFM 3 VBA , LRECL 3 137 , BLKSIZE 3 141 ) , 

// SPAC£ 3 (CYL,(3,1),RLSE) 

//•FT08F001 DO DDNAME 3 INPCRD SEA KEYWORD CARD INPUT 

//•INPCRD DD * 

//FT08F001 DD * SEA KEYWORD CARD INPUT 


EPOCHTIM 

3 


-1 

881101 . 

000000.0000 


USERSATO 

1 

11 

1 

6978140. 

0.0001 

28 8 

USERSATO 

1 

12 

1 

0. 

0. 

0. 

SATELITE 

2 

11 

0 

421667S0 

0 . 0004 

5.0 

SATELITE 

2 

12 

0 

358 

0 

0 

SATELITE 

3 

11 

0 

42163592.42 

0 0004 

5 0 

SATELITE 

3 

12 

0 

228. 

0 

0. 

SATELITE 

4 

11 

0 

42163592 42 

0.0004 

5 0 

SATELITE 

4 

12 

0 

113. 

0. 

0. 

STATIONS 

1 

0 

0 

323002 . 867 

2532329. 163 

1441 

EARTH 

15 

IS 




150.0 

SPCPARAM 

1 

0 

0 

0.00272 


2.0 

EBOUTPUT 


-1 

2 

1440.0 

30.0 

30.0 

CLKBIAS 

1 


1 

1000000. 



CLKDRIFT 

1 


1 

200.0 

0.000001 


USERDRAG 

1 


1 

.025 



GRAVCOEF 



-1 

0.25 



/* 







// DD UNIT 

•DISK , OSN=ZBBTF . ZBEXS . GEM 

. DATA(REC0EF2 ) , 

DISP=SHR 

// DD * 







SATSOLPR 

1 


-1 

1 



MEASBIAS 

2 


-1 

1000200010000. 

10.0 


MEASBIAS 

3 


-1 

1000200010000. 

.001 
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MEASBIAS 

2 

-1 

1000300010000 

10 0 



MEASBIAS 

3 

-1 

1000300010000. 

.001 



MEASBIAS 

2 

-1 

1000200040001 . 

10.0 



MEASBIAS 

3 

-1 

1000200040001 . 

.001 



EPHEMERR 

1 


25.0 

23 0 


40.0 

EPHERROR 


99 





CLKACCEL 

1 

-1 

.11574 




COVARANC 

1 






2 . 5D+05 

2. 

5D+05 2 . 5D+05 

1 OD+OO 

1 OD+OO 

1 OD+OO 

NOISECOV 

1 

1 0 






0 O. 0. 1 000Q0D-10 1 . 00000D- 10 1.00000D-10 

/* 

//*FT09F001 DD SYSOUT®* , DCB= ( RECFM=VBA, LRECL=137 , BLKSIZE =3990 ) , 

//* DSN-&&NOMOUT SATELLITE NOMINAL TRAJECTORY OUTPUT 

//FT09F001 DD SYSOUT®*, 

// DCB®(RECFM=VBA,LRECL=137,BLKSIZE=141 ) , 

// SPACE® (CYL, (2,1) , RLSE ) 

/ /*FT10F001 DD SYSOUT=* , DCB= ( RECFM® VBA, LRECL® 137 , BLKSIZE =39 90 ) , 

//« DSN s &&MEASOUT PROCESSED MEASUREMENTS PRINTER OUTPUT 

//FT10F001 DD SYSOUT®*, 

// DCB® ( RECFM=VBA, LRECL=137 , BLKSIZE® 141 ) , 

// SPACE = ( CYL , ( 2 , 1 ) , RLSE ) 

//FT20F001 DD DISP® (NEW, PASS) .SPACE® (CYL, ( 1 . 1 )) . 

// DCB® ( RECFM=VBS , LRECL® 144 , BLKSIZE=2884 ) , UNIT=DISK , 

// DSN®&&FESUM SEQUENTIAL ERROR BUDGET SUMMARY 

//FT21F001 DD DISP® ( NEW , PASS ), SPACE® ( CYL ,( 1 , 1 )) , 

// DCB® ( RECFM=VBS , LRECL® 144 , BLKSIZE=2884 ) , UNIT=DISK , 

// DSN=&&SESUM BACKWARD SMOOTHER ERROR BUDGET SUMMARY 

//FT27F001 DD UNIT=DISK,DCB=(RECFM=VBS,LRECL=196,BLKSIZE=1964, 

// BUFNO® 1 ) , DISP® ( NEW , PASS ) , SPACE® ( CYL , ( 1 , 1 ) ) . 

// DSN-&&AKKHIN VISIBLE TRACKING DATA SCHEDULE 

//FT40F001 DD UNIT-DISK, DCB® ( RECFM»FB , LRECL® 1688 . BLKSIZE® 18568 ) . 

// DISP® ( NEW , DELETE ) , SPACE = ( CYL .(2,1), RLSE ) , 

// OSN-&&SMSTOR SMOOTHER INFORMATION STORAGE FILE 

//GO . SYSUDUMP DD DUMMY 

//*SYSUDUMP DD SYSOUT®® , SPACE® (TRK, 1 ) , 

/ / * DCB ® ( RECFM® VBA , LRECL® 1 37 , BLKSIZE = 1 922 ) 

//GO. DAT AS DD * 

3 861101000000.0000 -1 


100020001 

2 180.0 

1 . 66667 


0.00 

2.7 

95.3333 

132.0 

200.3333 

235.7 

304 . 6666 

339 7 

407 . 6666 

443.7 

510.0 

547.0 

612.3333 

650.0 

715.3333 

753 0 

819.6666 

856.0 

924.3333 

959 7 

1028.3333 

1064.0 

1131.3333 

1167.7 

1233.6666 

1271.0 

1336.3333 

1374.4 

1439 6666 

1440.0 



100020001 

3 180.0 

0.00166667 


0.00 

2.7 

95.3333 

132.0 

200.3333 

235.7 

304 . 6686 

339 7 

407 6566 

443.7 

510.0 

547 0 

612.3333 

650.0 

715.3333 

753 0 

819.6666 

- 856.0 

924.3333 

959 7 

1028.3333 

1064.0 

1131 3333 

1167 7 

1233.6666 

1271 0 

1336 3333 

1374 4 

1439.6666 

1440 0 



100030001 

2 180.0 

1 66667 


63.6666 

95.0 

166.3333 

200 0 

268.6666 

304.4 

371.6666 

407 4 

475 3333 

509.7 

580.0 

612.0 

684.3333 

715.0 

787 6666 

819.4 

890.3333 

924.0 

992.6666 

1028.0 

1095.6666 

1131 0 

1199.6666 

1233.4 

1304.0 

1336.0 

1408.6666 

1439.4 

100030001 

3 180 0 

0.00166667 


63 6666 

95.0 

166.3333 

200 0 

268.6666 

304.4 

371.6666 

407 4 

475 . 3333 

509.7 

580.0 

612.0 

684.3333 

715.0 

787 6666 

819 4 

890.3333 

924.0 

992.6666 

1028.0 

1095.6666 

1131.0 

1199.6666 

1233.4 


A-44 



1304.0 


1336.0 


1408 6666 


1439.4 


1000200040001 2 

180.0 

1 66667 


3 00 

63.4 


132 3333 

166.0 

236.0 

268.4 


340.0 

371.4 

444 0 

475.0 


547.3333 

579.7 

650.3333 

684.0 


753 3333 

787.4 

856.3333 

890.0 


960 0 

992 4 

1064.3333 

1095.4 


1168.0 

1199.4 

1271.3333 

1303.7 


1374.6666 

1408.4 

1000200040001 3 

180 0 

0.00166667 


3.00 

63.4 


132.3333 

166.0 

236.0 

268.4 


340 0 

371 4 

444.0 

475 0 


547 3333 

579.7 

650.3333 

684 0 


753.3333 

787.4 

856.3333 

890.0 


960.0 

992.4 

1064.3333 

1095.4 


1168 0 

1199.4 

1271.3333 

1303.7 


1374.6666 

1408 4 

100020001 

3 

1800.0 

10000.0 


0.00 

1440.0 




100030001 

3 

1800.0 

10000 0 


0.00 

1440 0 




1000200040001 3 

1800.0 

10000.0 


0.00 

1440.0 





/* 

//NTSO EXEC PGM=NOTIFY , COND=EVEN 

// 
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APPENDIX B - UTILITIES FOR GEOPOTENTIAL ERROR MODEL STUDIES 


This appendix describes two utilities adapted from the Se- 
quential Error Analysis (SEA) program for geopotential error 
model studies. Both utilities compute the gravitational 
acceleration errors that result from a given geopotential 
error model. The first utility outputs the magnitudes of 
the acceleration errors as a function of geocentric latitude 
and longitude at a fixed altitude. The second utility out- 
puts the magnitudes as well as the acceleration components 
in the radial, alongtrack, and crosstrack directions as the 
satellite travels in its orbit. These two utilities are 
described in more detail in Sections B.l and B.2. 

B.l GEQPOTENTIAL ACCELERATION ERROR WORLD-MAP UTILITY 

This SEA utility computes the uncertainty in the gravita- 
tional acceleration due to nonspherical effects of the Earth 
at a specified altitude for a given geopotential error 
model. The utility computes and prints out the global dis- 
tribution of gravitational acceleration uncertainties in 
meters per second squared for latitudes between -89 deg and 
+89 deg 1 with 5-deg increments and longitudes varying be- 
tween 0 deg and 360 deg with 10-deg increments. 

The utility also generates a global distribution plot of the 

- 5 

gravitational acceleration uncertainties in milligals (10 m/ 
2 

sec ) rounded to the nearest integer, with longitude as 
abscissa and latitude as ordinate. 


^The geopotential computation subroutine is adapted from 
the corresponding GTDS subroutine, which has a singularity 
at the poles (+90 deg) . The use of +89 deg avoids this pro- 
gram deficiency. 
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B.1.1 INPUT DATA SETUP 


This utility can be executed using the SEA keyword cards 
setup, and no special keyword cards are required. The geo- 
potential error model is provided by the LUMPGEOP keyword 
card followed by the corresponding set of RECOEF cards if 
other than the default model--GEM-9 formal uncertainties--is 
desired. The coefficients of any geopotential error model, 
except the default model, can be scaled using the first real 
field of the LUMPGEOP keyword card. The altitude for which 
geopotential acceleration uncertainties are computed is 
based on the orbital elements provided on the USERSATO card. 

B.1.2 MODIFICATIONS TO SEA PROGRAM 

Modifications are made in subroutines RDKEYS and SPART to 
adapt SEA to function as described in Section B.l. In sub- 
routine RDKEYS, which reads in keyword cards, modifications 
are made so that the first real field (Rl) on the LUMPGEOP 
keyword card can be used as a scaling factor (multiplication 
factor) to the sine and cosine coefficients uncertainties of 
an input geopotential error model. This provides a con- 
venient means to scale down or scale up the coefficients to 
study the effects of an error model. 

Subroutine SPART, which normally computes the gravitational 
accelerations due to a geopotential model, is modified to 
compute the acceleration uncertainties resulting from a geo- 
potential error model instead. A nested DO-loop for lati- 
tudes and longitudes is added so that the subroutine 
computes and outputs the global distribution of gravita- 
tional acceleration uncertainties. These uncertainties are 
calculated based on the orbital radius of the user satellite 
as provided by the orbital elements in the USERSATO card. 

The subroutine also converts the acceleration uncertainties 
from meters per second squared to milligals (rounded to the 
nearest integer) and generates a global distribution plot of 
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gravitational acceleration uncertainties. After printing 
out the above results, the program execution is allowed to 
terminate instead of executing the rest of the normal SEA 
program. 

B.1.3 REQUIRED CODING UPDATES TO SEA PROGRAM 

The following is the listing of Panvalet updates to the 
Panvalet SEA source program stored in the data set 
GCDEV. MVT. SEA. PANLIB. DATA. To create an updated SEA program 
load module, the updates that follow must be compiled and 
linked to the existing SEA load module GCDEV. SEAMVS . LOAD. 
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♦♦UPDATE RDKEYS. 1 .TEMP 
++C 229.229 

DPMAP! I LUMP ) = 1 ODO 
IF ( R1 NE. 0. ODO ) DPMAP( I LUMP ) = R1 
++C 251 ,252 

CSZGdl. 12+1 ) = R1 * DPMAPI I LUMP ) 

CSIG( 31-11, 33-12 ) = R2 * DPMAPI I LUMP ) 

♦♦WRITE WORK, RDKEYS 
♦♦UPDATE SPART , 1 , TEMP 
♦♦C 69,69 

COMMON/CONST / AE , GM, 0UM(3), THETG , DTR 
++C 73 73 

COMMON/FMOOEL/CSIG! 30 , 33 ) , C( 30 . 33 ) , MAXDEG , MAXORD . NMAX , MMAX 

++C 76 

LOGICAL* 1 HORZ( 115) .VERT, SLASH 
INTEGER*2 PLOT( 40 . 40 ) , A1 ( 40 ) , A2 ( 40 ) 

DATA DLAT/5 DO/, DL0N/10 DO/, FACTOR/ 1.D5/ 

DATA HORZ/115*'-'/, VERT/' + V, SLASH/'I'/ 

C 

DO 50 ICOL = 1,112,3 
HORZ(ICOL) = SLASH 
50 CONTINUE 
++D 93,94 

WRITEfS, 10IR 

10 FORMAT ( ' 1 ' ///,' *** ENTER SPART TO GENERATE WORLD MAP ', 

* '- USE GEOP UNCER ' ' S TO CALC ACCELERATIONS ***',///. 

* ' THE SIZE OF THE GEOPOTENTIAL UNCERTAINTY IS SPECIFIED BY 

* 'THE "EARTH" KEYWORD',//. 

* ' THE ORBITAL RADIUS FOR WHICH THIS WORLD MAP IS BEING ', 

* 'CALCULATED IS R = ' D25 15,//////, 

* ' THE C ARRAY CONTAINING THE GEOPOTENTIAL UNCER"S IS.' 

C 

C PRINT OUT THE C ARRAY CONTAINING THE GEOPOTENTIAL UNCER 'S 
C 

CALL OUTCOF ( C , MAXDEG , MAXORD , 1 ) 

C 

INDO = 180 DO/DLAT + 1 . 

IND1 = 360 . DO/DLON + 1 
C 

C ENTER LATITUDE LOOP 
C 

WRITE! 6, 390) 

390 FORMAT! '1 NON-SPHERICAL ACCELERATION AS A FUNCTION', 

* ' OF LATITUDE AND LONGITUDE') 

C 

ALAT = (DLAT ♦ 90 DO ) 

DO 1000 1=1, INDO 

IF ( MOD! 1,2) .NE 0 AND. I GT 1 ) WRITE(6.401) 

401 FORMAT ( ' 1 ' ) 

WRITE! 6 , 400) 

400 FORMAT!/) 

ALAT = ALAT - DLAT 
An I) = ALAT 
ALAT 1 » ALAT 

IF ( ALAT GE. 90. DO ) ALAT1 = ALAT - 1 0 
IF ( ALAT LE -90 DO ) ALAT 1 = ALAT + 1 O 
SINP = OSIN! DTR*ALAT1 ) 

COSP = OCOS (DTR* ALAT 1 ) 

++C 112,112 
C 

C ENTER LONGITUDE LOOP 
C 

ALON = - ( DLON + 0 DO) 

DO 2000 J = 1 , IND1 
ALON = ALON + DLON 
A2(d) = ALON 
ALAM = DTR*ALON 

-+C 205 

PLAMDA = PLAMDA / ( R * COSP) 

++C 209 

PPSI = PPSI/R 
C 

C COMPUTE ACCELERATION 
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c 

ACC = DSQRT ( PS*PS + PLAMDA*PLAMDA + PPSI*PPSI) 

C 

WRITE! 6 ,500) ALAT1 , ALON, ACC 

S00 FORMAT! ' ALAT=',G14 8.' AL0N=',G14 8,' ACC=',G14 8) 

LOAD PLOT ARRAY 
PLOT ( I , d ) = ACC* FACTOR + 0.5 
2000 CONTINUE 
1000 CONTINUE 

******** GENERATE PLOT ******** 

WRITE (6, 700) 

700 FORMAT! '1' ,//,T36, ' UNCERTAINTY IN NONSPHERICAL GRAVITATIONAL', 
* ' ACCELERATION ( MGAL ) ' , / , 2X , ' <LAT> ' 3X,/) 

C 

DO 800 1=1, INDO 

WRITE (6,750) Aid), VERT , ( PLOT! I , J ) , J« 1 , IND1 ) 

750 F0RMAT(3X, I3,3X,A1 ,4013) 

800 CONTINUE 
C 

WRITE ( 6 , 770 ) H0R2 , ( A2 ( I ) . I = 1 . IND 1.2) 

770 FORMAT! /.T10, 115A1,//,2X, '<LONG>' ,2X,I3. 1916 I 
C ********************************** 

IF ( 5 NE 0 ISTOP 
♦♦WRITE WORK , SPART 
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B.2 UTILITY FOR COMPUTING GRAVITATIONAL ACCELERATION UNCER- 
TAINTIES ALONG THE TRAJECTORY OF A USER SATELLITE 


This SEA utility computes the gravitational acceleration 
uncertainties that result from a given geopotential error 
model as a user satellite travels in its orbit. The utility 
computes and prints out a table of radial, crosstrack, and 
alongtrack components; the magnitudes of the gravitational 
acceleration uncertainties; and the latitude and longitude 
of the ground trace of the satellite, for each error budget 
time requested. 

The data are then presented in five separate plots. The 
first four plots show the radial (height), crosstrack, and 
alongtrack components and the magnitudes of the gravita- 
tional acceleration uncertainties, respectively, as a func- 
tion of time from epoch. The last plot shows the magnitudes 
of the gravitational acceleration uncertainties as the 
satellite moves in its orbit as a function of latitude and . 
longitude . 

B.2.1 INPUT DATA SETUP 

This utility can be executed using the SEA keyword cards 
setup, and no special keyword cards are required. The geo- 
potential error model can be set up the same way as de- 
scribed in Section B.1.1. 

B.2. 2 MODIFICATIONS TO SEA PROGRAM 

In subroutine RDKEYS , which reads in keyword cards, modifi- 
cations are made so that the first real field (Rl) on the 
LUMPGEOP keyword card can be used as a scaling factor 
(multiplication factor) to the sine and cosine coefficients 
uncertainties of an input geopotential error model. This 
provides a convenient means to scale down or scale up the 
coefficients to study the effect of an error model. 
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In subroutine FORCE, which calculates force per unit mass 
acting on the satellite for the integration of the nominal 
trajectory, additional codings are added to print out a 
table of time from epoch; latitude; longitude; radial, 
crosstrack, and alongtrack components; and the magnitudes of 
the gravitational acceleration uncertainties whenever cur- 
rent time equals error budget time. This subroutine is 
modified so that, for each error budget time, SPART is 
called twice: the first time for computing the geopotential 

acceleration uncertainties due to the geopotential error 
model; the second time to help in computing the nominal tra- 
jectory of the user satellite. The acceleration uncertain- 
ties computed are stored in separate arrays for later 
plotting purpose. To be able to recognize the error-budget 
time, the calling sequence of FORCE is also modified to in- 
clude the input of the variable EBTIME. 

Subroutine FORCE is called from subroutine INTAG for each of 
the integration steps for the propagation of the satellite 
trajectory. Coding modification is therefore required in 
INTAG whenever FORCE is invoked using the modified calling 
sequence. The calling sequence (argument list) of INTAG is 
also modified with the addition of the variable EBTIME, to 
be able to transfer the error-budget time to FORCE. 

Subroutine SEQUEN, which controls the forward sequential 
filter computation, specifies the variable EBTIME, which is 
passed to subroutine INTAG through its calling sequence. 

Modifications are made to suppress the error-budget and 
standard deviation correlation output reports. In addition, 
a new subroutine, GTPLOT , is added at the end of the filter 
computation to generate five separate plots (Section B.2). 
GTPLOT uses the arrays previously loaded in subroutine FORCE 
to generate the required plots. 
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Minor modifications are also made in subroutine SEA01 to 
suppress the error summary reports and error plots. In ad- 
dition to the nonsphencal gravitational accelerations or 
acceleration uncertainties, the calling sequence of subrou- 
tine SPART is modified to pass the latitude and longitude of 
the user satellite to the calling subroutine, FORCE. 

Subroutine THCL is also modified to output the 3 by 3 trans- 
formation matrix for transforming the gravitational acceler- 
ation uncertainties from inertial coordinates to height, 
crosstrack, and alongtrack coordinates. 

B.2.3 REQUIRED CODING UPDATES TO SEA PROGRAM 

The following is the listing of Panvalet updates to the 
Panvalet SEA source program stored in the data set 
GCDEV.MVT. SEA. PANLIB. DATA. To create an updated SEA program 
load module, the updates that follow must be compiled and 
linked to the existing SEA load module GCDEV.SEAMVS.LOAD. 
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-+UPDATE FORCE, 1, TEMP 
++C 1 , 1 

SUBROUTINE FORCE ( T , ISAT , D2Y , ACCEL , EBTIME ) 

++C 40 

LOGICAL* 1 LTOOUT 
++C 45.45 

COMMON/CONST / AE, GM, DUMMY ( 5 ) . RTD 
COMMON/ FMODEL/C( 30.33), CSIG( 30,33) 

COMMON/GTRACK/GTMIN( 96 ) . GALATf 96 ) , GALAMf 96 ) , GGRAV( 96 , 3 ) , GG( 96 ) , 

* IGCNT 

DIMENSION CTEMP( 30,33) 

DATA OEBTIM / -1 ODO / 

++C 50 

DIMENSION GRAV ( 3 ) , TR (3,3) 

C LTOOUT - TRACKING ORIENTED OUTPUT FLAG 
LTOOUT * .FALSE 

C IS EBTIME A NEW ERROR BUDGET TIME 

IF ( EBTIME NE. OEBTIM ) LTOOUT = TRUE 

++D 72 
C 

C GENERATE THE TRACKING ORIENTED OUTPUT ONLY DURING ERROR 
C BUDGET REQUEST TIMES 

C 

IF ( T .NE EBTIME ) GO TO 550 
IF ( NOT. LTOOUT ) GO TO 550 
C 

C FIND THE ACCELERATION DUE TO GEM UNCERTAINTIES CONTAINES IN CSIG 

C - BUT FIRST STORE THEM IN ARRAY C BECAUSE THIS IS THE ARRAY ACTED 

C UPON BY SUBROUTINE SPART 

C 

DO 510 IROW = 1,30 
DO 515 ICOL =1,33 

CTEMP(IROW.ICOL) = C( IROW, ICOL) 

C( IROW, ICOL) * CSIG( IROW, ICOL) 

515 CONTINUE 

• 510 CONTINUE 

CALL SPART ( X , XDD , ALAT , AL AM ) 

DO 530 IROW = 1,30 
DO 535 ICOL *1,33 

C( IROW, ICOL) = CTEMP( IROW, ICOL) 

535 CONTINUE 

530 CONTINUE 

CALL THCL 1 ( STATE (1,1 SAT ),TR) 

CALL MATMUL(TR, XDO, GRAV, 3,3, 1 ) 

G = DSQRT( GRAV( 1 )*GRAV( 1 )+GRAV( 2 )*GRAV(2 )+GRAV( 3 )*GRAV( 3 ) ) 

ALAT * RTD * ALAT 

ALAM = RTD * ALAM 

IF ( ALAM GE 0 ODO ) GO TO 537 

IMULT = 1.0D0 + DABS (ALAM / 360 DO) 

ALAM = ALAM + DFLOAT( IMULT ) * 360 ODO 
537 TMIN = T / 60 ODO 
C 

C LOAD ARRAYS USED LATER FOR PLOTTING 
C 

IGCNT = IGCNT + 1 
GTMIN( IGCNT) * TMIN 
GALAT( IGCNT) = ALAT 
GALAM( IGCNT) = ALAM 
GGRAV( IGCNT, 1) * GRAV( 1 ) 

GGRAV ( I GCNT , 2 ) = GRAV(2) 

GGR AV ( I GCNT , 3 ) = GRAV(3) 

GG( IGCNT) = G 
C 

WRITE ( 6 , 500 ) TMIN , ALAT , ALAM , GRAV , G 
500 FORMAT ( ' TMIN=',G17 10,' ALAT=',G14 8,' ALAM=',G14 8, 

+ ' GRAV* ' , 3 ( G 1 4 8, IX),' G=',G14 8) 

OEBTIM = EBTIME 

550 CALL SPART (X.XOD, ALAT, ALAM) 

++WRITE WORK, FORCE 
++UPDATE INTAG , 1 , TEMP 
++C 1 , 1 

SUBROUTINE INTAG( T 1 , T2 , PXX , TXU , PXZ , D2Y , YDY , YDY 1 , RK , ACCEL , EBTIME ) 
++C 244,244 
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CALL FORCE ' T , I SAT , D2Y . ACCEL , EBTIME 1 
— C 2SS 295 

CALL FORCE i T , I SAT , C2Y , ACCEL , EBTIME ) 

++WRITE WORK, INTAG 
++UPDATE RDKEYS , 1 ,TEMP 
++C 229,229 

DPMAP( I LUMP ) = 1 ODO 

IF ( R1 NE O ODO ) DPMAP( I LUMP ) = R1 
++C 251,252 

CSIG( 1 1 , 12+1 ) = R1 * DPMAPf I LUMP ) 

CSIG(31-I1, 33-12 ) = R2 * DPMAP(ILUMP) 

++WRITE WORK, RDKEYS 
++UPDATE SEQUEN, 1 ,TEMP 
++C 167 

COMMON/GTRACK/GTMIN( 96 ) , GALAT ( 96 ) , GALAMI 96 ) , GGRAV( 96 , 3 ) , GG( 96 ) , 

* IGCNT 
++C 176,177 
++C 207 

IGCNT = 0 
++C 292,292 

95 CALL INTAG( TIME . TTO , PXX , TXU , PXZ , D2Y , YDY , YDY 1 , RK, ACCEL , EBTIME ) 
++D 343,346 
++D 412,413 
++C 458,459 
C 

C GENERATE GEOPOTENTIAL TRACK ERROR PLOTS 

C 

800 CALL GTPLOT ( EBSTOP ) 

++WRITE WORK. SEQUEN 
++UPDATE SEA01 , 1 , TEMP 
++C 216,216 

C SUPPRESS ERROR SUMMARY REPORT AND PLOT, BECAUSE THIER EXECUTION 

C IN COMBINATION WITH THESE UPDATED, RESULTS IN A S0C1 ABEND 

C 

C CALL OUTRSS 

++C 220.220 
C CALL PLTHLC 

++WRITE WORK . SEA01 
++UPDATE SPART , 1 , TEMP 
++C 1 , 1 

SUBROUTINE SPART ( X, XDD, ALAT, ALAM ) 

++C 93 

ALAT = DARSIN(SINP) 

++WRITE WORK, SPART 
++UPDATE THCL.1.TEMP 
++C 4,4 

SUBROUTINE THCL1 ( X, T ) 

++R 17. , /6 , 6/3 , 3/ 

++D 35,37 
++D 55,61 
++WRITE WORK.THCL 
++INSERT WORK 

SUBROUTINE GTPLOTf EBSTOP ) 

C 

IMPLICIT REAL* 8 ( A-H, 0--Z ) 

C 

COMMON/GTRACK/GTMIN ( 96 ) . GALAT ( 96 ) , G ALAM ( 96 ) , GGRAV ( 96 , 3 ) , GG t 96 ) , 

* IGCNT 
C 

DIMENSION TP LOT ( 53 , 103) ,XLINE( 101 ) ,YLINE(53) 

DIMENSION XLAB (11), SYM( 4 ) 

DIMENSION PLOT ( 40 , 111 ) ,DIGITS(36) 

LOGICAL* 1 HORZ( 113) , VERT, SLASH 
INTEGER*2 A1 ( 40 ) , A2 ( 40 ) 

DATA DLAT/5 DO/, DLON/IO DO/, FACTOR/ 1 D5/ 

DATA HORZ/113*'-'/, VERT/'I'/, SLASH/' '/ 

C 

DATA SYM /'H'. 

* 'C' . 

* 'L' , 

* 'R'/ 

DATA DOT /' '/, 

* OASH /'-'/, 
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* ANILET ! ' I / , 

* BLANK ' / , 

» QVERFL /'*'/ 


DATA 

DIGITS/ 'O' , '1', '2' 

, '3' 

'4 

* 

'5' , '6' , '7' 

'8' 

'9 

* 

'A' ,’ 'B' , 'C' 

'D' 

'E 

* 

'F' , 'G' , 'H' 

'I ' 

' J 

* 

'K' , 'L' , 'M' 

'N' 

'0 

% 

'P' , 'Q' , 'R' 

'S' 

'T 

% 

'U' 'V' 'W' 

'X' 

'Y 

* 

'Z'/ 




DO 5 ICOL = 4,112,6 
HORZ(ICOL) = SLASH 
5 CONTINUE 

HORZ! 1) = VERT 
HORZ( 113) = VERT 
C 

DO 7 1=1,40 
DO 8 J* 1,111 

PLOT ( I , J ) = BLANK 
8 CONTINUE 
7 CONTINUE 

BUILD XLINE AND YLINE 

DO 10 1=1 . 101 
10 XLINE (I ) = DASH 
DO 20 L=1 , 101 , 10 
20 XLINE ( I ) = DOT 
DO 30 1=1,53 
30 YLINE(I) = ANILET 
DO 40 1=2,52,5 
- 40 YLINEd ) = DOT 

INITIALIZE ARRAY TPLOT AND LOAD AXIES 

DO 50 IC0L=2 , 102 
DO 60 IR0W=2 , 52 

TPLOT (I ROW, ICOL) = BLANK 
60 CONTINUE 
50 CONTINUE 

DO 80 1=2,102 

TPLOT ( 1 ,1) = XLINE! 1-1) 

TPLOT (53,1) = XLINE! 1-1 ) 

80 CONTINUE 

DO 90 1=1,53 

TPLOT! I, 1 ) = YLINEd) 

TPLOT ( I , 103) = YLINEd) 

90 CONTINUE 

FIND MAX AND MIN OF THE GEOPOTENTIAL ACCELERATIONS 

AMAX = -1 0D-50 
AMIN = 1 00+50 

DO 200 1=1 , IGCNT 

IF ( GGd) GT. AMAX ) AMAX = GG(I) 

DO 210 ISET= 1 , 3 

IF ( GGRAV! I , ISET) LT AMIN ) AMIN = GGRAV! I , ISET ) 
210 CONTINUE 
200 CONTINUE 

WRITE ( 6 , 4000 ) AMAX . AMIN 

000 FORMAT! '0*»»* AMAX = ' , D25 15,' AMIM = ',025 15./) 

FIND SCALING FOR THE X AXIES 

IEBED = EBSTOP/60 ODO + 5D0 

IEBST =00 

IXSCAL = ( ( DFLOAT! IEBED) - DFLOAT ( IEBST ) )/ 100 ODO ) + 
XSCAL = IXSCAL 
IXMIN = IEBST 

WRITE(6, 4010 )EBSTOP, IEBED, IEBST. XSCAL, IXMIN 
4010 FORMAT ( '0==*** EBSTOP =',D14 7,' IEBED = ',16,' IEBST =', 


9999900 


16 , 
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* ' XSCAL = ' ,014 7,' IXMIN =',I6,/I - 

FIND SCALING FOR THE Y AXIES 

YRANGE = ( AMAX - AMIN ) * 1 0D5 
YSCAL = YRANGE / 50 0D0 
YTEMP = YSCAL * 5 0D0 
ROUND YTEMP UP TO NEAREST 1 MG 

IYTEMP = ( YTEMP * 10 ODO ) + 99999DO 

YTEMP = DFLOAT( IYTEMP) / 10 DO 
YSCAL = YTEMP / S 000 

YMIN = ( ( AMIN * 1.0D5 ) / YRANGE ) * ( YSCAL * 50 ODO ) 

ROUND YMIN DOWN TO NEAREST . 1 MG 
IYTEMP = YMIN * 10 DO 

IF ( YMIN LT 0 ODO ) IYTEMP = ( YMIN * 10 DO ) - 99999DO 

YMIN = DFLOAT( IYTEMP) / 10 DO 
YMAX = YMIN + 50 ODO * YSCAL 
WRITE ( 6 , 4020 ) YRANGE , YSCAL , YMIN , YMAX 
,020 FORMAT ( '0**** YRANGE =',D14 7,' YSCAL =',D14 7,' YMIN =',D14 7. 
*' YMAX =' ,D14 7,/) 

LOOP OVER ALL FOUR PLOTS 

WHERE. ISET = 1, HEIGHT PLOT (H) 

2, CROSS TRACK PLOT (C) 

3, ALONG TRACK PLOT (L) 

4, RSS PLOT (R) 

DO 500 ISET =1.4 

LOAD INDIVIDUAL POINTS IN TO ARRAY TPLOT 


YCOR = YMIN / YSCAL 
WRITE ( 6 , 4040 ) ISET, YCOR 

4040 FORMAT ( '0**** ISET =',13, ' YCOR =',D14 7,/) 

DO 600 1=1, IGCNT 

IXCOL = ( GTMIN(I) / XSCAL ) + 5D0 

IF ( ISET LT 4 ) 

* IYROW = ( (GGRAV( I , ISET) * 1 0D5 ) / YSCAL ) - YCOR + 5 

IF ( ISET EQ. 4 ) 

* IYROW = ( ( GG ( I ) * 1 ODS) / YSCAL ) - YCOR + 5 

WRITE( 6 , 4050 ) IYROW , IXCOL , SYM( ISET ) 

4050 FORMAT ( ' ***« IYROW =',I5,' IXCOL ='.15,' USING SYMBOL-> ' , A1 ) 
TPLOT ( 52-IYROW , IXCOL+2 ) = SYM(ISET) 

600 CONTINUE 

PRINT OUT ARRAY TPLOT ONE ROW AT A TIME 
WRITE (6, 2000) 

2000 FORMAT( '1' ,T10, 'PLOT OF TRACK ACCELERATIONS FROM THE 
* 'GEOPOTENTIAL UNCERTAINTY, IN MGAL ' , // ) 

YVAL = YMAX + 5.0 * YSCAL 
DO 700 IROW = 1,53 

WRITE (6, 2010) (TPLO T ( I ROW , ICOL) , IC0L=1 , 103) 

2010 F0RMAT( ' ',T14,103A1) 

IF ( MOD( IROW-2,5) NE 0 ) GO TO 800 
YVAL = YVAL - 5 0 * YSCAL 
WRITE( 6 , 2050 )YVAL 
2050 FORMAT ( ' + ' , T3 , F 10 . 2 ) 

800 CONTINUE 

LATER ON ONE MAY WISH PLACE HERE THE CODE TO GIVE 
Y-AXIS LABELS 
700 CONTINUE 

IXVAL = IXMIN - 10 * IXSCAL 
DO 900 1=1,11 

IXVAL = IXVAL + 10 * IXSCAL 
XLAB(I) = IXVAL 
900 CONTINUE 

WRITE ( 6 , 2 100 )XLAB 
2100 FORMAT ( ' ',T8,11F10 2,//,T47, 

* 'TIMES FROM EPOCH IN MINUTES') 


BEFORE PRINTING NEXT PLOT, BLANK OUT TPLOT ARRAY 


B-I2 



ooono ooo 


IF < I SET EQ 4 i GO TO 500 
DO 1000 IC0 l=2 , 102 
DO 1010 IR0W=2 , 52 

TPLOT ( IROW, I COL ) = BLANK 
1010 CONTINUE 

1000 CONTINUE 
500 CONTINUE 

LOAD PLOT ARRAY TO GENERATE LONGITUDE LATITUDE PLOT 

DO 1020 11=1,1 GCNT 
I = 19 5D0 - ( GALAT(II) / DLAT ) 

ALAMT = GALAM(II) 

JCOL = ( ALAMT / 3 333300 ) + 3.5D0 
AMGAL = GG(II) * FACTOR +05 
IMGAL = AMGAL 

IF ( IMGAL LE 36 ) GO TO 1018 
PLOT (I, JCOL) = OVERFL 
GO TO 1020 
1018 CONTINUE 

WRITE ( 6 , 4060 ) I I , GTMIN ( I I ) , I I , GALAT (II), ALAMT , AMGAL , I , JCOL 
4060 FORMAT (' **** GTMIN ( ' , 12 , ' ) = ' , F7 . 2 , ' GALAT (', 12 , 

* D14 . 7 , ' ALAMT = ' , D 1 4 7,' AMGAL= ' , F6 2,' I=',I3,' JC0L=',I3) 

PLOT ( I , JCOL ) = DIGITS ( IMGAL + 1 ) 

1020 CONTINUE 

GENERATE LATITUDE/LONGITUDE PLOT 

INDO = 180 ODO / DLAT + 1 
IND1 = 360 ODO / DLON + 1 
ILONGV = 0 
ILATV = 90 
DO 1025 1=1,40 
Aid) « ILATV 
ILATV = ILATV - 5 
A2( I ) = ILONGV 
ILONGV = ILONGV + 20 
1025 CONTINUE 

WRITE(6, 1030) HORZ 

1030 FORMAT ( ' 1 ' , // , T3S , ' UNCERTAINTY IN NONSPHERICAL GRAVITATIONAL', 

* ' ACCELERATION (MGAL ) ' , / , 2X, '<LAT>' , 3X , / , 

* 9X.115A1) 

C 

DO 1040 1=1, INDO 

WRITE ( 6 , 1050) A1(I ), VERT, (PLOT(I.J) ,J=1 . Ill) , VERT 
1050 FORMAT (3X,I3,3X,A1,111A1,A1) 

1040 CONTINUE 
C 

WRITE ( 6 , 1060 ) HORZ , ( A2 (I ) , I = 1 , 1 9 I 
1060 FORMAT ( ' ' , T10 , 1 13A1 , // , 10X . 13 , 1816 , // , 60X , 

* ' <LONG>' ) 

RETURN 

END 
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